16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Jochen  Marschall 

a.  REPORT 

U 

b.  ABSTRACT 

UU 

c.  THIS  PAGE 

SAR 

133 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

650-859-2667 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


Final  Report  •  November  201 1 


Testing  and  Modeling  Ultra-High  Temperature  Ceramic 
(UHTC)  Materials  for  Hypersonic  Flight 


SRI  Project  Number:  PI  8570 
AFOSR  Contract  No.  FA9550-08-C-0049 
MP  011-078 


Prepared  by: 

Dr.  Jochen  Marschall 
Molecular  Physics  Laboratory 
SRI  International 


Prepared  for: 

Dr.  Ali  Sayir,  Program  Manager 
High  Temperature  Aerospace  Materials 

Directorate  of  Aerospace,  Chemistry  and  Materials  Science  Air  Force  Office  of  Scientific  Research 
875  North  Randolph  Street 
Suite  325,  Room  3026 
Arlington,  Virginia  22203 


llntornationy 


333  Ravenswood  Avenue  •  Menlo  Park,  California  94025-3493  •  650.859.2000  •  www.sri.com 


INTRODUCTION 


This  report  summarizes  the  activities  and  results  of  experimental  and  modeling  work 
performed  for  the  Air  Force  Office  of  Scientific  Research  (AFOSR)  under  Contract  FA9550-08- 
C-0049,  “Testing  and  Modeling  Ultra-High  Temperature  Ceramic  (UHTC)  Materials  for 
Hypersonic  Flight.”  The  period  of  performance  was  08/01/2008  to  1 1/30/2011.  This  project 
focused  on  two  main  research  areas:  i)  thermal  and  electrical  transport  properties;  and  ii) 
oxidation/volatilization  behavior  in  simulated  re-entry  environments. 

The  first  area  involved  experimental  studies  of  the  thermal  and  electrical  properties  of 
UHTC  materials  as  a  function  of  temperature  and  the  analysis  of  these  results  in  terms  of 
effective  conductivity  and  Wiedemann-Franz  models.  Thermal  diffusivity  was  measured  using  a 
photothermal  radiometry  technique,  and  the  electrical  resistivity  and  Hall  coefficient  were 
measured  using  a  van  der  Pauw  geometry.  Studies  were  performed  on  various  ZrB2-  and  HfB2- 
based  UHTC  materials  manufactured  using  conventional  and  spark  plasma  sintering. 

The  second  area  involved  experiments  under  simulated  re-entry  heating  conditions 
conducted  at  the  1 .2  MW  Plasmatron  facility  of  the  von  Karman  Institute  for  Fluid  Dynamics  in 
Rhode-Saint-Genese,  Belgium.  ZrB2-based  UHTCs  manufactured  by  pressureless  sintering 
techniques  at  the  Missouri  University  of  Science  and  Technology  were  tested  over  a  range  of 
surface  temperatures  exceeding  2000  °C.  A  temperature-jump  phenomenon  was  observed  at 
threshold  heating  levels,  at  which  the  surface  temperature  spontaneously  increased  by  hundreds 
of  degrees  under  constant  free-stream  test  conditions.  This  temperature  jump  is  thought  to  be 
related  to  an  abrupt  transition  in  surface  chemistry  at  the  gas-surface  interface. 

In  addition,  research  has  been  performed  in  the  areas  of  high-velocity  impact  resistance, 
O-atom  oxidation  of  UHTC  materials,  material  property  requirements  for  UHTC  component 
design,  and  in  situ  optical  diagnostics  for  re-entry  heating  test  characterization  and 
documentation.  The  results  of  this  research  have  led  to  5  presentations  at  various  conferences 
and  meetings,  and  have  been  written  up  in  6  technical  publications  (5  published  and 
1  submitted).  Three  further  manuscripts  are  under  preparation. 


SUMMARY  OF  EFFORTS  AND  ACCOMPLISHMENTS 

Research  efforts  and  accomplishments  are  summarized  in  abstract  form  below. 
Comprehensive  details  are  found  in  the  attached  Appendices  (A1-A4),  which  incorporate  edited 
versions  of  published  or  submitted  papers. 


Al)  Temperature  Jump  Phenomenon  during  Plasmatron  Testing  of  ZrB2-SiC  Ultra-High 
Temperature  Ceramics  (submitted  to  the  AIAA  Journal  of  Thermophysics  and  Heat  Transfer, 
August  201 1 ) 

The  performance  of  ZrB2-SiC  ultra  high  temperature  ceramic  composites  in  subsonic 
high-enthalpy  dissociated  air  flows  was  investigated  in  the  1 .2  MW  Plasmatron  facility  at  the 
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von  Karman  Institute  for  Fluid  Dynamics  (VKI).  Samples  manufactured  by  pressureless  sintering 
at  the  Missouri  University  of  Science  and  Technology  (Missouri  S&T)  were  sequentially 
exposed  to  low  heat  flux  and  high  heat  flux  conditions  during  a  single  test  run  by  adjusting  the 
Plasmatron  power.  During  some  high  heat  flux  conditions,  we  observed  a  rapid  and  spontaneous 
rise  in  sample  surface  temperature  under  constant  freestream  conditions.  We  associate  this 
temperature  jump  with  a  transition  in  surface  chemistry  that  leads  to  the  loss  of  protective  silica 
glass  and  substantially  increases  the  chemical  component  of  heat  flux  delivered  to  the  surface. 
These  tests  provide  the  first  set  of  experimental  observations  that  characterize  a  temperature 
jump  phenomenon  in  a  diboride/silica  former  UHTC  composite. 


A2)  Thermal  and  Electrical  Transport  Properties  of  Spark  Plasma  Sintered  HfB2  and  ZrB2 
Ceramics  (The  Journal  of  the  American  Ceramic  Society,  Vol.  94(8),  2011,  pp.  2562-2570.) 

The  thermal  and  electrical  transport  properties  of  various  spark-plasma-sintered  HfB2- 
and  ZrB2-based  polycrystalline  ceramics  were  investigated  experimentally  over  the  300-700  K 
temperature  range.  Measurements  of  thermal  diffusivity,  electrical  resistivity,  and  Hall 
coefficient  are  reported,  as  well  as  the  derived  properties  of  thermal  conductivity,  charge  carrier 
density,  and  charge  carrier  mobility.  Hall  coefficients  were  negative,  confirming  electrons  as  the 
dominant  charge  carrier,  with  carrier  densities  and  mobilities  in  the  3-5  x  10  cm'  and 
100-250  cm2  V"1  s'1  ranges,  respectively.  Electrical  resistivities  were  lower,  and  temperature 
coefficients  of  resistivity  were  higher  than  those  typically  reported  for  HfB2  and  ZrB2  materials 
manufactured  by  conventional  hot  pressing.  A  Wiedemann-Franz  analysis  confirms  the 
dominance  of  electronic  contributions  to  heat  transport.  The  thermal  conductivity  was  found  to 
decrease  with  increasing  temperature  for  all  materials.  Results  are  discussed  in  terms  of  sample 
morphology  and  compared  to  data  previously  reported  in  the  scientific  literature. 


A3)  High-Enthalpy  Test  Environments,  Flow  Modeling  and  In  Situ  Diagnostics  for 
Characterizing  Ultra  High  Temperature  Ceramics  (The  Journal  of  the  European  Ceramic 
Society,  Vol.  30,  2010,  pp.  2323-2336.) 

Ultra-high  temperature  ceramic  materials  and  composites  under  development  as  nose-tip 
and  wing  leading  edge  components  for  hypersonic  flight  vehicles  must  operate  in  extreme 
aerothermal  heating  environments.  The  performance  of  UHTCs  for  this  application  is  ultimately 
evaluated  using  high-enthalpy,  long-duration  flow  facilities  that  simulate  the  reactive  gas 
environment  encountered  in  hypersonic  flight.  In  this  paper,  we  describe  the  test  environments 
generated  by  two  types  of  these  ground  test  facilities — subsonic  inductively-coupled  plasma 
tunnels  and  supersonic  arc -jet  tunnels — and  discuss  the  important  roles  of  computational  fluid 
dynamics  modeling  and  in  situ  optical  diagnostics  for  interpreting  test  results  from  a  materials 
science  perspective. 
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A4)  Material  Property  Requirements  for  Analysis  and  Design  of  UHTC  Components  in 
Hypersonic  Applications  ( The  Journal  of  the  European  Ceramic  Society,  Vol.  30,  2010,  pp. 
2239-2251.) 

Analytical  modeling  of  thermal  and  mechanical  response  is  a  fundamental  step  in  the 
design  process  for  UHTC  components,  such  as  nose  tips  and  wing  leading  edges  for  hypersonic 
applications.  The  purpose  of  the  analyses  is  to  understand  the  response  of  test  articles  to  high- 
enthalpy  flows  in  ground  tests  and  to  predict  component  performance  in  particular  flight 
environments.  Performing  these  analyses  and  evaluating  the  results  require  comprehensive  and 
accurate  physical,  thermal,  and  mechanical  properties.  In  this  paper,  we  explain  the  nature  of  the 
analyses,  highlight  the  essential  material  properties  that  are  required  and  why  they  are  important, 
and  describe  the  impact  of  property  accuracy  and  uncertainty  on  the  design  process. 
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APPENDIX  1 


Temperature  Jump  Phenomenon  during  Plasmatron  Testing  of 
ZrB2-SiC  Ultra-High  Temperature  Ceramics* 


Jochen  Marschall  and  Dusan  A.  Pejakovic 
SRI  International,  Menlo  Park,  California  94025,  USA 

William  G.  Fahrenholtz  and  Greg  E.  Hilmas 
Missouri  University  of  Science  and  Technology,  Rolla,  Missouri  65409,  USA 

Francesco  Panerai  and  Olivier  Chazot 

von  Karman  Institute  for  Fluid  Dynamics,  Rhode-St-Genese,  1640,  Belgium 


Abstract 

The  performance  of  ZrB2-SiC  ultra-high  temperature  ceramic  composites  in  subsonic  high- 
enthalpy  dissociated  air  flows  was  investigated  in  the  1 .2  MW  Plasmatron  facility  at  the  von 
Karman  Institute  for  Fluid  Dynamics.  Samples  were  sequentially  exposed  to  low  heat  flux  and 
high  heat  flux  conditions  during  a  single  test  run.  Under  certain  high  heat  flux  conditions,  a 
rapid  and  spontaneous  rise  in  sample  surface  temperature  was  observed  under  constant 
freestream  conditions.  The  temperature  jump  seems  associated  with  a  transition  in  surface 
chemistry  that  leads  to  the  loss  of  protective  silica  glass  and  substantially  increases  the  chemical 
component  of  heat  flux  delivered  to  the  surface. 


Introduction 

Ultra-high  temperature  ceramic  (UHTC)  composites  of  hafnium  diboride  (HfB2)  and 
zirconium  diboride  (ZrB2)  with  a  silica  former,  most  commonly  silicon  carbide  (SiC),  have  been 
studied  extensively  over  the  last  decade  as  materials  for  leading  edge  and  control  surface 
components  on  hypersonic  vehicles.1'3  Such  components  experience  extreme  aerothermal 
heating  in  chemically  aggressive,  partially  dissociated  air  environments.  Promising  aspects  of 
diboride-based  UHTC  composites  include  the  very  high  melting  points  of  HfB2  and  ZrB2  and 
their  refractory  oxides  hafnia  (Hf02)  and  zirconia  (Zr02),  as  well  as  the  high  thermal 
conductivities  of  HfB2  and  ZrB2  that  enables  efficient  heat  conduction  away  from  stagnation 
point  regions.4  Because  it  is  lighter  and  less  expensive  than  HfB2,  ZrB2  has  some  advantages  as 
an  aerospace  material. 

The  oxidation  of  ZrB2  produces  both  zirconia  and  boron  oxide  (B2O3).  Significant  oxidation 
of  ZrB2  in  atmospheric  air  begins  at  about  1050  K.  The  softening  temperature  for  amorphous 


This  is  an  edited  version  of  a  manuscript  submitted  to  the  AIAA  Journal  of  Thermophysics  and 
Heat  Transfer,  August  2011. 
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B2O3  is  in  the  range  of  830  K  to  900  K5;  below  about  1500  K  the  oxide  scale  consists  of  a  porous 
ZrC>2  network  filled  with  liquid  B2O3  that  acts  as  an  effective  oxygen  diffusion  barrier.6"7 

O 

However,  the  vapor  pressure  of  B2O3  increases  rapidly  with  temperature,  resulting  in  rapid  loss 
of  B2O3  above  1500  K.  The  residual  porous  zirconia  scale  provides  little  resistance  to  inward 
oxygen  transport  and  further  oxidation,9"10  making  the  oxidation  resistance  of  pure  ZrB2 
insufficient  for  high-temperature  hypersonic  vehicle  applications. 

The  addition  of  a  silica  former  to  ZrB2  improves  its  oxidation  resistance.11"14  Compositions 
containing  from  10  to  30%  (by  volume)  SiC  have  generally  been  found  to  be  optimal  in  this 
regard.  The  virgin  ZrB2-SiC  surfaces  oxidize  through  parallel  reactions  that  generate  ZrC>2,  B2O3 
and  SiC>2.  Liquid  B2O3  mixes  with  amorphous  SiC>2  to  form  a  borosilicate  glass  that  seals  the 
ZrC>2  scale.15  With  increasing  temperature,  boron  oxide  evaporates  preferentially  from  the 
surface,  and  the  scale  becomes  silica  rich.  As  the  silica-rich  oxide  layer  thickens,  it  limits  the 
inward  diffusion  of  oxygen  to  the  virgin  material  below,  slowing  further  oxidation.16  The 
oxidation  rates  of  ZrB2-SiC  materials  up  to  -1900  K  at  atmospheric  pressure  are  generally 
proportional  to  the  square  root  of  time  (i.e.,  parabolic  kinetics)  consistent  with  diffusion-limited 
oxidation.4 

The  melting/softening  temperature  of  silica  is  much  higher,  and  the  vapor  pressure  is  much 
lower  than  that  of  B2O3.  Nevertheless,  as  temperatures  exceed  -1900  K,  silica  volatilization  too 
becomes  significant,  again  leaving  an  outer  oxide  scale  of  porous  ZrC>2.  In  this  regime,  oxidation 
proceeds  by  rapid,  linear  kinetics.17  Since  one  goal  of  UHTC  research  is  to  push  the  operating 
temperatures  of  leading  edge  and  control  surface  components  above  2273  K  (2000°C),  it  is 
widely  recognized  that  reliance  on  silica  formers  alone  for  oxidation  protection  is  not  feasible.18 
A  number  of  efforts  are  currently  underway  to  enhance  the  oxidation  protection  properties  of  the 
ZrC>2  scale  through  non-Si-containing  additives.19"21 

A  related  issue  with  SiC-based  high  temperature  composites  is  the  transition  between  passive 
and  active  SiC  oxidation.22  The  transition  occurs  at  threshold  combinations  of  temperature  and 
oxidant  pressure  at  which  the  thermodynamically-favored  oxidation  reaction  switches  from  one 
that  produces  condensed  silica  to  one  that  produces  gaseous  silicon  products  or  vice  versa. 
Active-to-passive  (AP)  and  passive -to-active  (PA)  oxidation  transitions  have  been  studied 
extensively  for  SiC,22"28  as  well  as  other  silica  formers  including  silicon23,29  and  silicon 
nitride.  ’  Most  experiments  were  performed  in  furnace  environments  in  which  temperature, 
pressure,  flow  speed,  and  gas  composition  are  easily  controlled  and  thermochemical  equilibrium 
prevails.  Very  limited  experimental  results  in  partially  dissociated  air  and  oxygen  seem  to 
indicate  that  the  boundary  between  passive  and  active  SiC  oxidation  moves  to  lower  oxygen 
pressures  at  a  given  temperature,  "  but  results  depend  on  the  type  of  SiC  used  in  the 
experiment.32 

The  PA  transition  is  clearly  a  concern  for  any  high-temperature  thermal  protection  system 
(TPS)  containing  a  silica  former  and  relying  on  silica  for  oxidation  protection  (and,  hence,  its 
stability  and  reusability).  Passing  through  the  PA  boundary  means  accelerated  mass  loss,  a 
changing  outer  mold  line,  and  possible  component  failure.  Because  high-temperature  TPS 
systems  will  encounter  partially  dissociated  air  in  chemical  non-equilibrium  with  the  TPS 
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surface,  the  PA  transition  has  been  studied  for  a  variety  of  carbon/silicon  carbide  (C/SiC) 
ceramic  matrix  composites  in  induction-coupled  plasma  (ICP)  '  and  arc-jet  facilities.  '  ’ 
During  these  tests,  spontaneous  jumps  in  surface  temperature  of  several  hundred  degrees  Kelvin 
have  been  observed.  These  have  been  associated  with  passing  through  the  PA  boundary  of  SiC, 
although  neither  the  kinetic  details  triggering  the  PA  transition  in  dissociated  oxygen  and  air,  nor 
the  sources  of  the  additional  heat  flux  necessary  to  explain  the  jump,  are  fully  understood. 

on 

Recently  Glass  suggested  that  Si-containing  UHTC  materials  should  also  exhibit  such 
temperature  jumps,  and  cited  as  a  possible  example  data  collected  during  arc-jet  tests  at  NASA 
Ames  Research  Center  for  the  SHARP  B-l  program  as  reported  by  Kolodziej  et  al.40  In  one 
instance,  as  a  ZrB2-SiC  sample  was  heating  up  in  the  arc-jet  stream,  its  surface  temperature 
suddenly  increased  sharply  from  3800°F  to  4600°F  (2360  K  to  2810  K)  over  the  course  of  about 
15  seconds.  This  temperature  jump  was  rationalized  as  resulting  from  the  lower  emittance  and 
lower  thermal  conductivity  of  a  rapidly  forming  oxide  layer.  Since  this  report,  the  results  of 
numerous  arc-jet  and  ICP  test  campaigns  on  UHTC  composites  containing  a  diboride  and  a  silica 
former  have  been  published,  but,  to  our  knowledge,  none  have  clearly  identified  a  similar 
spontaneous  temperature  jump.41"47 

This  paper  reports  test  results  that  demonstrate  a  temperature  jump  phenomenon  when  ZrB2- 
SiC  composites  are  exposed  to  certain  high-temperature,  low-pressure  dissociated  air 
environments.  The  tests  were  performed  in  the  1 .2  MW  Plasmatron  facility  at  the  von  Karman 
Institute  for  Fluid  Dynamics  (VKI)  on  specimens  manufactured  by  pressureless  sintering  at  the 
Missouri  University  of  Science  and  Technology  (Missouri  S&T).  Specimens  were  sequentially 
exposed  to  low  heat  flux  and  high  heat  flux  conditions  during  a  single  test  run  by  adjusting  the 
Plasmatron  power.  During  some  high  heat  flux  conditions,  we  observed  a  rapid  and  spontaneous 
rise  in  sample  surface  temperature  under  constant  ffeestream  conditions.  We  associate  this 
temperature  jump  with  a  transition  in  surface  chemistry  that  leads  to  the  loss  of  protective  silica 
glass  and  substantially  increases  the  chemical  component  of  heat  flux  delivered  to  the  surface. 
These  tests  provide  the  first  set  of  experimental  observations  that  characterize  the  temperature 
jump  phenomenon  in  a  diboride/silica  former  UHTC  composite. 


Materials  and  Characterization 
A.  Test  Specimen  Fabrication 

Composites  were  prepared  with  the  following  target  formulations:  1)  ZrB2  with  30  volume 
percent  SiC;  2)  (Zr,W)B2  with  4  atomic  percent  tungsten  substitution  for  zirconium;  and  3) 
(Zr,W)B2  with  4  atomic  percent  tungsten  substitution  and  30  volume  percent  SiC.  These  three 
compositions  are  referenced  as  ZrB2-30SiC,  (Zr,4W)B2,  and  (Zr,4W)B2-30SiC  throughout  this 
paper.  Test  specimens  of  these  compositions  are  identified  by  1.x,  2.x,  and  3.x,  respectively, 
where  “x”  is  the  sample  number. 

All  composites  were  made  from  commercially  available  ZrB2  powder  (Grade  B,  H.  C.  Starck, 
Germany),  which  had  a  starting  particle  size  of  ~2  pm.  ZrB2-30SiC  and  (Zr,4W)B2-30SiC 
materials  also  used  SiC  powder  (Grade  UF-10,  H.C.  Starck,  Germany),  which  had  a  starting 
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particle  size  of  ~1  pm;  (Zr,4W)B2  and  (Zr,4W)B2-30SiC  materials  also  used  tungsten  carbide 
(WC)  powder  (74R-0601,  Inframat,  Manchester,  Connecticut),  which  had  a  starting  particle  size 
<  1  pm.  Particle  sizes  are  as  reported  by  the  manufacturer.  Superadditions  of  a  combination  of 
2  wt%  B4C  (Grade  HS,  H.C.  Starck,  Germany)  and  1  wt%  carbon  derived  by  decomposition  of 
phenolic  resin  (GP  2074,  Georgia  Pacific,  Atlanta,  Georgia)  were  used  as  sintering  aids  for  all 
three  compositions. 

Powder  batches  were  produced  by  dispersing  ZrB2  with  SiC,  WC,  and/or  B4C  in  methyl  ethyl 
ketone  using  a  dispersant  (BYK  1 10,  BYK  Chemie  Co,  Willingford,  Connecticut).  Batches  were 
ball  milled  for  24  hours  using  WC  media  in  a  high-density  polyethylene  bottle.  After  24  hours  of 
milling,  1  wt%  polypropylene  carbonate  binder  (QP AC-40,  Empower  Materials,  Newark, 
Delaware)  and  phenolic  resin  were  added,  and  the  slurries  were  ball  milled  for  an  additional 
24  hours. 

After  milling,  the  slurries  were  dried,  ground,  and  sieved  (-50  mesh)  to  prepare  for 
compaction.  Powders  were  pressed  into  cylinders  approximately  40  mm  in  diameter  and  20  mm 
tall  using  a  uniaxial  compression  pressure  of  18.6  MPa.  Next,  the  billets  were  cold  isostatically 
pressed  at  3 15  MPa.  After  compaction,  the  billets  underwent  binder  burnout  in  an  atmosphere  of 
Ar-5%  H2,  and  following  binder  burnout,  they  were  machined  using  a  WC  tool  on  a  standard 
metal  working  lathe  to  produce  near  net-shaped  parts  approximating  the  desired  “mushroom¬ 
shaped”  test  specimens  (Fig.  la). 

The  near  net  shape  specimens  were  then  sintered  according  to  the  following  schedule. 
Specimens  were  first  heated  at  10  °C  min'1  from  room  temperature  to  1250  °C  under  mild 
vacuum  (-20  Pa).  After  holding  at  1250  °C  for  1  hour,  specimens  were  heated  at  10  °C  min1  to 
1450°C  (ZrB2-30SiC  and  (Zr,4W)B2-30SiC)  or  1500°C  ((Zr,4W)B2)  and  held  for  1  hour.  After 
the  second  hold,  the  atmosphere  was  switched  to  flowing  argon  at  a  nominal  pressure  of  ~105  Pa. 
Specimens  were  heated  at  20  °C  min'1  to  the  densification  temperature  and  held  there  for 
2  hours.  The  densification  temperatures  were  2000°C  for  ZrB2-30SiC  and  (Zr,4W)B2-30SiC 
materials  and  2075°C  for  (Zr,4W)B2.  After  the  densification  hold,  specimens  were  cooled  to 
room  temperature  at  the  natural  furnace  cooling  rate  (~20  °C  min"1  on  average). 
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Fig.  1  (a)  UHTC  sample  geometry;  (b)  stagnation  point  test  configuration  with  sample, 
specimen  holder,  and  sting  arm. 


The  sintered  ceramic  specimens  were  fine -machined  by  a  commercial  vendor  using  diamond 
grinding  (Bomas  Machine  Specialties,  Somerville,  Massachusetts)  in  an  attempt  to  reach  the 
target  dimensions  shown  in  Fig.  la.  While  most  dimensions  could  be  achieved  with  good 
accuracy,  it  was  not  possible  to  maintain  the  7  mm  thickness  of  the  “mushroom  caps”  for  the 
majority  of  the  specimens,  resulting  in  thinner  caps  (by  up  to  2  mm)  and  larger  flat  face 
diameters  (by  up  to  9  mm  larger). 

B.  Specimen  Characterization 

The  densities  of  the  machined  UHTC  specimens  were  measured  by  the  Archimedes  method 
using  a  Mettler  Toledo  XP105  analytical  balance  with  ±  0.01  mg  accuracy.  The  average 
densities  for  the  ZrB2-30SiC,  (Zr,4W)B2,  and  (Zr,4W)B2-30SiC  formulations  were  found  to  be 
4.88  ±  0.15  g  cm'3,  5.72  ±  0.28  g  cm'3,  and  4.75  ±  0.21  g  cm'3,  respectively.  Specimen  mass 
changes  were  determined  from  pre-  and  post-plasma  exposure  mass  measurements  at  VKI  using 
an  Adventurer  SL  precision  balance  with  ±  0. 1  mg  accuracy. 

Optical  microscopy  of  post-test  specimen  surfaces  was  performed  using  a  KEYENCE  VHX- 
600  Digital  Microscope  (Itasca,  Illinois).  Post-test  specimen  surface  composition  was  assessed 
with  energy  dispersive  X-ray  spectroscopy  (EDS)  instrumentation  (ED AX,  Mahwah,  New 
Jersey)  integrated  into  a  field-emission  scanning  electron  microscope  (JEOL  6100,  JEOL,  Ltd. 
Tokyo,  Japan).  Crystal  structure  was  determined  by  X-ray  diffraction  (XRD)  using  a  Phillips 
apparatus  (New  York,  New  York)  with  a  CuXa  source  and  a  Ge  monochromator.  Elemental 
analysis  of  the  oxide  scale  was  performed  by  a  commercial  analytic  laboratory  (ATI  Wah  Chang, 
Albany,  Oregon). 
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Plasmatron  Testing  Setup  and  Procedures 

Air  plasma  exposure  experiments  were  performed  in  the  1 .2  MW  Plasmatron  facility  at 
VKI.48'49  A  solid-state  power  supply  (400  kHz,  1 .2  MW,  2  kV)  was  used  to  generate  a  high 
enthalpy  subsonic  gas  flow  by  inductive  coupling  in  a  160  mm  diameter  plasma  torch.  The 
energetic  gas  stream  was  directed  into  a  large  vacuum  chamber,  where  test  specimens  and 
diagnostic  probes  for  the  dynamic  pressure,  Pdyn ,  and  the  cold- wall  heat  flux,  qcw ,  could  be 

swung  into  and  out  of  the  flow  on  three  separate  water-cooled  sting  arms.  The  distance  between 
the  torch  exit  and  the  specimen  or  probe  face  was  44.5  cm. 

Test  specimens  were  mounted  in  graphite  holders  attached  to  a  water-cooled  sting- arm,  as 
shown  in  Fig.  lb.  Specimens  were  held  in  place  with  pins  that  captured  the  grooves  on  the 
specimen  stems  (Fig.  la.)  In  this  campaign,  the  Plasmatron  was  always  operated  with  an  air 
mass  flow  rate  of  16  g  s_1  and  a  static  chamber  pressure  of  10  kPa.  The  freestream  enthalpy  was 
varied  over  the  range  of  about  10  to  28  MJ  kg'1  by  adjusting  the  Plasmatron  power  from  170  to 
390  kW. 

The  cold- wall  stagnation  point  heat  flux  was  measured  using  a  water-cooled  copper 
calorimeter  installed  flush  with  the  surface  of  a  water-cooled  copper  probe.  Both  the  calorimeter 
and  probe  had  polished  surfaces.  Heat  flux  values  were  determined  from  the  temperature  change 
and  mass  flow  rate  of  the  water  used  to  cool  the  calorimeter.  Two  heat  flux  probes  of  different 
size  and  external  shape  were  used  during  the  course  of  these  experiments.  One  probe  was 
geometrically  similar  to  the  test  specimen  (flat-faced  cylinder,  30  mm  body  diameter  and  7  mm 
comer  radius)  so  that  the  flow  environment  around  the  specimen  was  closely  reproduced. 
Unfortunately,  this  probe  was  damaged  early  in  the  campaign  and  a  second  larger  probe  having 
the  standard  European  Space  Agency  (ESA)  stagnation  point  test  specimen  geometry  was  used 
instead  (flat-faced  cylinder,  50  mm  body  diameter,  1 1  mm  corner  radius). 

Cold- wall  heat  flux  measurements  serve  two  purposes:  i)  they  are  used  to  characterize  and 
reproduce  heating  conditions  between  different  test  runs;  and  ii)  they  enable  the  rebuilding  of 
boundary  layer  edge  conditions  and  the  estimation  of  gas  composition  at  the  specimen  surface 
using  computational  fluid  dynamics  (CFD)  codes.  Probe  geometry  is  not  critical  in  the  first 
application,  in  which  only  measurement  consistency  is  important.  However,  probe  geometry  has 
a  large  influence  in  the  second  application  because  aeroconvective  heating  rates  are  strongly 
affected  by  probe  shape  and  dimensions.  In  the  present  instance,  the  large  probe  experiences  a 
lower  cold-wall  heat  flux  than  the  small  probe  for  the  same  Plasmatron  freestream  conditions 
because  of  differences  in  the  boundary  layer  thickness  and  the  velocity  gradient  as  the  flow 
approaches  the  probe  surface.  The  thermal  and  dynamic  boundary  layers  are  thicker  and  the 
velocity  gradient  is  smaller  for  the  large  versus  the  small  probe. 

Additional  tests  comparing  heat  flux  measurements  with  the  two  different  probes  were  carried 
out  after  the  UHTC  test  campaign  concluded.  These  tests  were  performed  by  swinging  the  two 
calorimeter  probes  mounted  on  separate  sting  arms  into  the  flow  consecutively,  at  the  same  static 
pressure,  the  same  air  mass  flow  rate,  and  over  the  same  range  of  Plasmatron  powers  as 
employed  in  the  UHTC  tests.  The  resulting  data  are  shown  in  Fig.  2.  The  linear  correlation 
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established  by  these  measurements,  qcw  =  1.69 qcw  lp  - 17.32 ,  enabled  conversion  of  the  cold-wall 

heat  fluxes  measured  with  the  large  probe  to  values  appropriate  for  the  small  probe  (and  thus  to 
the  UHTC  sample  geometry.)  Both  heat  flux  measurements  should  agree  at  zero  heat  flux,  and  a 
linear  correlation  with  a  zero  intercept  could  also  be  used  to  fit  the  data.  However,  in  the 
absence  of  measurements  at  lower  heat  fluxes,  the  suitability  of  a  linear  fit  over  the  entire  range 
extending  to  zero  cannot  be  confirmed,  and  it  is  safer  to  limit  the  correlation  to  the  range  used  in 
our  tests. 

The  test  chamber  was  fitted  with  several  windows  to  allow  for  optical  diagnostics,  including 
two-color  and  wide-band  infrared  radiometry  of  the  sample  surface,  and  emission  spectroscopy 
of  the  gas  in  front  of  the  sample.  In  addition,  tests  were  also  documented  by  different 
combinations  of  digital  photography  and  videography. 

A  two-color  radiometer  (Marathon  Series  MR1SC,  Raytek  Corp.,  Santa  Cruz,  California)  was 
used  to  record  sample  surface  temperatures  at  a  rate  of  1  Hz.  The  radiometer  measures 
temperatures  in  the  range  1000  to  3000  °C,  using  two  overlapping  infrared  wavelength  bands  at 
0. 8-1.1  and  1. 0-1.1  pm.  Thermal  emission  was  collected  through  a  quartz  window  at  an  angle 
of -35  degrees  relative  to  the  surface  normal,  from  an  elliptical  region  centered  on  the  specimen 
stagnation  point  (major  axes:  ~5.6  mm;  minor  axis  -4.6  mm).  The  radiometer  was  previously 
calibrated  with  the  same  window  in  place  using  a  black-body  radiation  source  (LANDCAL 
R1500T,  Land  Instruments  International,  Dronfield,  United  Kingdom);  temperature 
measurements  are  believed  to  be  accurate  to  ±  10  °C. 

A  broadband  infrared  radiometer  (KT19.XX,  Heitronics  Infrarot  Messtechnik  GmbH, 
Wiesbaden,  Germany)  with  a  temperature  range  of  0  to  3000  °C  was  used  to  collect  thermal 
emission  from  the  sample  surface  over  the  0.6-39  pm  wavelength  range.  Thermal  emission 

measurements  were  recorded  at  a  frequency 
of  1  Hz.  Emission  was  collected  through  a 
KRS5  window  at  -47  degrees  relative  to  the 
surface  normal  from  an  elliptical  region 
centered  on  the  stagnation  point  (major  axes: 
-14.8  mm;  minor  axis  -10.1  mm).  The 
broadband  radiometer  was  calibrated  with  the 
KSR5  window  in  place  using  the  same 
blackbody  source  as  the  two-color  radiometer. 

With  this  geometric  arrangement  of  the 
two  radiometers,  and  with  their  respective 
fields  of  view,  the  emitting  area  seen  by  the 
two-color  radiometer  lies  within  the  area  seen 
by  the  broadband  radiometer,  which  in  turn 
lies  completely  on  the  flat  face  of  the  sample 
away  from  the  rounded  edges. 

Tests  were  filmed  with  a  high-definition 


Fig.  2.  Cold-wall  heat  flux  probe 
measurement  comparison;  data  and 
linear  correlation. 
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digital  dual  camera  (YPC-GH1,  Sanyo  Europe  Ltd.,  Watford,  United  Kingdom).  The  camera 
was  positioned  to  view  the  sample  through  a  quartz  window  ~35  degrees  off  the  surface  normal. 
The  camera  has  a  1920  by  1080  C-MOS  sensor,  a  double  range  focal  lens,  and  a  video  shutter 
speed  that  ranges  from  1/30  to  1/10000  seconds. 

Gas-phase  emission  spectrometry  was  performed  through  a  quartz  window  viewing  the  test 
specimen  and  impinging  gas  flow  from  the  side.  Emission  was  collected  through  a  variable 
aperture,  and  focused  by  a  300  mm  focal  length  converging  quartz  lens  into  a  600  pm  diameter 
optical  fiber  leading  to  a  spectrometer  with  a  200-1100  nm  wavelength  range  (HR-4000,  Ocean 
Optics,  Dunedin,  Florida.)  Spectra  were  recorded  at  1  Hz  with  spectral  resolution  ranging  from 
0.27  nm  to  0.23  nm.  The  optical  axis  of  the  collection  system  was  aligned  parallel  to  the  sample 
surface  at  its  stagnation  point,  with  the  ~0.5  mm  diameter  emission  collection  volume  about 
2  mm  in  front  of  the  surface.  A  relative  intensity  calibration  of  the  system  was  performed  with  a 
deuterium  lamp  in  the  ultraviolet  range  (<  400  nm)  and  with  a  ribbon  tungsten  lamp  in  the  visible 
(>  400  nm).  The  estimated  relative  uncertainty  of  the  total  integrated  emission  of  the  radiative 
signature  is  less  than  10%. 


Computation 

Computational  fluid  dynamics  simulations  of  the  plasma  ffeestream  and  the  boundary  layer 
around  the  specimen  and  holder  were  used  to  rebuild  the  boundary  layer  edge  conditions  and  to 
compute  the  gas  composition  and  the  catalytic  component  of  heating  at  the  specimen  surface. 

The  VKI  Boundary  Layer  Code  solves  the  boundary  layer  equations  for  a  two-dimensional, 
steady,  laminar  flow  of  chemically  reacting  gas  over  a  catalytic  surface,  including  chemical  non¬ 
equilibrium.50'51  The  VKI  ICP  Code  computes  the  Plasmatron  magnetohydrodynamic  field 
(plasma  torch  plus  vacuum  chamber)  under  local  thermodynamic  equilibrium  and  provides 
nondimensional  parameters  that  characterize  the  boundary  layer  over  the  test  article.  '  Both 
codes  use  the  PEGASE  library  to  obtain  thermodynamic  and  transport  properties  of  the  gas 
mixture.54 

The  total  heat  flux  experienced  by  a  material  in  the  Plasmatron  stream  has  both  convective 
and  chemical  components:  qw  =  qconv  +  qchem .  Convective  heating  arises  from  the  large 

temperature  difference  between  the  hot  boundary  layer  gases  and  the  cooler  specimen  surface, 
while  chemical  heating  results  from  exothermic  reactions  on  the  specimen  surface  such  as 
catalytic  recombination  and  oxidation.  The  Boundary  Layer  Code  computes  the  chemical 
component  of  heat  flux  assuming  the  two  independent  catalytic  surface  reactions,  O  +  O  — ►  O2 
and  N  +  N  —*■  N2,  and  restricts  them  to  have  the  same  recombination  efficiency  with  complete 
accommodation  of  the  exothermic  recombination  energy.  The  recombination  efficiency,  y ,  is 

defined  as  the  fraction  of  atom  collisions  with  the  surface  that  result  in  atom  loss.  The  maximum 
chemical  heating  component  is  computed  with  y  =  1  and  the  minimum  ( qchem  =  0 )  with  y  =  0 . 

Under  steady-state  conditions,  the  input  heat  flux  is  balanced  by  radiative  energy  loss  from  the 
surface  to  the  environment  and  by  heat  conduction  into  the  interior:  qw  =  qmd  +  qcond  .  During 

calorimeter  measurements,  water  cooling  keeps  the  surface  temperature  low  making  qmd 
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negligible,  and  the  cold-wall  heat  flux  is  measured:  q cw  =  qcond  .  During  specimen  tests,  the 

surface  reaches  a  high  steady-state  temperature  and  the  hot-wall  heat  flux  is  approximately 
qhw  =  qmd  ,  where  it  is  presumed  that  qcond  is  much  smaller  than  qrad  .  The  radiative  flux  is 

computed  from  the  radiance  measured  by  the  broadband  radiometer,  or  equivalently  from  the 
measured  surface  temperature,  Tw ,  as  qrad  =  soT ,  where  cr  is  the  Stefan-Boltzmann  constant 

and  s  is  the  total  hemispherical  emittance  of  the  surface  derived  from  the  two-color  and 
broadband  radiometer  measurements,  as  detailed  in  Section  V  below. 

The  Boundary  Layer  Code  and  the  ICP  Code  were  used  together  to  rebuild  the  flow 
conditions  at  the  boundary  layer  edge.  First,  the  non-dimensional  boundary  layer  parameters 
were  computed  by  the  ICP  Code  for  each  test  condition  -  mass  flow  rate,  static  pressure,  and 
power  coupled  into  the  gas  -  given  the  geometries  of  the  30  mm  probe  and  the  Plasmatron 
facility.  With  the  assumption  of  local  thermodynamic  equilibrium  at  the  boundary  layer  edge, 
the  computed  non-dimensional  parameters  then  served  as  inputs  to  the  Boundary  Layer  Code.  A 
five-species  model  of  dissociated  air  (O,  O2,  N,  N2,  and  NO)  and  the  gas-phase  reaction  rates  of 
Dunn  and  Kang55  were  used.  The  copper  calorimeter  was  assigned  a  catalytic  recombination 
efficiency  of  0.01  and  a  surface  temperature  of  350  K.  The  Boundary  Layer  Code  was  iterated 
over  the  boundary  layer  edge  temperature  until  the  numerical  and  experimental  cold- wall  heat 
fluxes  matched.  The  final  outputs  of  this  procedure  were  the  five  species  concentrations,  the 
temperature,  Te ,  the  velocity,  Ve ,  and  the  enthalpy,  He  at  the  boundary  layer  edge. 

Once  the  edge  conditions  were  determined,  the  Boundary  Layer  Code  was  used  to  compute  an 
abacus  of  heat  flux  vs.  surface  temperature  for  different  values  of  recombination  efficiency.  For 
each  point  on  the  abacus,  the  Boundary  Layer  Code  concurrently  computes  the  gas-phase 
concentrations  of  O,  O2,  N,  N2,  and  NO  in  contact  with  the  hot  wall  at  the  stagnation  point.  The 
intersection  of  the  measured  surface  temperature  and  hot-wall  heat  flux  on  this  abacus  identified 
the  recombination  efficiency  necessary  to  reproduce  the  experiment,56  and  the  particular  CFD 
solution  used  to  estimate  the  composition  of  the  gas  impinging  on  the  specimen  surface. 


Experimental  Results 

A.  Observation  of  the  Temperature  Jump  Phenomenon 

Figure  3  shows  the  surface  temperature  histories  measured  by  the  two-color  radiometer  during 
Plasmatron  testing  of  the  various  UHTC  specimens.  Tests  were  performed  by  inserting  the 
specimen  into  the  flow  at  a  low  Plasmatron  power  setting,  letting  the  sample  temperature 
equilibrate  for  a  given  time  (typically  5.5  minutes),  then  increasing  the  power  and  holding  at  the 
new  heating  condition  for  an  additional  length  of  time  (typically  2  to  4.5  minutes).  Tests  were 
terminated  by  either  shutting  off  the  torch  or  removing  the  specimen  from  the  flow;  the  former 
method  produces  the  gradual  temperature  decays  seen  in  Fig.  3.  The  dynamic  pressures  and  the 
cold-wall  heat  fluxes  were  measured  at  the  low-power  and  high-power  settings  during  each  test 
prior  to  sample  insertion.  Table  1  lists  the  Plasmatron  powers,  the  measured  cold- wall  heat  flux 
and  dynamic  pressure  (rounded  to  the  nearest  5  Pa)  and  the  derived  boundary  layer  edge 
conditions  for  each  test.  The  Plasmatron  powers  listed  correspond  to  an  approximate  value 
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determined  from  an  image  of  the  full  power  delivered  by  the  high  frequency  generator;  the 
power  transferred  into  the  gas  stream  is  lower. 


0  1 


3456789  10  11 
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Fig.  3.  Surface  temperatures  measured  during  Plasmatron  testing  of  ZrB2-30SiC  (top 
panel),  (Zr,4W)B2-30SiC  (middle  panel),  and  (Zr,4W)B2  (bottom  panel)  UHTC 
specimens.  Open  diamonds  indicate  when  the  Plasmatron  power  was  increased  and  open 
triangles  indicate  the  onset  of  spontaneous  surface  temperature  jumps.  Specimens  that 
underwent  a  temperature  jump  are  identified  in  bold  text  in  the  figure  legends. 
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Table  1  Plasmatron  operating  conditions.* 


Sample* 

Power, 

kW 

QcW  ’ 

W  cm  2 

P 

1  dyn  ? 

Pa 

He, 

kJg  1 

T 

±e  ’ 

K 

ve, 

m  s  1 

1.4 

180 

-f— 

OO 

o 

25 

12.2 

5048 

92 

290 

185' 

60 

17.7 

5656 

159 

1.5 

200 

109 

30 

12.0 

5003 

100 

353 

202* 

75 

18.6 

5731 

180 

1.6 

205 

109 

30 

12.0 

5003 

100 

387 

236 

90 

21.7 

5944 

207 

1.7 

210 

108 

25 

12.2 

5048 

92 

370 

219 

80 

20.3 

5849 

191 

1.8 

212 

109 

30 

12.0 

5003 

100 

292 

185 

60 

17.7 

5656 

159 

1.9 

200 

111 

30 

12.2 

5034 

101 

2.4 

195 

111 

30 

12.2 

5034 

101 

2.6 

198 

111 

30 

12.2 

5034 

101 

330 

245 

90 

22.9 

6015 

210 

2.10 

172 

57 

20 

7.2 

3464 

65 

261 

152 

45 

15.1 

5427 

131 

2.11 

170 

54 

20 

6.9 

3380 

64 

3.2 

195 

109 

30 

12.0 

5003 

100 

310 

194* 

70 

18.0 

5681 

172 

3.3 

224 

111 

30 

12.2 

5034 

101 

339 

202 

75 

18.6 

5731 

180 

3.4 

172 

111 

30 

12.2 

5034 

101 

3.5 

203 

111 

30 

12.2 

5034 

101 

402 

253 

95 

23.5 

6055 

218 

3.9 

211 

109 

30 

12.0 

5003 

100 

295 

185 

60 

17.7 

5656 

159 

*  p 

stat 

=  104  Pa  and 

m  =  16  g  s 

1  for  all  tests; 

*  Estimated  based  on  similar  tests; 

^Samples  that  underwent  a  temperature  jump  are  indicated  in  bold  text. 


Upon  sample  insertion,  the  surface  temperature  rises  and  reaches  a  relatively  steady  value 
determined  by  the  energy  balance  between  the  hot-wall  heat  flux  to  the  specimen,  heat 
conduction  losses  into  the  interior,  and  radiation  to  the  environment.  The  low  power  setting  for 
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the  majority  of  tests  was  adjusted  to  produce  a  cold- wall  heat  flux  of  ~1 10  W  cm'  .  Under  these 
conditions  all  ZrB2-30SiC  and  (Zr,4W)B2-30SiC  specimens  reached  steady  surface  temperatures 
between  1800  and  1900  K.  The  (Zr,4W)B2  material  reached  a  considerably  higher  surface 
temperature  of  about  2130  K  under  these  same  test  conditions  (specimens  2.4  and  2.6  in  Fig.  3.) 

Additional  tests  were  run  on  samples  2.10  and  2.1 1  with  the  low-power  setting  decreased  to 
produce  a  cold- wall  heat  flux  near  55  W  cm'2.  This  lower  heating  condition  resulted  in  surface 
temperatures  similar  to  those  observed  for  the  ZrB2-30SiC  and  (Zr,4W)B2-30SiC  materials  at 
1 10  W  cm'1.  The  lower  surface  temperatures  experienced  by  the  SiC-containing  materials  versus 
the  (Zr,4W)B2  specimens  are  likely  due  to  the  development  of  a  silica-rich  oxide  scale  covering 
the  ZrC>2  formed  during  testing,  which  retards  the  rate  of  oxidation.  Additionally,  silica  surfaces 
are  known  to  have  relatively  low  catalytic  efficiencies  for  oxygen  recombination57  compared  to 
zirconia  surfaces,58  and  so  will  experience  lower  chemical  heating  contributions  to  the  heat  flux 
for  a  given  freestream  condition. 

The  transition  between  the  low-power  and  high-power  settings  during  each  test  run  is 
indicated  in  Fig.  3  by  open  diamonds  on  the  temperature  traces.  An  increase  in  Plasmatron 
power  leads  to  an  increased  cold-wall  heat  flux  (Table  1),  and  is  expected  to  correspondingly 
increase  the  hot-wall  heat  flux  to  the  specimen  and  drive  its  surface  temperature  to  a  higher 
steady-state  value.  Examples  of  this  anticipated  behavior  are  seen  for  samples  1.4,  1.8,  3.2,  and 
3.9,  where  the  high-power  setting  has  been  adjusted  to  produce  cold- wall  heat  flux  of 
~1 85  W  cm'2  and  all  four  specimens  settle  to  a  steady-state  surface  temperature  near  21 10  K. 
Similar  behavior  is  seen  for  sample  2.10  when  the  cold- wall  heat  flux  is  raised  from  57  to  152  W 
cm'2.  Although  the  temperature  measurement  shows  large  fluctuations,  it  varies  around  a  mean 
steady-state  temperature  of  about  2400  K  under  the  higher  power  conditions. 

However,  five  of  the  SiC-containing  specimens  -  1.5,  1.6,  1.7,  3.3  and  3.5  -  showed  a  very 
different  behavior.  In  addition  to  the  temperature  rise  expected  from  the  increased  Plasmatron 
power,  each  specimen  underwent  a  further  spontaneous  surface  temperature  jump  of  several 
hundred  degrees.  The  times  at  which  these  temperature  jumps  initiated  are  indicated  by  open 
triangles  on  the  temperature  plots  of  Fig.  3.  Because  these  temperature  jumps  occurred  under 
constant  Plasmatron  operating  conditions,  they  must  be  associated  with  a  changing  energy 
balance  at  the  surface,  not  an  increase  in  the  freestream  enthalpy. 

Figure  4  provides  an  expanded  view  of  the  temperature  jump  region  for  the  two  different  SiC- 
containing  formulations;  two  specimens  have  been  removed  from  each  panel  for  clarity.  The 
cold- wall  heat  flux  for  samples  1.5  and  3.3  under  the  high-power  setting  was  -202  W  cm'  ,  just 
17  W  cm'2  higher  than  for  samples  1.8  and  3.9.  Figure  4  shows  that  both  samples  1.5  and  3.3 
approached  a  new  steady-state  surface  temperature  near  2215  K  to  accommodate  this  higher  heat 
input.  Then,  after  about  30  to  45  seconds  at  this  temperature,  the  surface  temperature  rose 
sharply  by  over  400  K  within  a  period  of  20  to  30  seconds. 
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Fig.  4.  Surface  temperatures  in  the  temperature  jump  region  measured  during  Plasmatron 
testing  of  ZrB2-30SiC  specimens  1.5, 1.6,  and  1.8  (top  panel)  and  (Zr,4W)B2-30SiC 
specimens  3.3, 3.5,  and  3.9  (bottom  panel). 


Similar  behavior  is  seen  for  samples  1.6  and  3.5,  but  with  the  temperature  jump  occurring 
sooner  after  the  switch  to  high-power  conditions.  For  specimens  1.6  and  3.5,  the  high-power 
settings  produced  cold-wall  heat  fluxes  of  236  and  253  W  cm'  ,  respectively,  compared  to  the 
-202  W  cm'  experienced  by  specimens  1.5  and  3.3.  These  higher  heating  conditions  cause  the 
temperatures  of  specimens  1.6  and  3.5  to  rise  more  rapidly  to  higher  levels,  and  the  time  to 
trigger  the  spontaneous  temperature  jump  is  shortened. 
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Fig.  5.  Still  images  of  the  temperature  jump  on  the  surface  of  (Zr,4W)  B2-30SiC  specimen 
3.3.  Images  are  3  seconds  apart;  the  measured  surface  temperatures  corresponding  to  the 
image  times  are  indicated  by  red  symbols  on  the  plot. 


The  temperature  jump  phenomenon  was  visually  evident  during  testing  as  a  large  increase  in 
light  emission  from  the  sample  surface.  Figure  5  shows  still  images  taken  from  a  video  of 
sample  3.3  showing  the  progression  of  the  temperature  jump  over  the  surface.  The  heat  flux  is 
largest  near  the  edges  of  the  test  specimen  because  of  the  abrupt  acceleration  of  the  flow  around 
the  specimen,  which  thins  the  boundary  layer  in  this  region.59  The  temperature  jump  initiates  on 
the  lower  left-hand  edge  of  the  specimen  and  progresses  around  the  edge  in  both  directions, 
completely  encircling  the  specimen  after  about  18  seconds  and  encompassing  the  entire  surface 
after  about  24  seconds.  This  same  progression  was  observed  for  the  other  four  specimens  that 
underwent  a  surface  temperature  jump,  although  the  times  to  complete  the  transition  varied. 

Note  that  the  camera  was  in  automatic  mode,  continuously  adjusting  to  accommodate  the 
increased  intensity  radiating  from  the  sample  surface.  Thus,  the  contrast  in  brightness  within 
each  image  is  a  better  indication  of  temperature  differences  than  comparisons  of  brightness 
between  different  images.  Since  temperature  is  measured  at  the  center  of  the  specimen  face  by 
the  two-color  pyrometer,  the  temperature  at  the  specimen  edge  as  the  transition  progresses  is 
likely  higher  than  the  value  plotted  in  Fig.  5. 
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B.  Surface  Oxidation 


Figure  6  shows  post-test  optical  micrographs  of  the  ZrB2-30SiC  specimens  1.9,  1.8,  1.5,  and 
1.6  and  (Zr,4W)B2-30SiC  specimens  3.4,  3.9,  3.3,  and  3.5.  Samples  1.9  and  3.4  were  exposed  to 
low-power  conditions  only  and  experienced  maximum  surface  temperatures  below  1900  K. 

Both  specimens  are  blackish  in  color  and  have  a  smooth,  shiny  surface.  Samples  1.8  and  3.9 
were  exposed  to  high-power  conditions,  reaching  surface  temperatures  near  2100  K,  but  did  not 
undergo  a  temperature  jump.  The  appearance  of  their  surfaces  is  very  similar  to  samples  1.9  and 
3.4,  still  relatively  smooth,  but  more  grayish  in  color.  Samples  1.5  and  3.3  experienced  the 
temperature  jump  phenomenon,  as  have  samples  1.6  and  3.5;  however  the  former  specimens 
spent  less  time  at  a  lower  post-jump  surface  temperature.  The  specimen  surfaces  are  whitish  and 
rougher,  especially  around  the  periphery.  The  peripheral  whitish  scale  of  sample  1.7  became 
detached  after  testing.  Poor  adherence  and  decreased  thermal  contact  of  this  peripheral  scale  to 
the  underlying  material  may  contribute  to  the  higher  temperatures  at  the  sample  edge  after  the 
temperature  jump  (as  indicated  by  the  images  in  Fig.  5). 


Fig.  6.  Optical  micrographs  of  ZrB2-30SiC  (top  panel)  and  (Zr,4W)B2-30SiC  (bottom 
panel)  specimens  after  Plasmatron  testing;  the  specimen  diameter  is  30  mm. 
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Fig.  7.  a)  XRD  spectrum  of  the  rim  of  sample  1.7  showing  the  crystal  structure  of 
monoclinic  ZrC>2;  b)  phosphorescence  of  sample  3.5  after  exposure  to  ultraviolet  light  from 
a  mercury  vapor  pen  light  indicating  the  presence  of  ZrC>2. 


Portions  of  scale  detached  from  sample  1.7  were  used  for  XRD  analysis  and  elemental 
analysis  for  Zr,  O,  and  Si  at  ATI  Wah  Chang.  XRD  analysis  confirmed  that  the  whitish  scale 
contains  monoclinic  zirconia  (Fig.  7a),  and  the  elemental  composition  reported  by  ATI  Wah 
Chang  is  consistent  with  ZrC>2  and  residual  Si  at  the  2900  parts-per-million  level  by  mass.  The 
presence  of  zirconia  was  also  indicated  by  the  observation  of  phosphorescence  after  ultraviolet 
light  illumination  (Fig.  7b).  Phosphorescence  is  a  known  property  of  both  Zr0260"61  and  HfC^.62 
Finally,  EDS  analysis  showed  the  presence  of  Zr,  O,  and  Si  in  various  relative  amounts  in 
different  regions  of  the  surface. 

These  analyses  indicate  that  after  the  temperature  jump,  the  rougher,  whitish  peripheral  oxide 
scale  is  mainly  composed  of  ZrC>2  with  only  minor  silica  remnants,  whereas  the  oxide  scale 
covering  the  sample  interior  is  smoother  containing  ZrC>2  islands  in  a  glassy  background.  The 
zirconia  rich  oxide  scale  on  the  sample  periphery  correlates  with  the  region  experiencing  the 
highest  aerothermal  heating59  and  the  location  where  the  temperature  jump  begins  and 
propagates  (see  Fig.  5). 

C.  Emittance 

One  possible  driving  force  for  the  temperature  jump  is  that  rapid  surface  oxidation  changes 
the  composition  and  morphology  of  the  surface  in  such  a  way  that  it  lowers  the  emittance 
substantially.  With  a  lower  emittance,  a  higher  temperature  is  required  to  maintain  the  steady- 
state  energy  balance  at  the  surface.  The  required  change  in  the  total  hemispherical  emittance  that 
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would  be  necessary  to  explain  the  observed  temperature  jumps  for  samples  1.5  and  3.3  can  be 
estimated  from  the  expression; 

^H=[TjTHYeL  ,  (1) 

where  the  subscripts  L  and  H  indicate  values  before  and  after  the  temperature  jump.  Equation 
(1)  is  derived  from  the  surface  energy  balance  by  assuming  that  the  surface  heat  flux  and  the 
conduction  losses  into  the  interior  are  constant:  qhw  =  £LoT*  +  qcond  =  £HoT ^  +  qcond  .  Plugging 
in  representative  temperatures  ( TL  =  2206  K  and  TH  =  2642  K  for  sample  1.5),  we  find 
sH  =  0.49^ ,  indicating  that  the  emittance  must  decrease  by  a  factor  of  ~2  if  it  is  to  account  for 
the  observed  temperature  jump. 


While  it  is  not  possible  to  measure  the  total  hemispherical  emittance  in  situ,  changes  in  the 
average  directional  emittance  (over  the  0.6-39  pm  wavelength  range,  47°  off  the  surface  normal) 
can  be  followed  by  combining  the  data  of  the  two-color  and  the  broadband  radiometers.  This 
average  directional  emittance  is  defined  by 


f e{T,?L,4T)Eh{T,X)dA 
s{T,  AT06_39,47°)  =  =  M - _ - 

06-39  jX(r,2>/2 


(2) 


where  L'06_i9  is  the  measured  radiance,  L06_39  is  the  ideal  blackbody  radiance,  and  Eb(T,X)  is 

the  Planck  blackbody  function.  The  broadband  radiometer  provides  a  direct  measurement  of 
L'q  6_39 ,  while  L0  6_39  is  computed  by  integrating  the  Planck  function  over  the  0.6-39  pm 

wavelength  range  at  the  surface  temperature  measured  by  the  two-color  radiometer. 
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Fig.  8.  Directional  emittance  in  the  temperature  jump  region  measured  during 
Plasmatron  testing  of  ZrB2-30SiC  specimens  1.5, 1.6,  and  1.8  (top  panel)  and 
(Zr,4W)B2-30SiC  specimens  3.3, 3.5,  and  3.9  (bottom  panel). 

The  0.6-39  pm  wavelength  range  contains  more  than  96%  of  the  total  emissive  power  for  any 
surface  temperature  in  the  1000-3000  K  range,  therefore  L'0  6_39  and  Z06_39  are  very  good 

approximations  for  total  spectral  values.  For  dielectric  surfaces  (e.g.,  oxides)  with  refractive 
indices  between  1  and  4,  electromagnetic  theory  predicts  differences  of  no  more  than  ~5% 
between  hemispherical  and  normal  emittance,  and  between  normal  emittance  and  any  directional 
emittance  less  than  50  degrees  off-normal.63  While  surface  roughness  may  cause  deviations 
from  these  predictions,  the  measured  average  directional  values  are  still  the  best  approximations 
currently  available  for  the  total  hemispherical  emittance  of  the  IJHTC  specimens  during  testing. 

Figure  8  shows  the  directional  emittance  derived  in  the  above  manner  for  the  same  specimens 
and  test  time  interval  as  in  Fig.  4.  The  emittance  prior  to  switching  to  the  high-power  condition 
is  very  high,  between  0.9  and  0.95  for  all  specimens.  Such  high  emittances  are  consistent  with 
estimates  made  by  Monteverde  and  Savino  during  arc-jet  testing  of  ZrB2-15SiC  materials  using  a 
similar  technique,41  as  well  as  reflectance  measurements  performed  on  ZrB2-30SiC  specimens 
after  Plasmatron  tests  in  which  sample  temperatures  remained  below  1900  K.46  After  the  switch 
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Table  2  Initial  specimen  mass  and 
mass  change. 


Sample 

Am,g 

1.4 

32.485 

-0.040 

1.5 

33.003 

-0.742 

1.6 

33.097 

-1.518 

1.7 

33.310 

-1.046 

1.8 

32.762 

-0.043 

1.9 

32.759 

+0.010 

2.4 

39.773 

+0.125 

2.6 

39.281 

-0.616 

2.10 

41.692 

-0.418 

2.11 

41.594 

+0.053 

3.2 

34.669 

-0.020 

3.3 

34.623 

-0.503 

3.4 

33.970 

+0.021 

3.5 

33.540 

-0.801 

3.9 

33.507 

-0.018 

to  high  power,  the  emittance  settles  to  a  slightly  lower 
value  between  0.85  and  0.9  for  specimens  1.8  and  3.9, 
which  did  not  undergo  a  temperature  jump. 

For  specimens  1.5  and  3.3,  the  emittance  stabilizes 
near  0.9  after  the  switch  to  high  power  and  before  the 
temperature  jump.  After  the  temperature  jump,  the 
emittance  decreases  and  shows  greater  fluctuation, 
however  it  never  falls  below  0.75,  and  for  the  most 
part  remains  higher  than  0.8.  This  result  argues 
strongly  against  decreasing  emittance  as  the  driving 
force  for  the  temperature  jump  phenomenon,  since  a 
much  lower  post-jump  emittance  would  be  required; 
i.e.,  sH  =  0.49  x  0.9  =  0.44 .  Similar  conclusions  can 

be  drawn  from  the  pre-  and  post-jump  emittances  of 
samples  1.6  and  3.5. 

D.  Mass  Changes 

The  changes  in  sample  mass  after  Plasmatron 
testing  are  listed  in  Table  2.  Most  specimens  lost 
mass,  but  several  gained  mass,  most  prominently, 
samples  2.4  and  2. 1 1 . 


The  oxidation  of  ZrB2  to  Zr02  and  B2C>3  results  in 
mass  increases  of  70%  if  both  oxides  remain  and  of  9.2%  if  B2C>3  completely  volatilizes.  The 
oxidation  of  (Zr,4W)B2  to  Zr02,  WO3,  and  B2C>3  results  in  mass  increases  of  69%  if  all  oxides 
remain,  9.5%  if  all  B2C>3  volatilizes,  and  1.5%  if  both  WO3  and  B2C>3  volatilize.  The  observed 
positive  mass  change  for  samples  2.4  and  2. 1 1  (which  reached  temperatures  of  -2130  K  and 
-1990  K,  respectively)  is  consistent  with  oxidation  experiments  performed  on  the  same 
(Zr,4W)B2  composition  in  air  at  1873  K. 19-20  The  large  mass  loss  observed  for  specimens  2.6 
and  2.10  cannot  be  accounted  for  by  oxide  volatilization  alone  and  must  be  due  to  the 
mechanical  loss  of  oxide  scale  (i.e.,  spallation  of  Zr02)  under  more  aggressive  heating 
conditions. 


The  oxidation  of  ZrB2-30SiC  to  Zr02,  Si02,  B2C>3,  and  CO(gas)  results  in  a  mass  change  of 
+67%  if  all  oxides  remain,  +16.5%  if  all  B2C>3  volatilizes,  and  -10.5%  if  both  Si02  and  B2C>3 
volatilize.  Thus  an  oxidation-driven  net  mass  loss  can  occur  with  sufficient  volatilization  of  the 
silica  rich  oxide  scale  at  high  temperatures.  Similar  behavior  is  calculated  for  the  (Zr,4W)B2- 
30SiC  composition. 

The  SiC-containing  specimens  tested  at  the  lowest  temperatures,  1.9  and  3.4,  show  small  mass 
gains,  while  those  reaching  higher  temperatures  without  undergoing  a  temperature  jump  - 
samples  1.4,  1.8,  3.2,  and  3.9  -  show  small  mass  losses.  This  result  is  consistent  with  previous 
Plasmatron  tests  of  the  ZrB2-30SiC  material  that  showed  no  net  mass  loss  below  1873  K,46  and 
indicates  that  the  production  of  a  Si02  oxide  scale  exceeds  the  rate  of  Si02  volatilization  below 
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this  temperature.  All  specimens  that  underwent  a  temperature  jump  -  specimens  1.5,  1.6,  1.7, 

3.3,  and  3.5  -  show  significantly  higher  mass  loss.  While  some  of  this  is  probably  due  to 
mechanical  removal  of  surface  oxides,  significant  volatilization  of  SiC>2  is  likely  when 
temperatures  exceed  2500  K  as  they  have  here. 

E.  Evidence  of  Volatilization 

Emission  spectroscopy  of  electronically  excited  B,  Si,  and  W  atoms  was  used  to  track  the 
appearance  and  evolution  of  volatilized  products  in  the  gas  phase.  These  atoms  serve  as 
convenient  signatures,  as  they  can  only  originate  from  the  test  article  and  have  sharp  spectral 
features  that  make  them  stand  out  from  background  molecular  band  emissions.  The  spectrometer 
channels  centered  at  249.92  nm,  288.16  nm,  and  400.  99  nm  were  used  to  follow  the 
B  (251/2->2Pi;23/2),  Si  {p°— >lD2)  and  W  (7P4° ->75’3 )  transitions.64  The  surface  temperatures  and 

raw  signal  intensities  at  these  three  emission  wavelengths  are  shown  together  in  Fig.  9  for  the 
tests  of  specimens  1.5  and  3.3.  Raw  signal  intensities  are  plotted  over  identical  ranges  in  all 
panels  and  a  logarithmic  scale  is  used  to  highlight  low  intensity  features. 
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Fig.  9.  Surface  temperature  and  raw  signal  intensity  for  characteristic  emission  lines  of 
boron,  silicon,  and  tungsten  atoms:  ZrB2-30SiC  specimen  1.5  (left-hand  panel)  and 
(Zr,4W)B2-30SiC  specimen  3.3  (right-hand  panel). 

During  the  initial  heating  phase  at  low  power,  the  rise  and  the  eventual  decay  of  the  boron 
emission  signal  is  observed  for  both  specimens  1.5  and  3.3.  This  type  of  emission  transient  was 
seen  before  in  Plasmatron  testing  of  ZrB2-30SiC  65  and  is  consistent  with  a  simple  model  of 
oxide  scale  development  that  begins  with  B2O3  formation  at  low  temperatures  followed  by  boron 
oxide  volatilization  and  silica  scale  growth  at  higher  temperatures.66  No  Si  atom  emissions  were 
observed  in  these  previous  tests.  Here,  a  trace  of  Si  emission  can  be  detected  for  sample  1 .5  and 
a  very  clear  Si  emission  signal  is  seen  for  sample  3.3.  Sample  3.3  reaches  temperatures  above 
1900  K  during  the  initial  heat  up  phase,  whereas  sample  1.5  remains  below  1850  K.  No  tungsten 
atom  emissions  are  detected  for  sample  1.5  (as  expected)  and  only  a  trace  can  be  seen  for 
sample  3.3. 

Upon  switching  to  the  high-power  condition,  the  emission  intensities  at  the  boron  and  silicon 
wavelengths  rise  sharply  by  a  factor  of  2  to  3.  The  increased  boron  and  silicon  emission  intensity 
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is  consistent  with  increased  B2O3  and  SiC>2  volatilization  rates  as  the  samples  approach  a  higher 
steady-state  surface  temperature  around  2215  K  in  response  to  the  new  ffeestream  conditions.  A 
small  step  change  in  intensity  at  the  tungsten  wavelength  is  seen  for  both  specimens.  Inspection 
of  individual  spectra  indicates  a  baseline  intensity  offset  between  the  low-power  and  high-power 
test  conditions  at  this  wavelength,  not  the  presence  of  tungsten  atom  emission. 

An  additional  factor  of  2  to  3  increase  in  the  B  and  Si  wavelength  emission  intensity  is 
observed  as  the  temperature  jump  occurs  on  each  sample.  The  signal  intensity  at  the  W 
wavelength  also  rises  for  both  samples.  In  the  case  of  sample  3.3,  the  presence  of  atomic 
tungsten  is  confirmed  by  an  abundance  of  strong  tungsten  atom  lines  in  individual  emission 
spectra.  These  characteristic  tungsten  features  are  missing  in  the  spectra  of  sample  1.5,  and  the 
measured  emission  at  400.99  nm  is  due  to  a  strongly  fluctuating  background.  These  emission 
data  confirm  substantially  increased  oxide  volatilization  coinciding  with  the  temperature  jump 
phenomenon,  which  is  coincident  with  the  substantial  mass  loss  measured  for  these  specimens. 

After  roughly  30  seconds  around  2660  K,  the  surface  temperature  of  sample  3.3  decays 
steadily,  and  the  B,  Si,  and  W  atom  emissions  follow  suit.  Similar  temperature  decays  were  seen 
for  samples  1.6  and  1.7  (see  Fig.  3),  and  similarly  coincident  decreases  in  Si  and  B  emission 
intensity  were  observed. 


Computational  Results 

The  computed  gas  composition  at  the  boundary  layer  edge  is  shown  in  Fig.  10  for  the  low- 
and  high-power  test  conditions  experienced  by  specimens  1.5,  1.6,  1.8,  3.3,  3.5,  and  3.9.  Oxygen 
is  essentially  fully  dissociated  at  all  Plasmatron  powers  used  in  this  study,  but  the  level  of 
nitrogen  dissociation  increases  between  the  low-power  and  high-power  settings.  The  mole 
fractions  of  O2  and  NO  are  always  less  than  0.1%  and  1%,  respectively,  and  are  effectively 
absent  from  Fig.  10.  The  O-atom  mole  fraction  remains  relatively  constant  around  30  -  35%, 
while  the  N2  and  N  mole  fractions  change  substantially  between  the  low-  and  high-power 
conditions.  The  gas  density  will  increase  as  the  hot  gas  boundary  layer  edge  gas  approaches  the 
cooler  specimen  surface,  and  both  O-atom  and  N-atom  mole  fractions  will  decrease  between  the 
boundary  layer  edge  and  the  sample  surface  because  of  gas-phase  and  surface  recombination 
reactions. 
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Fig.  10.  Boundary  layer  edge  gas  composition  for  low  (L)  and  high  (H)  power  test 
conditions  for  specimens  1.5, 1.6, 1.8,  3.3, 3.5,  and  3.9. 


The  heating  conditions  for  test  specimens  1.5,  1.6,  1.8,  3.3,  3.5,  and  3.9  are  listed  in  Table  3. 
The  second  and  third  columns  of  this  table  list  the  Plasmatron  power  and  the  measured  cold-wall 
heat  flux;  the  final  three  columns  list  average  values  of  the  surface  temperature,  emittance,  and 
hot- wall  heat  flux  for  30  second  intervals  characteristic  of  the  low  power,  high  power,  and  post¬ 
temperature  jump  conditions.  The  surface  temperature  and  hot-wall  heat  flux  values  are  used  to 
locate  each  test  on  the  heat  flux  abacus  computed  by  the  Boundary  Layer  Code,  as  described  in 
Section  IV.  Table  4  shows  the  derived  recombination  coefficients,  the  ratio  of  hot- wall  to  fully- 
catalytic  heat  flux,  and  the  computed  air  species  number  densities  at  the  specimen  surface  for 
each  test  specimen  under  low-  and  high-power  conditions. 


A 1-23 


Table  3  UHTC  heating  conditions. 


Sample 

Power, 

kW 

W  cm  2 

T 

±  w  ’ 

K 

£ 

4  k*} 

W  cm  2 

1.5 

200 

109 

1826 

0.949 

59.8 

353 

202 

2206 

0.901 

121 

2642 

0.842 

233 

1.6 

205 

109 

1865 

0.946 

65.0 

387 

236 

2841 

0.791 

292 

1.8 

212 

109 

1893 

0.924 

67.3 

292 

185 

2102 

0.898 

99.5 

3.3 

224 

111 

1889 

0.946 

68.4 

339 

202 

2223 

0.891 

123 

2669 

0.830 

239 

3.5 

203 

111 

1872 

0.929 

64.6 

402 

253 

2803 

0.870 

305 

3.9 

211 

109 

1877 

0.925 

65.1 

295 

185 

2119 

0.886 

101 

^From 

4/nv  = 

+  <lcond  >  with  qcond  =  0 

Table  4  Catalytic  efficiencies,  the  ratios  of  hot-wall  to  fully-catalytic  heat  flux. 

and  the  species  number  densities  at  the  surface. 

29  _ 3 

Number  Density,  xlO  m 

Sample 

r 

[N2] 

[N] 

[O] 

[02] 

[NO] 

1.5 

0.0031 

0.48 

23.90 

1.174 

14.16 

0.358 

0.088 

0.0099 

0.47 

16.28 

5.333 

10.83 

0.313 

0.093 

0.24 

0.94 

17.15 

0.046 

2.283 

5.918 

2.030 

1.6 

0.0051 

0.52 

23.90 

0.608 

13.52 

0.691 

0.177 

0.32 

0.95 

15.80 

0.077 

1.895 

5.551 

2.175 

1.8 

0.0063 

0.54 

23.67 

0.403 

13.02 

0.929 

0.245 

0.0062 

0.43 

16.88 

5.816 

11.49 

0.217 

0.060 

3.3 

0.0062 

0.54 

23.64 

0.504 

13.13 

0.855 

0.222 

0.0106 

0.48 

16.32 

5.084 

10.75 

0.337 

0.102 

0.37 

0.96 

17.02 

0.027 

1.432 

6.499 

2.167 

3.5 

0.0047 

0.51 

23.58 

0.826 

13.58 

0.567 

0.143 

0.16 

0.89 

15.84 

0.274 

3.890 

4.081 

1.758 

3.9 

0.0052 

0.53 

23.73 

0.574 

13.39 

0.716 

0.185 

0.0067 

0.44 

16.91 

5.578 

11.41 

0.235 

0.065 
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The  recombination  coefficients  for  the  low-power  conditions  range  from  3.1  -  6.3  x  10'  , 
slightly  higher  but  of  the  same  magnitude  as  previously  reported  for  the  ZrB2-30SiC  composition 
for  temperatures  around  1850  K.  Recombination  coefficients  for  samples  1.8  and  3.9  under 
high-power  conditions  are  near  the  top  end  of  this  same  range.  Somewhat  higher  recombination 
coefficients  of  ~1  x  10'  are  derived  for  samples  1.5  and  3.3  in  the  high-power  region  before  the 
temperature  jump.  Under  all  of  these  conditions,  the  observed  hot- wall  heat  flux  is  in  the  range 
of  43-54%  of  the  fully-catalytic  heat  flux. 

After  the  temperature  jump,  extremely  high  catalytic  efficiencies,  ranging  from  0.16  to  0.37, 
are  necessary  to  reproduce  the  observed  hot-wall  heat  fluxes.  The  observed  hot-wall  heat  flux  is 
upwards  of  89%  of  the  fully-catalytic  value.  Because  of  these  large  increases  in  catalytic 
efficiency,  the  computed  gas  composition  at  the  sample  surface  after  the  temperature  jump  is 
very  different  from  that  before  the  temperature  jump  (see  Table  4).  For  sample  1 .5,  the  N-atom 
number  density  drops  by  a  factor  of  -1 15  and  the  O-atom  concentration  by  a  factor  of  ~5. 
Similarly,  for  sample  3.3,  the  N-atom  number  density  drops  by  a  factor  of  -190  and  the  O-atom 
concentration  by  a  factor  of -7.5. 


Discussion 

In  this  section,  we  examine  the  sensitivity  of  our  conclusion  that  the  temperature-jump 
corresponds  to  a  large  increase  in  the  chemical  heating  component  of  heat  flux;  then  we  discuss 
exothermic  oxidation  reactions  as  sources  of  additional  chemical  heating  at  the  surface.  Finally, 
we  comment  on  UHTC  arc -jet  and  ICP  test  conditions  in  the  context  of  published  SiC  oxidation 
experiments  in  partially-dissociated  oxygen. 

A.  The  Temperature  Jump  and  Increased  Chemical  Heating 

A  heat  flux  abacus  for  the  high-power  conditions  experienced  by  samples  1.5  and  3.3  is 
shown  in  Fig.  11.  The  dashed  vertical  arrows  between  zero  and  the  y  =  0  contour  indicate  the 
convective  heating  component  and  the  solid  vertical  lines  arrows  between  the  y  =  0  contour  and 
the  observed  hot-wall  heat  flux  indicate  the  chemical  heating  component.  It  is  evident  that 
within  the  present  computational  model  the  experimental  post-temperature-jump  temperature  and 
heat  flux  can  only  be  explained  by  a  large  increase  in  the  chemical  component  of  heating. 
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Fig.  11.  Heat  flux  abacus  for  Pstat  =  10  kPa  and  He  =  18.6  MJ  kg'1.  The  horizontal  arrow 

indicates  the  temperature  jump  for  specimens  1.5  and  3.3;  the  dashed  and  solid  vertical 
arrows  show  the  convective  and  chemical  contributions  to  the  computed  heat  flux, 
respectively. 
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Several  factors  can  influence  the  location  of  the  experimental  (Tw,qhw)  point  on  the  heat  flux 

abacus.  An  error  in  surface  temperature  measured  by  the  two-color  radiometer  can  move  the 
experimental  point  to  the  right  or  left  on  the  abacus,  but  this  does  not  change  the  magnitude  of 
the  chemical  heating  component  necessary  to  reproduce  the  increased  heat  flux.  An  error  in  the 
radiance  measured  by  the  broad-band  radiometer  will  shift  the  experimental  heat  flux  (computed 
as  qhw  =  qmd )  up  or  down  on  the  abacus,  but  this  error  is  unlikely  to  exceed  5%,  and  does  not 

change  the  result  that  the  chemical  heating  contribution  must  increase  at  least  threefold. 
Neglecting  heat  conduction  loss  in  the  estimate  of  hot- wall  heat  flux  implies  that  the  point 
(K’VhJ  is  actually  lower  on  the  abacus  than  it  should  be.  If  the  conduction  loss  exceeded  6% 
of  the  radiative  heat  flux,  the  true  (Tw ,  qhw  )  point  would  be  above  the  fully  catalytic  ( y  =  1 ) 
boundary  for  both  samples  1.5  and  3.3  in  Fig.  11. 

The  abacus  itself  is  dependent  on  the  edge  conditions  rebuilt  using  the  calorimeter 
measurements.  The  largest  uncertainty  here  is  the  catalytic  efficiency  assigned  to  the  copper 
calorimeter  surface;  we  have  used  a  value  of  y  =  0.01 .  Krassilchikoff  et  al.68  developed  a 

combined  experimental/computational  procedure  to  derive  appropriate  catalytic  recombination 
efficiencies  for  a  cold  copper  surface  and  found  values  between  0.014  and  0.067.  In  general,  as 
the  calorimeter  catalytic  efficiency  is  increased,  the  rebuilt  boundary  layer  edge  enthalpy 
decreases,  and  a  higher  sample  catalytic  efficiency  is  required  to  reproduce  an  experimental 
(r  , ,  qhw  )  point.  A  sensitivity  study  performed  for  the  high  heat  flux  conditions  of  sample  1.5 
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showed  that  any  reference  catalytic  efficiency  greater  than  0.015  would  place  the  experimental 
post-temperature  jump  ( Tw,qhw  )  point  above  the  fully  catalytic  boundary.  Obviously,  this  would 

also  be  the  case  if  the  common  assumption  of  a  fully  catalytic  copper  surface  were  made. 

While  the  outlined  combination  of  experimental  uncertainties,  computational  limitations,  and 
analysis  assumptions  makes  the  extraction  of  accurate  numerical  values  for  surface  reaction 
energetics  problematic,  they  do  not  alter  the  conclusion  that  a  very  large  increase  in  chemical 
heating  is  required  to  explain  the  temperature  jump.  Moreover,  it  seems  probable  that  analysis 
corrections  would  tend  to  increase  (rather  than  decrease)  the  already  large  estimates  of  chemical 
heating  necessary  to  explain  the  temperature  jump. 

B.  Oxidation  and  Volatilization 

We  emphasize  that  the  boundary-layer  computations  used  to  compute  the  heat  flux  abacus  are 
limited  to  5 -species  air  and  that  surface  chemistry  is  constrained  to  a  simple  model  of 
independent  O-atom  and  N-atom  recombination  reactions.  Clearly  additional  chemical 
mechanisms  are  in  play  on  UHTC  surfaces  involving  the  formation  and  volatilization  of  silicon 
and  boron  oxides. 

Surface  oxidation  processes  that  transform  ZrB2  and  SiC  into  condensed  and  volatile  oxides 
are  all  exothermic  and  release  more  energy  per  oxygen  atom  consumed  than  catalytic 
recombination  to  molecular  oxygen  (A Hr  =  -256  kJ  mol'1).  On  the  other  hand,  processes  that 

lead  to  the  volatilization  of  condensed  oxides  are  all  endothermic.  As  an  example,  consider  the 
four  reactions: 


O  +  l/3SiC(s)  - 

^  l/3Si02(l)+  1/3CO 

A Hr  =  -565  kJ  mol1 

(Rl) 

O  +  l/2SiC(s)  -> 

►  l/2SiO  +  1/2CO 

A Hr  =  -340  kJ  mol'1 

(R2) 

Si02(l)  +  l/2SiC(s)  -+  3/2SiO  +  1/2CO 

A Hr  =  674  kJ  mol'1 

(R3) 

Si02(l)  -*•  Si02 

A Hr  =  545  kJ  mol'1 

(R4) 

These  reactions  describe,  respectively:  Rl)  the  passive  oxidation  of  SiC;  R2)  the  active  oxidation 
of  SiC;  R3)  the  decomposition  of  Si02  at  the  Si02/SiC  interface;  and  R4)  the  congruent 
evaporation  of  SiC>2.  The  standard  reaction  enthalpies  listed  were  computed  with  the  HSC 
Chemistry  program69  for  2500  K,  but  values  do  not  change  much  with  temperature  in  the  2000  to 
3000  K  range.  Similar  reactions  can  be  written  for  ZrB2  oxidation  and  the  decomposition  and 
volatilization  of  boron  oxides. 

More  energy  is  released  (per  O-atom  consumed)  during  passive  than  active  SiC  oxidation. 
However,  if  the  Si02/SiC  interface  remains  stable  and  the  Si02  evaporation  rate  is  low, 
condensed  Si02  forms  an  effective  oxygen  diffusion  barrier,  slowing  further  oxidation,  and  in 
turn  reducing  the  rate  at  which  energy  is  released  at  the  reaction  interface.  Under  such  passive 
oxidation  conditions,  exothermic  oxidation  reactions  and  endothermic 

evaporation/decomposition  reactions  have  little  impact  on  the  heat  balance  at  the  surface,  and  the 
chemical  component  of  heat  flux  is  dominated  by  surface  catalytic  reactions.  This  situation 
characterizes  the  low-power  test  conditions  of  this  study. 
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If  the  sample  temperature  is  raised,  the  rates  of  both  the  oxidation  and  volatilization  reactions 
increase,  and  it  is  possible  to  achieve  a  quasi-steady-state  where  condensed-phase  oxide  is 
formed  and  volatilized  at  similar  rates,  and  where  the  net  contribution  of  these  reactions  to  the 
chemical  heat  flux  remains  largely  neutral.  This  regime  may  describe  the  behavior  of  samples 
1.4,  1.8,  3.2,  and  3.9  under  the  high-power  test  conditions.  These  samples  all  lost  some  mass, 
but  the  derived  recombination  coefficients  (for  samples  1.8  and  3.9)  do  not  differ  much  from 
samples  tested  under  the  low-power  condition. 

In  contrast,  under  the  high-power  test  conditions  of  samples  1.5  and  3.3,  emission 
spectroscopy  demonstrates  that  the  silicon  and  boron  species  volatility  is  high  prior  to  the 
temperature  jump  and  increases  even  further  after  the  temperature  jump.  A  plausible  scenario  is 
that  the  existing  borosilicate  oxide  scale  is  steadily  eroded  during  the  pre-temperature  jump 
induction  period  by  decomposition  reactions  at  the  scale/substrate  interface  and/or  evaporative 
processes  at  the  scale  surface.  Once  the  oxide  scale  is  sufficiently  eroded  -  whatever  the 
mechanism(s)  -  the  scale  no  longer  limits  oxygen  diffusion  effectively  causing  exothermic 
oxidation  reactions  to  accelerate.  Consequently,  their  contributions  to  the  chemical  heat  flux  will 
increase  and  the  surface  temperature  will  rise.  The  oxidation  reaction  rates  will  increase  further 
with  rising  temperature,  and  this  process  can  “run  away”  until  a  new  steady-state  is  reached 
where  oxidation  produces  primarily  Zr02  and  gaseous  silicon  and  boron  oxides,  and  where  the 
higher  chemical  heating  rate  is  again  balanced  by  radiation  losses. 

In  this  scenario,  the  observation  of  a  temperature  jump  is  directly  related  to  existence  of  a 
well  developed  silica-rich  oxide  scale  that  retards  the  response  time  of  the  UHTC  specimen  to 
the  switch  from  low-power  to  high-power  heating  conditions,  the  latter  of  which  favors  active 
over  passive  oxidation.  The  time  required  for  the  existing  oxide  scale  to  decompose/volatilize 
would  be  shorter  for  larger  power  increases,  consistent  with  the  results  for  samples  1.6,  1.7,  and 
3.5  versus  samples  1.5  and  3.3.  If  a  sample  were  exposed  directly  to  a  high-power  condition 
favoring  active  oxidation,  no  temperature  jump  would  be  expected,  since  the  sample  would 
transition  smoothly  into  the  active  oxidation  regime  during  heating.  This  latter  scenario  may 
also  explain  why  no  temperature  jump  was  observed  for  the  (Zr,4W)B2  samples  and  why  they 
reached  higher  surface  temperatures  than  the  ZrB2-30SiC  and  (Zr,4W)B2-30SiC  specimens  under 
similar  low-power  conditions;  under  the  low-power  conditions  of  this  study  boron  oxide 
volatilization  is  already  substantial  enough  that  no  oxidation  protection  is  afforded  by  a 
condensed  B2O3  scale,  allowing  active  exothermic  oxidation  reactions  (for  example, 

O  +l/5ZrB2(s)=  l/5ZrC>2(s)+  B2O3;  A Hr  =  -574  kJ  mol"1)  to  operate  effectively. 

It  is  certainly  possible  that,  in  addition  to  oxidation  reactions,  other  exothermic  surface  and 
gas-phase  reactions  help  sustain  the  high  surface  temperatures  observed  after  the  temperature 
jump.  One  possible  heat  source  is  enhanced  O-atom  and  N-atom  surface  recombination  on 
zirconia  versus  silica.  Balat-Pichelin  et  al.58  have  measured  the  O-atom  recombination 
coefficient  on  zirconia  surfaces  and  found  values  on  the  order  of  0.1 -0.2  between  2000  K  and 
2400  K.  (The  recombination  coefficient  for  N  atoms  on  zirconia  has  apparently  never  been 
measured.)  Increased  N-atom  recombination  on  SiC  versus  SiC>2  was  used  by  Herdrich  et  al.36  to 
explain  part  of  the  increased  heat  flux  to  the  surface  after  temperature  jumps  associated  with  the 
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PA  transition  on  SiC.  A  further  possibility  is  that  volatilized  Si  or  B  species  react  exothermically 
with  O  or  N  in  the  boundary-layer.  For  example,  Hald70  has  suggested  a  catalytic  cycle 
involving  the  exothermic  gas-phase  reactions  Si  +N  — ►  SiN  followed  by  SiN  +  N  — ►  N2  +  Si  to 
explain  high  surface  temperatures  after  temperature  jumps  observed  for  C/SiC  materials.  Such  a 
mechanism  would  raise  the  temperature  of  the  boundary  layer  and  hence  the  convective 
component  of  heating  to  the  surface.  Unraveling  these  different  possibilities  will  require  the 
development  of  appropriate  materials  response  and  CFD  models  that  include  all  relevant  gas- 
phase  and  surface  chemistries. 

C.  Test  Conditions  and  the  PA  Transition 

The  transition  boundary  between  the  passive  and  the  active  oxidation  regimes  of  SiC  is 
typically  discussed  as  a  function  of  temperature  and  oxygen  pressure.  Even  in  molecular 
oxygen,  the  precise  location  of  this  boundary  varies  significantly  between  different  researchers. 
Laboratory  tests  by  Balat  '  and  Rosner  and  Allendorf  in  low-pressure  partially  dissociated 
flows  suggest  that  the  boundary  is  pushed  to  lower  oxygen  pressures  at  a  given  temperature;  i.e., 
that  the  temperature-pressure  parameter  space  for  passive  oxidation  is  larger  in  the  presence  of 
atomic  oxygen. 

Rosner  and  Allendorf  conducted  experiments  in  N2-Ar  at  a  low  total  pressure  (-133  Pa)  and 
O-atom  partial  pressures  on  the  order  of  0. 1  to  1  Pa.  Balat32  conducted  experiments  in  partially- 
dissociated  air  at  total  pressures  of  600  to  5000  Pa  and  O-atom  partial  pressures  in  the  range  20 
to  2000  Pa.  Our  test  conditions  are  closer  to  those  of  Balat  than  those  of  Rosner  and  Allendorf. 
The  computed  O-atom  partial  pressures  at  the  sample  surface  are  about  3300  to  3600  Pa  for  all 
our  pre-temperature  jump  test  conditions.  After  the  switch  to  high  power,  samples  1 .4,  1 .8,  3.2 
and  3.9  reached  temperatures  near  21 10  K  and  remained  stable,  whereas  samples  1.5  and  3.3 
reached  temperatures  near  2215  K  and  eventually  underwent  the  temperature  jump.  A  transition 
temperature  near  2150  K  for  our  computed  O-atom  pressures  is  consistent  with  the  transition 
boundary  published  by  Balat32  for  tests  on  CVD  SiC,  but  is  significantly  higher  (by  -300  K)  than 
the  transition  boundary  for  sintered  SiC  published  by  the  same  author. 

Previous  tests  of  the  ZrB2-30SiC  composite  were  conducted  in  the  VKI  Plasmatron  facility 
using  the  same  mass  flow  and  static  pressure,  but  at  lower  powers  and  on  flat  faced  samples  held 
in  the  standard  ESA  50  mm  stagnation  point  test  fixture.46'47  Maximum  sample  temperatures 
remained  below  1900  K  in  these  tests  and  no  temperature  jumps  were  observed.  This  result  is 
also  consistent  with  Balat’ s  transition  boundary  for  CVD  SiC. 

It  is  difficult  to  compare  other  previous  UHTC  tests  with  predicted  SiC  PA  transition 
boundaries  because  the  oxygen  partial  pressures  at  the  specimen  surface  are  not  determined. 

The  example  cited  by  Glass,  of  a  sudden  temperature  jump  during  UHTC  arc -jet  testing  for 
the  SHARP  B1  program,40  occurred  over  about  15  seconds  with  the  temperature  rising  from 
from  2360  to  2810  K.  The  arc -jet  operating  conditions  and  the  ZrB2-SiC  composition  were  not 
specified,  but  the  temperature  profile,  the  temperature  at  which  the  jump  initiated,  and  the 
magnitude  of  the  temperature  jump  are  all  very  similar  to  those  observed  for  samples  1.6  and  3.5 
in  this  study. 
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Monteverde,  Savino,  and  colleagues  have  tested  ZrB2-15SiC  and  ZrB2-20SiC  composites  in 
the  form  of  hemispheres  and  cones  in  the  arc -jet  facility  at  the  University  of  Naples  using  step¬ 
wise  increases  in  power  resulting  in  stepwise  increases  in  measured  surface  temperature.41'44 
Surface  temperatures  approached  -2200  K  under  subsonic  atmospheric  pressure  conditions41  and 
-2100  K  in  supersonic  low  pressure  flows,42'44  but  no  evidence  of  a  temperature  jump  is  found  in 
published  experimental  temperature  profdes. 

Zhang  et  al.45  have  tested  flat-faced  ZrB2-20SiC  specimens  in  an  arc -jet  environment  at  cold- 
wall  heat  fluxes  of  170  W  cm'2  and  540  cm'2.  Specimens  reached  peak  temperatures  of -1930  K 
and  2603  K  under  these  conditions,  respectively,  and  no  temperature  jump  is  seen  in  their 
heating  profdes.  However,  the  two  maximum  temperatures  lie  below  and  above  the  -2200  K 
level  where  we  observe  the  phenomenon  in  our  experiments,  and  it  is  possible  that  the 
temperature  jump  is  simply  part  of  the  initial  temperature  rise  under  the  high  heat  flux  condition. 
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Gasch  et  al.  have  heated  HfB2-20SiC  specimens  in  an  arc-jet  to  temperatures  near  2100  K, 
and  documented  a  slow  but  steady  rise  in  surface  temperature  starting  after  about  6  minutes. 

This  temperature  rise  was  on  the  order  of  300  to  400  K  over  several  minutes,  and  was  attributed 
to  the  slow  volatilization  of  silica  from  the  scale  leaving  a  residual  Hf02  surface  of  lower 
emittance,  higher  catalycity,  and  lower  thermal  conductivity. 

The  temperature  jump  phenomenon  is  likely  a  function  of  multiple  material  and 
environmental  parameters:  the  composition  of  the  UHTC  material,  the  temperature,  pressure,  and 
composition  of  the  gas-phase,  sample  configuration  and  ffeestream  flow  characteristics  that 
determine  reactant  and  product  transport,  etc.  More  systematic  studies  that  isolate  and  vary 
specific  parameters  are  required  for  a  better  understanding  of  the  conditions  and  mechanisms  that 
trigger  the  temperature  jump.  Since  heat  and  mass  transfer  at  the  gas-surface  interface  result 
from  the  coupled  transport  and  thermochemistry  occurring  in  the  boundary  layer  at  the  reacting 
surface,  it  is  difficult  to  evaluate  various  potential  contributions  until  the  candidate  chemical 
mechanisms  and  rate  coefficients  are  incorporated  into  the  computational  tools  used  to  simulate 
the  experiment. 


Summary  and  Conclusion 

The  experimental  observations  can  be  summarized  as  follows.  The  ZrB2-30SiC  and 
(Zr,4W)B2-30SiC  UHTC  specimens  heat  up  and  oxidize  as  they  are  exposed  to  the  high-enthalpy 
air  flow  in  the  Plasmatron.  Below  a  certain  Plasmatron  power  level,  the  specimens  reach  a 
steady-state  temperature  at  which  they  form  a  stable  surface  oxide  that  slows  progressive 
oxidation.  Some  volatilization  of  boron  and  silicon  can  be  detected  by  emission  spectroscopy, 
but  net  sample  mass  changes  are  small  and  the  oxides  scales  retain  a  relatively  smooth, 
homogeneous  surface.  Above  a  certain  Plasmatron  power,  specimens  reach  a  temperature  at 
which  their  oxide  scale  is  no  longer  protective  or  stable.  For  the  current  tests,  conducted  at  a 
static  pressure  of  10  kPa,  oxides  at  steady-state  surface  temperatures  around  2100  K  remained 
stable,  while  those  heated  to  2200  K  became  unstable  after  a  short  induction  period.  The 
(Zr,4W)B2  specimens  reached  higher  surface  temperatures  than  the  ZrB2-30SiC  and  (Zr,4W)B2- 
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30SiC  materials  for  similar  Plasmatron  powers  and  did  not  exhibit  a  temperature  jump 
phenomenon. 

Emission  spectroscopy  during  testing  of  ZrB2-30SiC  and  (Zr,4W)B2-30SiC  materials  shows  a 
significant  increase  in  boron  and  silicon  atom  emissions  when  the  Plasmatron  power  is  switched 
from  the  low-  to  high-power  conditions,  reflecting  increased  oxide  volatility  above  2100  K.  A 
further  large  increase  in  emission  intensity  is  observed  as  the  temperature  jump  occurs.  Video 
images  show  that  the  temperature  jump  initiates  at  a  point  on  the  specimen  edge,  and  then 
progresses  around  the  circumference  and  fills  to  the  center  of  the  specimen  face. 

By  combining  the  measurements  of  the  two-color  and  broadband  radiometer,  the  total 
emittance  is  shown  to  decrease  slightly  when  the  temperature  jump  occurs,  however  this  change 
in  emittance  is  not  sufficient  to  be  the  driving  force  for  the  temperature  jump.  Aerothermal 
heating  computations  based  on  the  rebuilt  freestream  conditions  for  the  high-power  Plasmatron 
tests  indicate  that  the  observed  temperature  jumps  require  a  large  increase  in  the  chemical 
contribution  to  heat  flux  to  reproduce  the  experimental  combination  of  surface  temperature  and 
heat  flux. 

Taken  in  aggregate,  the  experimental  observations  and  computational  analysis  support  the 
hypothesis  that  the  temperature  jump  is  associated  with  a  transition  in  surface  chemistry 
involving  the  Si-containing  compounds  -  the  silica-rich  oxide  scale  and  the  silica  former  SiC  - 
and  is  likely  related  to  the  passive-to-active  oxidation  transition  observed  for  other  SiC- 
containing  composites  in  aerothermal  test  environments.  A  deeper  understanding  of  the 
temperature  jump  phenomenon  in  UHTC  materials  will  require  additional  testing  to  span  a  wider 
temperature,  pressure,  and  composition  parameter  space,  and  the  development  of  appropriate 
materials  response  and  CFD  models  to  capture  observed  oxidation  transition  behavior  in 
dissociated  gas  flows. 
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Abstract 

The  thermal  and  electrical  transport  properties  of  various  spark-plasma-sintered  HfB2- 
and  ZrB2-based  poly  crystalline  ceramics  were  investigated  experimentally  over  the  300-700  K 
temperature  range.  Measurements  of  thermal  diffusivity,  electrical  resistivity,  and  Hall 
coefficient  are  reported,  as  well  as  the  derived  properties  of  thermal  conductivity,  charge  carrier 
density,  and  charge  carrier  mobility.  Hall  coefficients  were  negative  confirming  electrons  as  the 
dominant  charge  carrier,  with  carrier  densities  and  mobilities  in  the  3-5x10  cm'  and  100-250 
cm2  V"1  s'1  ranges,  respectively.  Electrical  resistivities  were  lower,  and  temperature  coefficients 
of  resistivity  higher,  than  those  typically  reported  for  HfB2  and  ZrB2  materials  manufactured  by 
conventional  hot-pressing.  A  Wiedemann-Franz  analysis  confirms  the  dominance  of  electronic 
contributions  to  heat  transport.  The  thermal  conductivity  was  found  to  decrease  with  increasing 
temperature  for  all  materials.  Results  are  discussed  in  terms  of  sample  morphology  and 
compared  to  data  previously  reported  in  the  scientific  literature. 

I.  Introduction 

Ultra-high  temperature  ceramic  (UHTC)  composites  based  on  HfB2  and  ZrB2,  together 
with  minor  silica  forming  constituents  such  as  SiC,  SfAL,  TaSi2,  and  MoSi2,  are  under 
investigation  for  use  in  aerothermal  heating  environments  as  sharp  leading  edge  components  on 
future  generations  of  hypersonic  reentry  vehicles.1'3  Transition  metal  diborides  have  mechanical 
properties  and  brittle  fracture  behavior  typical  of  ceramics,  and  yet  their  electrical  and  thermal 
conductivities  are  more  characteristic  of  a  metal.4  The  ability  to  effectively  conduct  heat  is  a 
desirable  property  for  sharp-leading-edge  components,  as  it  improves  their  thermal  shock 
resistance  by  reducing  temperature  gradients  and  thermal  stresses  within  the  material,  and 
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transports  energy  away  from  the  stagnation  point  over  a  larger  component  surface  area  from 
which  it  can  be  efficiently  radiated  back  to  the  environment.5 

As  typical  of  polycrystalline  materials,  the  effective  thermal  and  electrical  conductivities 
of  UHTC  materials  are  affected  by  factors  such  as  chemical  composition,  grain  size,  and 
porosity.  The  densification  of  pure  ZrB2  and  HfB2  materials  by  the  conventional  hot  pressing  of 
compound  powders  is  difficult  because  of  their  hardness  and  very  high  melting  points.  Additives 
such  as  SiC,  B4C,  and  carbon  can  aid  the  sintering  process,  however  extreme  conditions  of 
pressure  and  temperature  are  generally  required.  Moreover,  processing  times  can  still  reach 
hours  which  encourages  grain  growth,  undesirable  from  a  mechanical  property  perspective. 
More  recently,  spark-plasma-sintering  of  elemental  and  compound  powders,  in  which  a  pulsed 
DC  current  is  applied  during  pressing,  has  been  shown  to  produce  dense  ZrB2  and  HfB2 
composites  at  less  extreme  temperature  and  pressure  conditions,  and  in  shorter  (tens  of  minutes) 
processing  times.6'8 

Here  we  investigate  the  thermal  and  electrical  transport  properties  of  nine  spark-plasma- 
sintered  (SPS)  HfB2-  and  ZrB2-based  poly  crystalline  ceramics  over  the  300-700  K  temperature 
range.  Thermal  diffusivity,  electrical  resistivity,  Hall  Effect  measurements,  and  microscopy  are 
performed  on  the  same  specimen  of  each  ceramic.  Thermal  conductivity,  charge  carrier  density, 
and  charge  carrier  mobility  are  derived  from  the  measured  data.  The  measurement  of  both 
thermal  and  electrical  properties  allows  the  analysis  of  electron  and  phonon  contributions  to 
thermal  transport. 

II.  Experimental  Procedure 
(1)  Materials  and  Processing 

Specimens  from  nine  different  UHTC  billets  were  tested;  seven  billets  were  HfB2-based 
materials  and  two  were  ZrB2-based  materials.  We  have  grouped  the  nine  samples  into  three 
characteristic  groups.  Group  (A)  consists  of  nominally  pure  diboride  materials;  group  (B) 
consists  of  the  composites  prepared  with  excess  elemental  constituents:  B1  (HfB2-l%  Ir),  B2 
(HfB2-4%  Hf  or  HfBi.9),  and  B3  (HfB2-2.5%  B  or  HfB2.i);  group  (C)  consists  of  composites  with 
SiC  additions:  Cl  (HfB2-5%  SiC),  C2  (HfB2-5%  SiC),  and  C3  (ZrB2-20%  SiC).  Volume 
percentages  are  used  throughout  this  paper.  Sample  labels,  target  compositions,  and  processing 
identifiers  are  summarized  in  the  second  and  third  columns  of  Table  1. 

The  diboride  phase  was  produced  from  the  reaction  of  elemental  metal  (Hf  or  Zr)  and 
boron  in  all  but  two  billets  (Cl  and  C3)  which  were  produced  from  milled  diboride  powders. 
The  following  raw  powders  were  used:  HfB2  (-325  mesh,  99.5%,  CERAC,  Milwaukee,  WI), 
ZrB2  (-325  mesh,  99.5%,  CERAC),  Hf  (-325  mesh,  99.8%,  CERAC),  Zr  (-325  mesh,  99.7%, 
CERAC),  B  (amorphous,  -325  mesh,  99%,  Alpha  Aesar,  Ward  Hill,  MA),  SiC  (UF-05,  99.8%, 
H.  C.  Starck,  Newton,  MA),  and  Ir  (-325  mesh,  99.95%,  Surepure  Chemetals,  Florham  Park, 
NJ). 

All  elemental  powders  were  used  as-purchased  without  additional  milling.  Powders  were 
combined  in  desired  proportions  and  hand  mixed  prior  to  pressing,  either  on  the  bench  top  or 
within  a  glove  box  to  minimize  humidity  (processing  method  “dry#”  in  Table  1).  SiC  powder 
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was  also  used  as-purchased  when  combined  with  elemental  powders,  as  for  specimen  C2.  For 
the  billets  made  using  diboride  powders  (processing  method  “/c-1”  in  Table  1),  the  raw  diboride 
and  SiC  powders  were  weighed  and  combined  in  the  desired  volumetric  ratio,  and  then  wet- 
milled  in  cyclohexane  with  WC  milling  media  in  a  planetary  mill  (Fritsch  Pulverisette  5, 
Germany).  The  milled  powders  were  carefully  dried  to  prevent  phase  segregation  between  the 
SiC  and  the  denser  diboride  components. 

Table  I.  Sample  groups,  target  volumetric  compositions,  processing  methods,  densities 
(experimental,  theoretical,  and  their  ratio),  mean  intercept  length  (L  ±cr)  from  log¬ 
normal  fitting,  and  grain  boundary  area  per  unit  volume  (  Sv ). 


Group 

ID 

Target 

Composition 

Processing 

ID 

P  exp 

(g  cm'3) 

P the 

(g  cm'3) 

P exp  ^  P the 

% 

L  ±<7 
(pm) 

Sv 

(pm'1) 

A1 

HfB2 

TCI/  — 

11.0 

11.21 

98.1 

10.7  ±4.9 

0.205 

A2 

HfB2 

TC3/  dry4 

10.3 

11.21 

91.9 

5.5  ±3.2 

0.386 

A3 

ZrB2 

TC6/  dry2 

5.66 

6.12 

92.5 

5.6  ±3.2 

0.368 

B1 

HfB2-l%Ir 

TC5/  — 

11.1 

11.33 

98.0 

10.3  ±4.2 

0.200 

B2 

HfB2-4%Hf  f 

TC7/  dryl 

11.0 

11.29 

98.1 

10.0  ±5.1 

0.218 

B3 

HfB2-2.5%B  f 

TC8/  dryl 

11.1 

10.99 

101.0 

10.9  ±5.8 

0.191 

Cl 

HfB2-5%SiC 

TC2/  c-1 

10.7 

10.81 

99.0 

5.4  ±2.2 

0.387 

C2 

HfB2-5%SiC 

TC4/  dryl 

11.0 

10.81 

101.8 

5.5  ±3.3 

0.415 

C3 

ZrB2-20%SiC 

TC9/  c-1 

5.50 

5.54 

99.3 

A:  10.2  ±5.5 
B:  3.9  ±  1.9 

AA:  0.215 
AB:  0.206 
Total:  0.421 

f  HfB2-4%Hf  =  HfBi.9  and  HfB2-2.5%B  =  HfB2.i 
A:  ZrB2,  B:  SiC 


All  billets  were  consolidated  by  spark-plasma  sintering  (Model  SPS-1050,  Sumitomo 
Heavy  Industries  Ltd.,  Japan)  at  the  University  of  California  Davis.  Powders  were  packed  into 
20  mm  diameter  graphite  dies,  loaded  to  105-135  MPa,  and  densified  at  temperatures  from  1700 
to  1900°C,  with  hold  times  of  5-10  minutes.  Heating  rates  during  the  SPS  process  were  100 
to300°C  per  minute.  After  pressing,  test  specimens  were  diamond-machined  from  the  billets  into 
1  mm  thick,  12.7  mm  diameter  disks.  The  thermal  and  electrical  property  measurements  were 
performed  on  the  same  specimen  from  each  billet. 

(2)  Sample  Characterization 

The  densities  of  the  test  specimens  were  measured  using  the  Archimedes  method.  The 
specimen  crystal  structure  was  characterized  using  an  X-ray  diffraction  (XRD)  apparatus 
(Phillips  Electronics,  New  York,  NY)  with  a  Cu  Ka  source  and  a  Ge  monochromator.  After  the 
completion  of  all  thermal  and  electrical  property  testing,  specimens  were  polished  to  a  1  -pm 
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finish  (Struers,  RotoPol-31,  Germany),  chemically  etched  with  molten  potassium  hydroxide 
(99.99%,  Alpha  Aesar,  Ward  Hill,  MA)  and  fractured  into  smaller  specimens.  The 
microstructure  of  both  the  etched  and  fractured  specimen  surfaces  were  characterized  using  a 
JEOL  6100  high-resolution  field-emission  scanning  electron  microscope  (JEOL  Ltd.,  Tokyo, 
Japan)  with  accompanying  setup  (EDAX,  Mahwah,  NJ)  for  elemental  analysis  by  energy 
dispersive  spectroscopy  (EDS). 


Grain  size  analysis  was  performed  on  scanning  electron  microscope  (SEM)  images  using 
the  ImageJ  program  with  the  Concentric  Circles  plug-in.9  Two  circles  of  known  radius  were 
superimposed  on  each  SEM  image  and  the  arc  length  crossing  each  grain  was  determined  by 
measuring  the  central  angle  of  the  arc  and  converting  it  into  an  intercept  length,  L  .  This 
procedure  is  similar  to  those  found  in  ASTM  standard  E  112-96. 10  Typically,  200-300  intercept 
length  measurements  were  obtained  for  each  etched  specimen.  From  these  measurements  an 
intercept  length  probability  density  histogram  with  1pm  binning  was  constructed  for  each 
sample,  and  the  histogram  was  fit  by  a  log-normal  distribution  function: 
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where  the  fitting  parameters  //  and  cr  are  the  mean  and  standard  deviation  of  In  L .  The  mean 
and  standard  deviation  of  L  are  calculated  from  //  and  a  by: 
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The  mean  grain  size  D  is  of  similar  magnitude  to  the  mean  intercept  length  L  ,  and  various 
geometric  approximations  are  available  to  make  this  conversion  if  desired.11 

The  grain  boundary  interface  area  per  volume,  Sv  ,  was  also  calculated  for  each 
specimen  from  the  stereological  relationship  Sv  =  2 NGB/Ltot  ,  where  NGB  is  the  total  number  of 
grain  boundaries  intersected  by  the  two  circles  and  Ltot  is  the  sum  of  the  two  circle 
circumferences.12 


(3)  Thermal  Diffusivity  Measurements 

Thermal  diffusivity  was  measured  using  a  photothermal  radiometry  technique  in  which 
the  front  face  of  a  thin  disk  was  periodically  heated  using  a  laser  and  the  harmonic  back-face 
temperature  response  was  recorded  and  modeled  to  extract  the  thermal  diffusivity.  Photothermal 
techniques  have  been  extensively  developed  and  are  widely  used  to  measure  the  thermal 

1  O  1C 

properties  of  bulk  and  composite  materials.  '  The  experimental  arrangement  used  in  our 
measurements  is  shown  in  Fig.  1. 

The  UHTC  specimen  was  mounted  by  its  edge  in  a  cylindrical  aluminum  holder  fit  with  a 
band  heater  capable  of  heating  specimens  to  about  450°C.  Specimen  temperature  was  measured 
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by  a  type  K  thermocouple  in  direct  contact  with  the  specimen  surface.  A  Synrad  Firestar  F400 
CO2  laser  (output  wavelength  10.6  pm),  coupled  to  a  Stanford  Research  Systems  DG535  digital 
delay  and  pulse  generator  for  output  modulation,  was  used  to  irradiate  the  front  face  of  the 
sample.  A  3x  beam  expander  (Haas  Laser  Technologies,  Flanders,  NJ)  was  used  to  expand  the 
infrared  beam  diameter  from  about  3.5  mm  at  the  laser  exit  to  about  11  mm  on  the  sample 
surface.  This  expansion  made  the  Gaussian-like  laser  beam  diameter  much  larger  than  the  1  mm 
sample  thickness,  assuring  that  heat  transfer  between  the  front  and  back  faces  of  the  sample  was 
quasi  one-dimensional  near  the  center  of  the  disk.  Thermal  emission  originating  from  the  center 
of  the  back  face  was  focused  by  a  CaF2  lens  onto  a  liquid-nitrogen  cooled  InSb  detector  (Infrared 
Laboratories  Inc,  Model  ND-2).  The  detector  output  voltage  was  processed  by  a  Stanford 
Research  Systems  SR830  lock-in  amplifier  to  extract  its  amplitude  and  phase  shift  as  a  function 
of  laser  modulation  frequency.  Since  the  amplitude  of  the  harmonic  temperature  component, 
STb  ,  was  small  (~0.1  K),  the  detected  emissive  power  variations  were  essentially  linearly 

proportional  to  STb ,  as  can  be  verified  from  a  blackbody  radiation  analysis. 


Fig.  1.  Schematic  of  the  photothermal  radiometry  setup  for  thermal  diffusivity  measurement. 


A  one-dimensional  heat  conduction  model  was  used  to  determine  the  thermal  diffusivity 
of  the  test  specimens  from  the  measured  thermal  emission.  The  harmonic  component  of  the 
back-face  temperature  response  to  periodic  front  face  heating  is  given  by  the  complex 
expression15 


A  exp  (icot) 

a  cr[exp(cr<f )  -  exp(-crc/)]  ’ 


(4) 


where  a  is  the  thermal  diffusivity,  d  is  the  sample  thickness,  t  is  time,  and  co  =  lKf  with  / 
the  laser  heating  frequency  in  Hz.  In  Eq.  (4),  a  =  \l  +  i]y[a>/2a  and  the  leading  constant 
A  =  /0(1  -  R)j pcp  ,  where  I0  is  the  laser  intensity,  R  is  the  reflectance  of  the  illuminated 
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surface,  p  is  the  sample  density  and  cp  is  the  sample  heat  capacity.  Through  algebraic 

manipulation  of  Eq.  (4)  the  magnitude  and  phase  shift  of  the  harmonic  temperature  response  are 
respectively 


1^1-  A 

1 

a/  2  aco  n 

_P2(M  -  N)2  +  Q2(M  +  N)2  _ 

(5) 


and 


O  =  tan 


P(M  -N)  +  Q(M  +  N ) 
D(M  -N)-  Q(M  +  N) 


(6) 


with  M  =  exp {d-jtoirja},  N  =  exp dyfamja^,  P  =  cos (dy[co7rJa^,  and  Q  =  sin [dyfamja). 


Equations  (5)  and  (6)  can  be  used  to  derive  the  thermal  diffusivity  by  fitting  the  measured 
magnitude  and  phase  shift  of  the  detected  thermal  emission  as  a  function  of  laser  modulation 
frequency.  Use  of  Eq.  (5)  requires  two  fitting  parameters,  A  and  a  ,  while  Eq.  (6)  requires  only 
a  .  In  practice  we  find  phase  shift  fitting  to  be  more  reliable  and  reproducible  than  amplitude 
fitting  because  the  amplitude  of  the  harmonic  temperature  response  becomes  very  small  at  high 
heating  frequencies  and  is  more  susceptible  to  drifts  in  laser  power  and  convective  heat  losses. 
Therefore,  we  report  thermal  diffusivity  values  derived  by  phase  shift  fitting.  Based  on  the 
repeatability  of  measurements  and  fitting  uncertainties,  we  estimate  that  reported  thermal 
diffusivity  values  have  uncertainties  of  about  5%.  A  standard  graphite  material,  IG-1 10  nuclear 
grade  from  Toyo  Tanso  USA,  was  used  to  verify  the  performance  of  our  photo  thermal 
radiometry  setup.  Disc  samples  with  thicknesses  of  1  and  2  mm  were  used.  The  measured 
thermal  diffusivity  at  about  299  K  was  1.02  ±  0.03  cm  s'  ,  in  very  good  agreement  with  the 
literature  value  of  1.03  cm2  s'1  at  the  same  temperature.16"17 

(4)  Electrical  Property  Measurements 

Figure  2  shows  the  experimental  arrangements  used  for  the  electrical  resistivity  and  the 
Hall  coefficient  measurements.  Both  measurements  use  a  standard  four-terminal  van  der  Pauw 
geometry18'19  in  which  a  current  is  driven  between  two  terminals  while  a  voltage  is  measured 
between  the  remaining  two  terminals.  A  DC  power  supply  (Global  Specialties  Instruments, 
model  1302B)  was  used  as  a  constant  current  source  and  a  multimeter  (Keithley,  model  2000) 
was  used  to  measure  voltages.  Silver  paint  was  used  to  form  the  electrical  contacts  between  the 
lead  wires  and  the  disk  samples.  The  Hall  mobility  measurements  were  performed  in  a  custom- 
built  electromagnet  with  a  magnetic  induction  of  0.50  Tesla  (as  measured  using  a  Bell  600 
Gaussmeter,  F.  W.  Bell  Inc.).  The  electrical  resistivity  was  measured  as  a  function  of 
temperature  by  placing  the  specimen  on  an  electrically-isolated  heating  stage.  Sample 
temperatures  were  measured  with  a  type  K  thermocouple. 
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Resistivity  Geometry  Hall  Geometry 


Fig.  2.  Experimental  geometries  for  electrical  property  measurements; 

the  arrows  indicate  the  direction  of  current  flow  for  a  particular 
test. 


The  sheet  resistance  for  a  van  der  Pauw  measurement  is  Rijkl  =  Vkl  / /y ,  where  the  first 

two  indices  denote  the  positive  (i)  and  negative  (/)  leads  of  the  current  input  and  the  second  two 
indices  denote  the  positive  (k)  and  negative  (/)  leads  of  the  voltmeter.  Figure  2  shows  the  van 

der  Pauw  geometry  for  measuring  ^2143-  In  total,  twenty-four  measurements  were  made  on  each 
sample  at  each  temperature.  These  consisted  of  three  different  current  levels  in  the  range  of 
1-2  A  for  each  of  the  permutations  ijkl  =  1234,  3412,  2143,  4321,  2341,  4123,  3214,  and  1432. 
The  three  sheet  resistances  for  each  permutation  were  averaged  and  the  electrical  resistivity,  r, 
was  obtained  by  numerically  solving  the  relation18'19 


exp 


7T\ i(i?1234  +  7?3412  +-^2143  ~*~-^432l) 

4  r 


\  f 

+  exp 
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4  r 
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The  performance  of  our  setup  was  verified  by  making  electrical  resistivity  measurements 
on  thin  sheets  of  nickel  (0.125mm,  Aldrich)  and  Constantan®  (0.25  mm,  Goodfellow).  Our 
measured  room  temperature  resistivities  for  nickel  (7.18  pQ-cm  )  and  Constantan®  (50.4 
pQ  •  cm)  are  within  a  few  percent  of  literature  values.20"21 

The  Hall  coefficient  for  a  particular  configuration  is  defined  by  R+Jjjkl  -Vkld  /  /  y  B  ,  where 

the  superscript  (  +  or  -  )  defines  the  direction  of  the  magnetic  field  perpendicular  to  the  sample 
disk;  the  Hall  geometry  3 142  is  shown  in  Fig.  2.  We  fixed  the  current  input  and  measured  the 

voltage  with  the  magnetic  field  turned  off  and  on.  In  total,  eight  measurements  were  performed 
at  a  fixed  current  level  (between  1  and  2  A)  on  each  specimen:  the  permutations  ijkl  =  1324, 
3142,  2431  and  4213  for  each  of  the  two  magnetic  field  directions  (+0.5  T  and  -0.5  T).  The  Hall 
coefficient  was  then  computed  as  the  average 


Rh  —  (7^1324  +  -^3142 


+  7^243 1 


+  -^4213  +  ^1324  +  -^3142 


+  7?243l  +7?4213)/8 


(8) 


The  bulk  carrier  density  and  Hall  mobility  where  computed  as  ne  =  1/ qeRH  and  jue  =  \Ru\l r , 
respectively,  where  qe  =  -1.602xl0“19  C  is  the  charge  of  an  electron. 
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III.  Experimental  Results 


(1)  Microstructure  and  Composition 

Theoretical  densities  were  computed  from  the  target  volume  fractions  and  the  following 
pure  component  densities  (in  g  cm'3):  pHfBi  =11.212  4,  pZrB  =6.085  4,  puf  =13.10  22, 

pB  =2.34  22,  pIr  =22.65  22,  and  pSiC  =3.214  23.  The  experimental  and  theoretical  densities  are 

listed  in  the  fourth  and  fifth  columns  of  Table  I.  The  sixth  column  of  Table  I  shows  that  most 
specimens  reached  at  least  98%  of  their  theoretical  density;  the  two  exceptions  are  A2  (pure 
HfB2)  and  A3  (pure  ZrB2),  both  of  which  are  around  92%  theoretical  density.  Two  samples,  B3 
(HfB2.i)  and  C2  (HfB2-5%SiC)  have  densities  higher  than  their  theoretical  values.  The  measured 
density  of  B3  is  1%  higher  than  the  theoretical  value,  suggesting  that  the  final  boron  content  did 
not  reach  the  targeted  2.5%  excess  (or  HfB2.i  composition).  The  measured  density  of  C2  is  ~2% 
greater  than  the  theoretical  value,  indicating  that  the  actual  SiC  content  is  lower  than  the  target 
5%  volume  fraction. 

The  sample  analysis  by  XRD  shows  the  characteristic  sharp  diffraction  peaks  of 
crystalline  ZrB2  and  HfB2,  as  well  as  a  variety  of  weaker  peaks  associated  with  impurities  and 
minor  constituents,  as  shown  in  Fig.  3.  Diffraction  peaks  associated  with  Hf02  are  found  for 
most  HfB2-based  specimens.  Diffraction  peaks  associated  with  HfC  can  be  clearly  identified  for 
A2  (HfB2),  Cl  (HfB2-5%SiC),  and  C2  (HfB2-5%SiC)  and  ZrC  peaks  can  be  identified  for  A3 
(ZrB2).  No  XRD  peaks  related  to  Ir  (Bl:  HfB2-l%Ir)  or  Hf  (B2:  HfB2-5%Hf)  were  found. 
Diffraction  peaks  associated  with  SiC  were  observed  for  sample  C3  (ZrB2-20%SiC),  but  could 
not  be  clearly  identified  for  samples  Cl  and  C2  containing  nominally  5%  SiC. 

Consistent  with  XRD  results,  SEM/EDS  specimen  characterization  confirmed  the 
presence  of  oxygen  and  carbon  containing  grains  in  most  materials,  and  could  not  identify  any  Ir- 
containing  grains  or  pure  Hf  grains  in  samples  B 1  and  B2,  respectively.  A  small  population  of 
boron  particles  was  found  in  specimen  C3  (ZrB2-20%SiC),  which  was  likely  inadvertently 
introduced  during  the  manufacturing  process.  SiC  grains  were  sparse  but  clearly  identifiable  in 
sample  Cl.  However,  no  SiC  grains  were  found  in  sample  C2.  Instead,  a  variety  of  distinct  HfC 
and  carbon-rich  grains  were  identified.  This  surprising  finding  is  consistent  with  the  presence  of 
the  strong  HfC  lines  in  the  XRD  spectrum  (Fig.  3),  as  well  as  the  aforementioned  discrepancy 
between  the  measured  density  and  theoretical  density  for  this  specimen.  Specimen  C2  was 
manufactured  using  elemental  Hf  and  B  powders,  and  the  reaction  Hf  +SiC  — ►  Si  +  HfC,  which 
could  lead  to  the  loss  of  SiC,  is  thermodynamically  favored.24  However  EDS  was  unable  to 
detect  any  Si-containing  grains  in  sample  C2. 
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f  A1:  HfBj  B1:Hfflj-5%lr  Cl:  HfB2-5%SiC  \  f  *HJ02  0  J&C  } 

|  A2:  HfBj  82:^,  *  C2:  HfB^5%SiC  !  j 

{  A3:  ZrBj  B3:  HtB, ,  C3:  ZrB,«20%SlC  j  \  vHfC  ♦  SIC  1 


Fig.  3.  X-ray  diffraction  scattering  intensity  versus  scanning  angle  20;  intensities  have  been 
normalized  to  the  <10 1>  peak  intensity. 


Grain  size  analyses  are  presented  in  Fig.  4  as  histograms  with  superimposed  log-normal 
curve  fits.  The  log-normal  function  provides  a  reasonable  approximation  to  the  experimental 
histograms.  The  values  of  L  and  <JL  derived  from  these  fits  are  listed  in  Table  I,  along  with  the 
derived  grain  boundary  interface  areas  per  unit  volume.  Values  of  L  fall  into  two  groups:  one 
group  around  5.5  pm  and  another  around  10.5  pm.  The  samples  with  L  =5.5  pm  have  roughly 
twice  the  grain  boundary  interface  area  per  volume  than  those  with  L  =  10.5  pm.  The  grain 
composition  was  not  considered  during  grain  size  analysis  of  samples  Cl  and  C2.  For  sample 
C3,  we  have  differentiated  the  two  grain  boundary  interfaces  of  ZrB2-ZrB2  and  ZrB2-SiC  during 
analysis,  and  the  individual  and  total  interfacial  areas  of  C3  are  listed  in  Table  I. 
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Fig.  4.  Grain  size  distributions  obtained  from  the  analysis  of  SEM  images;  the  mean  and 
standard  deviations  of  the  intercept  lengths  are  indicated  in  each  panel. 


(2)  Thermal  Properties 

The  measured  thermal  diffusivities  are  plotted  in  Fig.  5  over  the  298-700  K  temperature 
range.  The  thermal  diffusivity  of  all  materials  decreases  with  increasing  temperature.  The  Ir- 
containing  HfB2  composite  B1  has  the  highest  thermal  diffusivity  and  the  pure  HfB2  material  A2 
has  the  lowest  thermal  diffusivity  over  this  temperature  range;  this  difference  is  about  a  factor  of 
1.5  at  room  temperature  and  1.2  at  700  K. 


TEMPERATURE.  K 


Fig.  5.  Thermal  diffusivities  as  a  function  of  temperature. 
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Figure  6  shows  thermal  conductivity  values  computed  from  the  measured  thermal 
diffusivity  data  using  the  relation  X  -  apcp .  We  have  used  the  measured  room-temperature 

densities  of  each  sample,  considering  thermal  expansion  over  the  298-700  K  temperature  range 
insignificant.  Specific  heat  values  for  pure  components  were  computed  as  a  function  of 
temperature  using  the  expression 

cp(T)=  A  +  BT  +  CT~2  +DT2  (9) 

with  coefficients  derived  from  the  HSC  Chemistry  database  24;  these  coefficients  are  listed  in 
Table  II.  For  composite  samples,  specific  heat  values  were  mass-averaged  using  theoretical 
densities  and  nominal  volume  fractions.  We  estimate  that  the  uncertainty  in  our  computed 
specific  heat  values  does  not  exceed  10%.  The  thermal  conductivity  of  all  materials  also 
decreases  over  this  temperature  range,  although  not  as  strongly  as  the  thermal  diffusivity. 


Fig.  6.  Thermal  conductivities  calculated  from  the  measured  thermal  diffusivities. 


Table  II.  Specific  heat  function  parameters  of  Equation  (9).24 


A 

B 

C 

D 

HfB2 

366.52 

0.03910 

-1.150X107 

0 

ZrB2  (<600K) 

-34.70 

2.068 

-4.830xl05 

-1.737X10'3 

ZrB2  (>600K) 

583.45 

0.06455 

-1.417x  107 

5.131xl0'6 

Hf 

127.23 

0.04645 

2.527xl05 

1.681xl0'9 

B 

1483.16 

1.193 

-7.003xl07 

-2.992X10'4 

Ir 

118.58 

0.03218 

2.091xl05 

-l.OxlO'9 

SiC 

824.55 

0.5819 

-2.786xl07 

-1.884X10'4 
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(4)  Electrical  Properties 


Figure  7  shows  the  results  of  electrical  resistivity  measurements.  Electrical  resistivity 
increases  linearly  with  increasing  temperature  for  all  samples.  The  data  for  each  sample  were  fit 
by  the  equation 

r  =  r29S[\  +  a(T- 298)]  .  (10) 

The  resulting  values  of  the  room  temperature  resistivity  r298  and  average  temperature  coefficient 
of  resistivity,  a  ,  are  listed  in  Table  III.  The  values  of  r298  and  a  for  eight  of  the  nine  samples 

fall  within  the  ranges  6-9  pQ  •  cm  and  3.5-4.3x10'  K'  ,  respectively.  The  only  exception  is 
sample  C3  which  has  a  considerably  higher  room  temperature  resistivity  (17.0  pQ  •  cm )  and 
lower  average  temperature  coefficient  of  resistivity  (2.12x10'  K' )  than  all  the  other  samples 
because  of  its  high  SiC  content.  The  electrical  resistivity  of  SiC  is  very  high  compared  to  the 

9  S 

diborides,  on  the  order  of  10  Q-cm  for  undoped  material. 


Fig.  7.  Temperature  dependence  of  the  electrical  resistivity;  lines  represent  linear  fits  using 
Eq.  (10). 


The  final  three  columns  of  Table  III  list  the  measured  Hall  coefficient  and  the  derived 
charge  carrier  density  and  Hall  mobility  for  each  sample.  The  measured  Hall  coefficients  are  all 
negative  confirming  that  electron  transport  dominates  charge  flow  in  these  materials.  All  Hall 
coefficients  fall  between  -1.3xl0'3  and  -1.9xl0'3  cm3C_1,  with  C2  having  the  lowest  value  and  C3 
the  highest  value.  The  carrier  densities  are  in  the  range  3.3-4.8xl021  cm'3  and  Hall  mobilities  in 
the  range  100-220  cm2V'1s'1.  Samples  A1  and  Cl  have  the  highest  Hall  mobility.  Sample  C2 
has  the  highest  carrier  concentration,  while  C3  has  the  lowest  values  of  both  Hall  mobility  and 
carrier  concentration. 
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Table  III.  Electrical  properties:  resistivity  at  298  K  (r298),  temperature  coefficient  of 
resistivity  (a),  Hall  coefficient  ( RH ),  bulk  carrier  density  ( ne ),  and  Hall  mobility  ( jue ). 


Group 

ID 

r298 

(|uf2-cm) 

a 

(xl0‘3  K'1) 

Rh 

(x  1 0"3  cm3  C1) 

ne 

(xlO21  cm"3) 

/c 

(cm2  V"1  s'1) 

A1 

6.32 

4.18 

-1.60 

3.9 

220 

A2 

8.90 

3.50 

-1.64 

3.8 

170 

A3 

6.71 

4.34 

-1.55 

4.0 

200 

B1 

6.10 

4.21 

-1.76 

3.5 

250 

B2 

7.07 

4.03 

-1.63 

3.8 

200 

B3 

6.28 

4.10 

-1.46 

4.3 

210 

Cl 

6.37 

3.68 

-1.55 

4.0 

220 

C2 

6.63 

3.60 

-1.30 

4.8 

180 

C3 

17.0 

2.12 

-1.89 

3.3 

100 

(5)  Electronic  Contributions  to  Thermal  Transport 

The  low  electrical  resistivity  and  the  relatively  high  carrier  density  and  Hall  mobility  of 
these  diboride-based  materials  imply  that  electronic  transport  contributes  substantially  to  heat 
transport.  If  the  total  thermal  conductivity  is  taken  as  the  sum  of  electronic  and  phonon 
contributions,  A  =  Ae+Aph,  and  the  electronic  contribution  is  approximated  by  the  Wiedmann- 

Franz  relationship,  it  is  possible  to  estimate  the  contribution  of  Ae  and  Aph  from  the  following: 

A,  (: T )  =  2  (T)-Ae(T)  =  A  (T) - f — ^ -  (11) 

Ph\  1  expV  )  e\  )  expV  )  [l  +  (J  -  298)] 

The  results  of  this  estimation  are  shown  in  Fig.  8  where  the  ideal  Lorenz  Number, 
L0  =2.45xl0“8  W  Q  K~2,  is  used.  Table  IV  list  the  room-temperature  thermal  and  electrical 
conductivities  along  with  the  room-temperature  Wiedemann-Franz  ratio  (Aej  Aexp  =29^oL()  j  Aexp) 

which  is  higher  than  0.8  for  all  materials  except  C3.  Because  Eq.  (1 1)  is  based  on  a  number  of 
physical  simplifications,  and  additionally  propagates  the  experimental  uncertainties  associated 
with  both  our  thermal  and  electrical  property  measurements,  the  absolute  numerical  values  of 
2p/i(r)  have  a  large  uncertainty.  Nevertheless,  it  is  evident  from  Table  IV  and  Fig.  8,  that  in 

300-700  K  temperature  range,  electronic  contributions  to  the  thermal  conductivity  dominate  over 
phonon  contributions. 
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Fig.  8.  Total  (solid  square),  electronic  (open  circle),  and  phonon  (solid  triangle)  thermal 
conductivities  versus  temperature;  the  electronic  and  phonon  contributions  are 
calculated  from  Eq.  (11). 

Table  IV.  Values  of  thermal  conductivity  (extrapolated),  electrical  conductivity,  and 
the  Wiedemann-Franz  ratio  at  298  K. 


Group 

ID 

^298 

(W  m'1  K'1) 

°298 

(xlO7  S  m"1) 

\  /  \xp 

A1 

112 

1.58 

1.0 

A2 

89 

1.12 

0.92 

A3 

108 

1.48 

1.0 

B1 

138 

1.62 

0.86 

B2 

103 

1.40 

0.99 

B3 

125 

1.58 

0.92 

Cl 

137 

1.55 

0.83 

C2 

118 

1.50 

0.93 

C3 

87 

0.60 

0.50 

A2-14 


IV.  Discussion 


The  measured  thermal  conductivities  are  in  the  upper  range  of  values  reported  for 
polycrystalline  materials  with  similar  composition,  both  pure  diborides  and  their  mixtures  with 
SiC.26"28  The  effective  (measured)  thermal  and  electrical  conductivities  of  polycrystalline 
ceramics  depend  on  the  individual  constituent  conductivities  and  the  additional  resistance  to 
transport  due  to  porosity  and  grain  boundaries.  Insight  into  these  contributions  can  be  gained 
through  simple  analytic  models:  the  Brick  Layer  Model  (BLM)  for  interface  resistance  and  the 
Effective  Medium  Approximation  (EMA)30  for  averaging  constituent  properties.  The  BML 
formula  approximates  the  effect  of  interface  resistance  on  thermal  conductivity  by 


f 


A  = 


1  R 


\ 


+ 

v4,  a  j 


(12) 


where  Rx  is  the  thermal  boundary  resistance,  a  is  the  brick  (grain)  length,  and  AuU  is  the 

intrinsic  thermal  conductivity.  The  EMA  approach  mixes  constituent  conductivities  according  to 
the  implicit  relationship 


At  /Le# 

t  '  A  +  2  Aeff 


=  0  , 


(13) 


where  v;.  is  the  constituent  volume  fraction.  Analogous  equations  apply  for  electrical 
conductivity. 

For  a  composite  in  which  porosity  is  isolated  in  a  polycrystalline  matrix  of  one  major 
constituent,  the  two  analytic  models  can  be  combined,  with  the  BLM  used  to  model  the  effect  of 
grain  boundary  resistance  on  the  matrix  conductivity  and  the  EMA  model  used  to  incorporate  the 
effect  of  porosity.30"32  The  “intrinsic”  thermal  and  electrical  conductivities  of  ZrB2  and  HfB2  are 
not  known,  but  based  on  the  measurements  of  Kinoshita  et  al.  for  single  crystal  ZrB2, 
reasonably  representative  values  are  Amt  =  140  W  m"1  K"1  and  crint  =  2.17  xlO7  S  m'1.  With  these 

values  fixed,  setting  a  =  L  ,  and  the  porosity  s  - 1  -  pexpj  pthe ,  the  interfacial  resistances  can  be 

adjusted  until  the  computed  effective  conductivity  matches  the  measured  conductivity  for  the 
nominally  pure  diboride  materials. 

For  the  HfB2  materials  Al  and  A2,  the  experimental  room-temperature  conductivities 
(Table  IV)  are  matched  with  thermal  interface  resistances  of  1.63xl0~8  and  1.50xl0~8  m2  K  W"1, 
and  electrical  interface  resistances  of  1.65x10"  and  1.77x10"  m  S,  respectively.  Since  these 
interface  resistances  vary  little  between  Al  and  A2,  the  reduced  thermal  and  electrical 
conductivities  of  A2  can  be  satisfactorily  explained  by  its  higher  porosity  and  higher  grain 
boundary  area  per  volume.  A  similar  analysis  of  A3  gives  lower  thermal  and  electrical  interface 
resistances  of  0.60x10"  m  K  W"  and  0.76x10"  m  S.  The  influence  of  grain  size  distribution 
on  effective  conductivity  can  be  incorporated  into  the  bricklayer  model  by  setting 
a  =  L  exp(2.5(cr/Z.)2).34  The  factor  a S L  lies  within  the  range  0.4-0. 6  for  all  our  materials; 
making  this  substitution  for  a  increases  the  derived  interface  resistances  by  factors  of  1.67,  2.32, 
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and  2.25,  for  Al,  A2,  and  A3,  respectively.  Typical  thermal  grain  boundary  resistances  in 
polycrystalline  ceramics  fall  in  the  range  10"7  to  10~9  m2  K  W'1 29;  35'36. 

BLM/EMA  analyses  can  also  be  performed  for  materials  with  additional  constituents 
(e.g.,  Ir,  SiC,  excess  Hf,  excess  B).  Although  “intrinsic”  property  inputs  multiply,  such  modeling 
nevertheless  suggests  that  the  high  thermal  conductivities  of  Cl  (HfB2-5%SiC)  and  Bl(HfB2- 
1  %Ir)  are  better  explained  by  lower  thermal  boundary  resistances  than  as  the  result  of  high 
thermal  conductivity  second  phases. 

It  is  interesting  to  compare  samples  B2,  Al,  and  B3,  which  were  manufactured  from 
elemental  powders  mixed  to  the  stoichiometries  HfBi.9,  HfB2,  and  HfB2.i,  respectively.  All  three 
samples  have  similar  grain  sizes  and  grain  boundary  area  per  volume.  Sample  B2  has  the  lowest 
thermal  and  electrical  conductivities  of  the  group,  while  B3  has  a  highest.  However,  enhanced 
electronic  heat  transport  alone  does  not  account  for  the  entire  difference  in  thermal  conductivity 
between  samples  B2  and  B3;  the  phonon  contribution  to  thermal  conductivity  is  also  augmented 
in  the  higher  boron  content  material  (see  Fig.  8  and  Table  IV.)  A  similar  observation  can  be 
made  for  samples  Cl  and  Bl;  their  high  thermal  conductivities  seem  to  derive  in  part  from 
additional  phonon  contributions  to  heat  transport,  not  from  enhanced  electronic  contributions 
(Fig.  8  and  Table  IV.).  The  addition  of  minor  constituents  to  diboride  materials  can  exert  a  large 
influence  on  sintering  behavior  and  chemistry  at  grain  boundaries.  We  speculate  that  such 
additions,  and  the  interface  modifications  they  introduce,  may  influence  the  efficiency  of 
electrical  and  phonon  heat  transport  across  grains  boundaries  in  different  ways. 

Table  V  summarizes  some  published  electrical  property  data  for  HfB2-  and  ZrB2-based 
materials.  Similar  to  other  investigators,  we  find  high  carrier  mobilities  for  the  diborides,  with 
slightly  lower  Hall  coefficients  and  slightly  higher  carrier  densities  for  our  SPS  UHTC  materials 
than  typical  in  Table  V.  Our  measured  electrical  resistivities  are  lower  than  most  literature 
values  for  polycrystalline  diborides  manufactured  by  conventional  hot-pressing  methods.  As 
shown  in  Fig.  8  and  Table  IV,  a  Wiedemann-Franz  analysis  for  these  SPS  materials  indicates  that 
more  than  80%  of  their  total  room-temperature  thermal  conductivity  can  be  attributed  to 
electronic  contributions,  except  for  sample  C3  (ZrB2-20%SiC)  in  which  the  contribution  is  50%. 
Similar  analyses  presented  in  the  literature  vary  widely.  Zimmermann  et  al.  find  electronic 
contributions  of  50-60%  for  ZrB2  and  ZrB2-30%SiC;  Samsonov  et  al.37  report  -65%  for  ZrB2 
and  HfB2;  Tye  and  Clougherty  derive  contributions  ranging  from  40  to  90%  for  a  variety  of 
HfB2  and  ZrB2  materials. 

The  Wiedemann-Franz  relationship  also  predicts  an  increasing  electronic  contribution  to 
thermal  conductivity  with  increasing  temperature  if  the  temperature  coefficient  of  resistivity  is 
below  a  <0.003356  K'1  (i.e.,  <1/298  K;  see  Eq.  11.)  We  measure  a  <0.003356  K'1  for  all  our 
SPS  samples  except  C3,  whereas  a  <  0.003356  K'1  for  most  materials  in  Table  V.  Both 
increasing  [dA/dT  >  0)  and  decreasing  [dA/dT  <  0)  thermal  conductivities  with  increasing 

temperature  are  reported  in  the  literature  for  ZrB2-  and  HfB2-based  materials.  '  ’  '  We  find 
dA/dT  <0for  all  of  our  SPS  UHTC  materials.  Researchers  who  have  measured  dA/dT  >0  and 
have  also  made  temperature-dependent  resistivity  measurements  report  a  <  0.003356  .  27;  37 
However,  others  have  found  a  <0.003356  together  with  dA/dT  <0  28,  as  in  our  data  for  sample 
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C3.  Since  the  phonon  conductivity  always  decreases  with  increasing  temperature,  an  additional 
mechanism  is  required  to  explain  dA/dT>0  .  Gasch  et  al.  have  examined  the  closing  of 

microcracks  with  increasing  temperature  as  a  possibility,  but  have  found  that  an  unsatisfactory 
explanation  for  various  polycrystalline  HfB2  materials.26  It  seems  that  a  <  0.003356  may  be  a 

necessary  but  insufficient  condition  for  a  diboride  to  exhibit  increasing  thermal  conductivity  with 
increasing  temperature. 

Table  V.  Published  electrical  property  measurements  for  HfB2-  and  ZrB2-based  materials: 
resistivity  at  298  K  ( r298 ),  temperature  coefficient  of  resistivity  ( a  ),  Hall  coefficient  ( RH ), 
bulk  carrier  density  (  ne ),  and  Hall  mobility  (jue ). 


Material 

r298 

(|uf2*cm) 

a 

(x10"3K_1) 

Rh 

(xlO'3  cm3  C'1) 

ne* 

(xlO21  cm'3) 

Re* 

(cm2  V"1  s"1) 

Ref. 

(Polycrystalline) 

HfB2 

15.8 

-1.7 

3.7 

110 

41 

HfB2 

16. 6f 

— 

-1.70 

3.7 

102 

42 

HfB2 

10.4 

3.29 

-1.80 

3.5 

170 

43 

HfB2-5%SiC 

12.4 

3.37 

— 

— 

— 

28 

ZrB2 

7 

— 

-2.0 

3.1 

290 

41 

ZrB2 

16.6f 

— 

-1.76 

3.6 

106 

42 

ZrB2 

7.8 

1.3 

— 

— 

— 

44 

ZrB2 

9.6 

2.21 

-1.90 

3.3 

200 

43 

ZrB2 

11.9 

2.72 

— 

— 

— 

28 

ZrB2  (90%  dense) 

11.0 

3.73 

— 

— 

— 

28 

ZrB2 

22 

2.06 

— 

— 

— 

27 

ZrB2-20%SiC 

10.3 

4.42 

— 

— 

— 

28 

ZrB2-30%SiC 

24 

2.52 

— 

— 

— 

27 

(Single  Crystal) 

ZrB2 

2.9-3.2| 

-1.2",  -2.31 

45 

ZrB2 

4.6f 

— 

— 

— 

— 

33 

*  Computed  from  the  tabulated  r298  and  RH  values  assuming  electrons  are  the  charge  carriers, 
f  Room  temperature  value,  not  from  fitting. 

" 1  parallel  and  perpendicular  to  hexagonal  crystal  axis 
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V.  Conclusion 

We  investigated  the  thermal  and  electrical  transport  properties  of  spark-plasma- sintered 
HfB2-  and  ZrB2-based  polycrystalline  ceramics  with  different  microstructures  and  minor 
additives,  in  the  300-700  K  temperature  range.  The  materials  have  high  carrier  mobilities,  with 
somewhat  lower  Hall  coefficients  and  higher  carrier  densities  than  values  typically  reported  for 
diboride-based  polycrystalline  materials.  The  measured  electrical  conductivities  are  higher  than 
most  literature  values  for  diboride  ceramics  manufactured  by  conventional  hot-pressing  methods. 
The  thermal  conductivities  are  in  the  upper  range  of  values  typically  reported  for  pure  diboride 
and  diboride-SiC  compositions  and  all  show  decreasing  conductivity  with  increasing 
temperature.  A  Wiedemann-Franz  analysis  indicates  that  thermal  conductivity  is  dominated  by 
the  electronic  contribution  to  thermal  transport.  The  variations  in  thermal  and  electrical 
conductivities  between  different  HfB2  samples  are  explained  by  their  different  porosities  and 
grain  boundary  areas  per  volume.  The  high  thermal  conductivities  of  HfB2-5%SiC  and  HfB2- 
l%Ir  samples  are  attributed  to  their  low  thermal  boundary  resistances.  In  samples  of  HfBi.9, 
HfB2,  and  HfB2.i,  which  have  very  similar  microstructures,  both  thermal  and  electrical 
conductivities  increase  with  boron  content,  a  finding  that  deserves  further  investigation. 
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Abstract 

Ultra-high  temperature  ceramic  materials  and  composites  under  development  as  nose-tip 
and  wing  leading  edge  components  for  hypersonic  flight  vehicles  must  operate  in  extreme 
aerothermal  heating  environments.  The  performance  of  ultra-high  temperature  ceramics  for  this 
application  is  ultimately  evaluated  using  high-enthalpy,  long  duration  flow  facilities  that  simulate 
the  reactive  gas  environment  encountered  in  hypersonic  flight.  In  this  paper,  we  describe  the  test 
environments  generated  by  two  types  of  these  ground  test  facilities  -  subsonic  inductively- 
coupled  plasma  tunnels  and  supersonic  arc-jet  tunnels  -  and  discuss  the  important  roles  of 
computational  fluid  dynamics  modeling  and  in  situ  optical  diagnostics  for  interpreting  test  results 
from  a  materials  science  perspective. 


1.  Introduction 

Transition  metal  borides  and  carbides,  as  well  as  their  composites  with  various  sintering 
aids  and  glass  formers,  are  collectively  known  as  ultra  high  temperature  ceramics  (UHTCs).  The 
field  of  UHTC  research  has  expanded  substantially  over  the  last  decade,  motivated  by  the  unique 
potential  of  these  ceramics  for  applications  in  extreme  high-temperature,  reactive 
environments.1,2  A  central  driver  for  UHTC  research  is  the  need  for  sharp  leading  edge  and 
control  surface  components  for  future  generations  of  hypersonic  flight  vehicles.  '  Sharp 
leading  edges  enable  vehicles  with  flight  performance  at  hypersonic  speeds  (e.g., 
maneuverability,  extended  cross-range  capability)  that  cannot  be  achieved  with  blunt  body 
designs.6,7  The  main  hurdle  for  sharp  vehicle  designs  is  the  severe  aerothermal  heating 
environment  that  acts  on  sharp  leading  edges.  The  concentrated  combination  of  high 
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temperature  and  aggressive  chemistry  causes  most  materials  to  fail  by  melting, 
vaporization/sublimation,  oxidation,  ablation,  spalling,  or  some  combination  of  these  processes. 

While  a  battery  of  conventional  mechanical  and  thermal  tests  aid  in  the  development  of 
UHTC  materials,  the  performance  of  UHTCs  intended  for  leading  edge  applications  must 
ultimately  be  demonstrated  in  a  representative  aerothermal  test  environment.  Such  test 
environments  are  provided  by  arc-jet  or  inductively-coupled  plasma  (ICP)  wind  tunnels  which 
expose  test  articles  to  high-enthalpy  reactive  gas  flows.  Tests  of  UHTC  materials  in  arc -jet 
facilities  have  been  reported  in  the  scientific  literature  by  Kaufman3,  Metcalfe  et  al.8,  Wuchina 
and  Opeka9,  Opila  et  al.10,  Gasch  et  al.11,  Chamberlain  et  al.12,  Savino  et  al.13,  Monteverde  and 
Savino14,  and  Zhang  et  al.15  Tests  of  UHTC  materials  in  ICP  facilities  have  been  reported  by  Ito 
et  al.16,  Marschall  et  al.17,  and  Playez  et  al.18 

In  this  paper  we  discuss  the  test  environments  provided  by  these  two  types  of  plasma 
wind  tunnel  facilities.  We  highlight  the  coupling  between  the  high-enthalpy  reactive  flow  stream 
and  the  test  specimen  in  determining  UHTC  material  response  to  these  environments.  The 
importance  of  computational  fluid  dynamics  (CFD)  modeling  of  the  free-stream  and  boundary 
layer  flows  is  discussed,  and  the  value  of  in  situ  diagnostics  for  documenting  test  conditions  and 
monitoring  UHTC  evolution  during  testing  is  emphasized. 


2.  High-Enthalpy  Test  Environments 

2.1.  Plasma  Wind  Tunnel  Facilities 

The  need  for  long-duration,  high-enthalpy  gas  flows  for  the  characterization  and 
qualification  of  thermal-protection  system  (TPS)  materials  and  components  led  to  the 
development  of  plasma  wind  tunnel  facilities  for  aerospace  applications.19  A  variety  of  different 
designs  have  been  investigated  and  implemented  for  plasma  wind  tunnels  but  in  general  all 
facilities  consist  of  an  electrical  power  supply,  an  input  gas  supply,  an  arc-discharge  or 
inductively-coupled  heater  section,  a  test  section,  and  an  exit  gas  handling  system. 

Arc  discharge  devices  in  which  a  current  passes  directly  through  a  flowing  gas  were 
developed  extensively  in  the  USA  because  of  their  large  energy  deposition  rates.  Stable  arc 
discharges  are  achieved  in  large  length-to-diameter  ratio  heaters  constructed  of  segmented  metal 
elements  and  fitted  with  magnetically-spun  electrodes.  Various  arc-jets  are  now  in  operation 
around  the  world,  including  facilities  in  the  USA,  Europe,  and  Asia.19  Direct-arc  heating 
facilities  can  typically  operate  with  large  test  gas  mass  flows  and  with  pressures  on  the  order  of 
1-2  MPa  within  the  heater,  enabling  supersonic  flow  conditions  in  the  test  section.  Most 
constricted-arc  plasma  facilities  use  conical  nozzles  to  produce  supersonic  flows  and  can  operate 
stably  for  long  test  times.  The  United  States  Air  Force  (USAF)  arc -jet  facilities  have  been 
developed  to  simulate  aero-heating  for  ballistic  atmospheric  trajectories,  while  the  National 
Aeronautics  and  Space  Administration  (NASA)  arc-jet  facilities  have  been  tailored  to  simulate 
lower  pressure  trajectories  associated  with  upper  atmospheric  and  planetary  entry.19 

In  ICP  devices  energy  is  coupled  into  the  test  gas  flow  inductively  using  high-frequency, 
high-power,  high-voltage  electrical  supplies.  Gases  flow  through  a  quartz  tube  (typically)  in  the 
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heater  section  and  direct  contact  between  the  gas  stream  and  hot  metal  surfaces  is  avoided.  One 
advantage  of  the  ICP  design  for  materials  science  studies  is  that  the  free-stream  is  usually  free  of 
metal  contaminants  (e.g.  copper)  often  found  in  arc-jet  flows.  Such  contaminates  may  interfere 
with  important  chemical  processes  like  oxidation,  reactive  volatilization  and  surface  catalysis,  or 
change  important  surface  properties  like  the  emittance.  Although  these  devices  were 
investigated  in  the  USA,20  they  were  developed  extensively  in  Russia  and  large-scale  ICP 
facilities  are  now  also  in  operation  in  Europe  and  Asia.19  These  devices  can  generate  either 
subsonic  or  supersonic  flows,  however  because  of  their  more  limited  capability  for  supporting 
high  gas  pressures  in  the  heater  section,  most  ICP  testing  is  done  in  the  subsonic  regime. 

The  test  environments  generated  by  these  two  types  of  facilities  differ  in  significant  ways, 
both  from  each  other  and  from  actual  flight  environments.  No  ground  facility  can  reproduce  all 
features  of  hypersonic  flight  accurately,  because  of  physical  constraints  and  operating  envelope 
limitations.  Choices  for  test  conditions  are  usually  made  to  match  particularly  important 
trajectory  parameters  like  the  anticipated  peak  heat  flux,  peak  pressure,  maximum  heat  load,  etc. 
The  extrapolation  of  ground  test  conditions  to  flight  environments  is  an  on-going  area  of  concern 
and  research21,22  and  will  not  be  addressed  further  here.  Rather,  we  focus  on  differences  between 
supersonic  arc -jet  and  subsonic  ICP  environments  as  they  impact  TPS  materials  testing  in 
general  and  UHTC  testing  in  particular. 


2.2.  Test  Conditions 

Each  plasma  facility  has  its  unique  envelope  of  operation  that  limits  the  range  of  free- 
stream  enthalpies  that  can  be  obtained.  Together  with  the  model  size,  shape,  and  chemical 
composition,  the  available  free-stream  enthalpy  determines  the  heat  flux  that  can  be  applied  to  a 
test  specimen.  Typically,  the  power  delivered  to  the  heater,  the  mass  flow  rate  of  the  test  gas, 
and  the  test  section  static  pressure  can  be  adjusted  within  prescribed  ranges.  The  stagnation 
point  conditions  achieved  for  a  particular  combination  of  facility  settings  are  then  measured 
using  a  calibration  probe  (or  probes)  of  the  same  geometry  as  the  test  sample,  but  fitted  with  a 
calorimeter  to  measure  heat  flux  and  a  Pitot  gauge  to  measure  pressure.  These  calibration  probes 
are  usually  made  of  copper  and  are  water  cooled. 

Figure  1  illustrates  the  salient  features  of  the  flow  interactions  with  a  test  specimen  for 
the  each  type  of  facility.  An  important  difference  between  supersonic  arc-jet  and  subsonic  ICP 
flows  is  the  thermochemical  state  of  the  gas  approaching  the  test  article.  Because  of  the  relative 
fluid  dynamic  and  chemical  relaxation  time  scales  in  arc -jet  facilities,  gases  exiting  the  nozzle 
and  flowing  into  the  test  chamber  are  typically  in  a  chemically  frozen  or  non-equilibrium  state. 
This  means  that  the  chemical  composition  of  the  free-stream  gas  may  differ  considerably  from 
the  thermodynamic  equilibrium  composition  associated  with  the  translational  gas  temperature. 

In  supersonic  arc-jet  flows  the  interaction  of  the  free-stream  with  the  test  article  produces  a  bow 
shock  in  which  the  gas  undergoes  adiabatic  compression.  High  energy  intermolecular  collisions 
heat  the  gas,  depositing  energy  into  a  manifold  of  excited  electronic,  vibrational,  and  rotational 
states,  and  driving  ionization  and  molecular  dissociation  processes.  Behind  the  shock  front, 
where  the  flow  is  subsonic,  the  residence  time  is  now  longer  than  typical  chemical  reaction  times 
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and  gases  undergo  thermal  and  chemical  relaxation  as  they  approach  the  boundary  layer  edge. 
The  gas  temperature  at  the  boundary  layer  edge  may  be  1000’s  of  degrees  hotter  than  the  surface 
of  the  test  specimen.  As  the  gas  flow  passes  through  the  boundary  layer  edge  to  the  surface, 
chemical  reactions  change  the  gas  composition  further.  However,  the  chemical  and  flow  time 
scales  in  the  boundary  layer  may  be  such  that  the  gas  interacting  directly  with  the  specimen  is 
not  in  chemical  equilibrium  at  the  surface  temperature. 
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Fig.  1.  Flow  features  for  arc-jet  and  ICP  free-streams  interacting  with  a  test  specimen. 


In  contrast,  in  ICP  facilities  the  flow  path  and  timescales  are  typically  such  that  gases  can 
relax  towards  their  thermochemical  equilibrium  state  in  the  free-stream  before  reaching  the 
boundary  layer  edge  at  the  test  article,  making  the  prediction  of  the  boundary  layer  edge 
composition  much  more  straight  forward.22  In  a  subsonic  ICP  flow,  no  bow  shock  is  produced 
when  the  flow  encounters  the  test  specimen,  but  a  reactive  boundary  layer  is  still  formed.  There 
is  still  a  large  temperature  difference  between  the  boundary  layer  edge  and  the  specimen,  and 

24 

similar  drivers  for  chemical  non-equilibrium  in  the  gases  at  the  surface  exist. 

In  supersonic  arc-jet  flows  the  pressure  at  the  stagnation  point  of  a  test  article  is  largely 
determined  by  the  dynamic  pressure  of  the  flow  stream,  while  in  subsonic  ICP  flows  the 
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stagnation  pressure  is  largely  determined  by  the  static  pressure  in  the  test  section  because  the 
dynamic  pressure  of  the  flow  stream  is  small.  The  aerodynamic  shear  stresses  exerted  on  model 
surfaces  can  be  much  higher  in  supersonic  arc -jet  flows  than  in  subsonic  ICP  flows.  While 
UHTC  materials  should  have  sufficient  mechanical  strength  to  withstand  aerodynamic  shear 
stresses,  this  may  not  be  true  for  developing  surface  oxides.  At  sufficiently  high  temperatures  in 
oxygen-bearing  environments,  many  UHTC  materials  form  glassy  oxide  phases  (e.g.,  B2O3  and 
SiC>2)  that  permeate  and  seal  porous  poly-crystalline  oxides  (like  ZrC>2  and  Hf02)  creating  a 
composite  scale  that  acts  as  an  efficient  oxygen  diffusion  barrier.25'28  Aerodynamic  shear 
stresses  may  redistribute  low-viscosity  glassy  phases  and  alter  the  formation  of  such  composite 
oxide  scales. 

Survey  probe  heat  flux  measurements  are  often  reported  as  part  of  UHTC  test  conditions, 
but  it  is  important  to  recognize  that  this  heat  flux  is  a  cold  wall  heat  flux  to  a  highly  catalytic 
surface,  which  is  not  the  same  as  the  heat  flux  delivered  to  a  UHTC  test  specimen  under  identical 
test  conditions.  The  convective  heat  flux  to  a  specimen  is  driven  by  the  temperature  difference 
between  the  specimen  surface  and  the  hotter  gas  at  the  boundary  layer  edge.  A  UHTC  surface 
will  reach  very  high  temperatures  during  an  arc-jet  or  ICP  test  because  it  is  not  actively  cooled. 
Because  of  the  smaller  temperature  difference,  the  hot  wall  heat  flux  to  a  hot  sample  surface  is 
always  lower  than  the  cold  wall  heat  flux  measured  by  a  water-cooled  survey  probe.  In  addition, 
a  UHTC  surface  may  be  less  efficient  than  a  copper  surface  at  catalyzing  the  recombination  of 
atomic  species  in  the  flow  and  will  therefore  experience  less  chemical  heating  from  exothermic 
surface  recombination  reactions  like  O  +  O  — ►  O2,  N  +  N  — >  N2,  and  O  +  N  — >  NO. 

The  primary  utility  of  cold  wall  heat  flux  measurements  is  to  provide  experimental 
calibration  data  to  aid  in  the  reconstruction  of  the  free-stream  enthalpy  using  CFD  computations 
(discussed  further  below).  Cold  wall  heat  flux  measurements  are  also  useful  for  confirming  the 
reproducibility  of  test  conditions,  but  caution  must  be  used  when  using  free-stream  enthalpy  or 
cold  wall  heat  flux  to  compare  test  conditions  from  supersonic  arc-jet  tests  with  subsonic  ICP 
tests.  Previous  investigations  have  shown  that  three  parameters  must  be  replicated  to  match 
stagnation  point  convective  heating  conditions  in  ground  test  facilities:  1)  stagnation  point 
enthalpy;  2)  stagnation  (or  impact)  pressure;  and  3)  velocity  gradient  at  the  boundary  layer 
edge.22,29  The  velocity  gradient  is  the  rate  that  the  velocity  increases  as  gas  flows  around  the  test 
article.  Thus,  the  convective  heat  transfer  to  the  same  test  specimen  exposed  to  arc -jet  and  ICP 
flows  with  similar  free-stream  enthalpy  and  total  pressure  values  can  be  substantially  different, 
owing  to  the  large  difference  in  velocity  gradients  between  subsonic  and  supersonic  test 
facilities.  When  differences  in  the  velocity  gradients  in  the  two  types  of  facilities  are  properly 
considered  (for  example,  by  scaling  test  article  dimensions  appropriately)  similar  boundary  layer 
edge  compositions  and  surface  heat  flux  values  can  be  obtained.22,29 

We  note  that  the  characterization  of  high-enthalpy  flow  test  environments  remains  an 
active  research  area  and  that  comparison  of  flow  configurations  and  boundary  layer 
environments  can  also  be  made  on  the  basis  of  non-dimensional  similarity  parameters  such  as  the 
Reynolds,  Schmidt,  and  Damkohler  numbers.  The  use  of  non-dimensional  similarity  parameters 
is  a  valid  approach,  but  one  practical  aspect  worth  mentioning  is  that  (with  some  work)  the 
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stagnation  enthalpy,  stagnation  pressure;  and  velocity  gradient  can  be  measured  more  easily  than 
most  non-dimensional  similarity  parameters  in  arc-jet  and  ICP  test  facilities. 

From  a  materials  science  perspective,  and  particularly  for  studies  of  UHTC  oxidation,  the 
most  relevant  metrics  for  comparing  different  tests  are  the  surface  temperature,  the  stagnation 
point  pressure,  and  the  gas  composition  at  the  sample  surface.  This  information  is  vital  for 
associating  physical  changes  in  test  specimens  with  specific  environmental  conditions,  and  for 
developing  models  of  thermal  and  chemical  material  response.  The  first  two  quantities  are 
routinely  measured  during  arc -jet  and  ICP  tests,  but  the  gas  composition  at  the  surface  -  which 
can  be  influenced  both  by  chemical  and  transport  kinetics  -  is  not  currently  experimentally 
accessible. 


3.  CFD  Modeling 

Since  the  chemical  nature  of  the  gas  at  the  test  specimen  surface  is  not  measured,  a 
combination  of  facility  data,  calibration  probe  measurements  and  CFD  modeling  must  be  used  to 
compute  the  gas  composition  at  the  specimen  surface  for  each  test  run.  This  process  involves 
two  steps.  First,  CFD  model  inputs  are  adjusted  until  ffee-stream  flow  conditions  result  in 
computed  heat  flux  and  pressure  values  consistent  with  those  measured  using  calibration  probes; 
then  the  free-stream  flow  conditions  are  held  fixed  and  the  CFD  model  is  used  to  compute  the 
temperature  drop  and  the  changing  gas  composition  through  the  boundary  layer  to  the  sample 
surface. 

This  procedure  was  described  in  detail  for  a  recent  series  of  UHTC  oxidation  tests  run  in 
the  1.2  MW  Plasmatron  facility  at  the  von  Karman  Institute  for  Fluid  Dynamics  (VKI).17  The 
CFD  codes  used  for  this  procedure  are  the  VKI  Boundary  Layer  Code  24,30  and  the  VKI  ICP 
Code  31’32,  both  of  which  use  the  PEGASE  library  to  perform  thermodynamic  and  transport 
property  calculations.  The  ICP  code  solves  the  time-averaged  magneto-hydrodynamic 
equation  at  low  Mach  and  low  magnetic  Reynolds  numbers,  assuming  axisymmetric  flow  and 
local  thermodynamic  equilibrium,  to  simulate  the  flow  inside  the  plasma  torch  and  around  the 
test  article  in  the  vacuum  chamber.  The  Boundary  Layer  Code  solves  the  boundary  layer 
equations  for  an  axisymmetric  or  two-dimensional,  steady,  laminar  flow  of  chemically  reacting 
gas  over  a  catalytic  surface,  including  thermal  and  chemical  non-equilibrium.  The  temperature 
and  velocity  at  the  boundary  layer  edge  are  adjusted  until  the  computed  and  measured  heat  flux 
agrees. 

At  NASA  Ames  Research  Center,  arc -jet  flows  are  computed  using  the  Data  Parallel  Line 
Relaxation  (DPLR)  code.34  DPLR  is  a  parallel  multiblock  finite-volume  code  that  solves  the 
Navier-Stokes  equations  including  finite -rate  chemistry  and  the  effects  of  thermal  non¬ 
equilibrium.  The  code  is  used  to  compute  nonequilibrium  expanding  flow  in  the  arc-jet  nozzle 

'J  rj  o  n 

and  supersonic  jet  entering  the  test  section,  as  well  as  the  reacting  flow  around  test  articles. 
Thermodynamic  properties  are  taken  from  NASA  Glenn  curve  fits38.  Transport  properties  are 
computed  using  the  self-consistent  effective  binary  diffusion  method  and  expressions  and 
mixing  rules  presented  by  Gupta  et  al.40,  together  with  collision  integrals  compiled  by  Wright  et 
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al.41’42  Simulations  are  typically  started  at  the  nozzle  throat  assuming  thermochemical 
equilibrium  flow  properties  and  some  radial  profile  of  enthalpy  and  mass  flux.  Facility  and 
calibration  data  include  measurements  of  the  chamber  pressure,  mass  flow  rate,  and  test  section 
pressure,  calorimeter  probe  heat  flux  and  pressure  measurements.  The  total  enthalpy  of  the  arc- 
jet  flow  and  its  radial  distribution  are  inferred  from  facility  data  and  previous  survey 
measurements.  In  the  CFD  reconstruction,  these  inferred  values  and  distributions  serve  as  the 
starting  point  for  iterative  adjustments  to  match  computational  predictions  to  the  heat  flux 
measured  by  the  calibration  probe. 

The  energy  balance  at  the  surface  of  a  test  article  or  calibration  probe  is  the  key 
relationship  that  couples  the  gas-phase  and  the  solid  material.  For  an  environment  free  of  gas 
radiation,  the  energy  balance  for  any  non-ablating/non-pyrolyzing  material  surface  can  be  written 
as 


Qchem  drad  dlcond 


(1) 


where  the  terms  on  the  left  hand  side  account  for  convective  and  chemical  heating  by  the  gas 
stream  and  those  on  the  right  hand  side  represent  cooling  by  radiation  and  in-depth  heat 
conduction.  With  some  simplifying  assumptions  (a  sensible  gas  enthalpy  convective  transfer 
coefficient  representation,  independent  surface  recombination  of  O  and  N  atoms  with  equal 
catalytic  efficiencies,  negligible  gas  radiation)  this  equation  can  be  expanded  as 


CH[H,(T,)-H,(T,)]+r'(T,)bE, 


RTS 

8  nM  r 


+/{ts)ae 


N2nN  i 


RT-  =^T;-k(T;dT 


8ttM 


dx 


(2) 


where  CH  is  the  convective  transfer  coefficient;  He  and  Hs  are  the  gas  enthalpies  at  the 
boundary  layer  edge  and  at  the  surface;  A E0  and  A EN  are  molecular  dissociation  energies;  nQ 
and  nN  are  atom  number  densities  above  the  surface;  M0  and  M N  are  molar  masses;  R  is  the 

universal  gas  constant;  T  is  temperature;  cr is  the  Stefan-Boltzmann  constant;  and  x  is  the 
coordinate  into  the  surface.  Equation  (2)  also  shows  how  the  surface  energy  balance  depends 
explicitly  on  temperature-dependant  surface  and  bulk  material  properties:  the  total  catalytic 
efficiency  y' ,  the  emittance  s  ,  and  the  thermal  conductivity,  k  .  The  total  catalytic  efficiency  is 

defined  as  y'  =  yfi ,  where  y  is  the  species  recombination  efficiency  (the  fraction  of  collisions 
with  the  surface  that  result  in  atom  loss)  and  ft  is  the  energy  accommodation  coefficient  (the 
fraction  of  exothermic  reaction  energy  transferred  to  the  surface). 
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In  calibration  measurements  the  right-hand  side  of  Eq.  (2)  is  measured  directly  using 
water-cooled  calorimeters.  Then  CFD  inputs  are  adjusted  to  reproduce  this  calibration  heat  flux 
measurement,  usually  with  values  of  y’  fixed  at  1  (a  fully  catalytic  surface)  and  Ts  fixed  to  a 

low  temperature  like  300  K  (a  cold  wall).  This  process  derives  the  convective  transfer 
coefficient  CH  and  the  boundary  layer  edge  temperature  Te  for  a  particular  test  condition, 

quantities  that  are  then  kept  fixed  in  simulations  seeking  to  reproduce  surface  temperatures 
subsequently  measured  during  UHTC  tests. 

When  qcond  is  negligible,  and  with  CH ,  Te ,  and  Ts  fixed,  Eq.  (2)  shows  that  the  energy 
balance  is  determined  by  y'  and  s  through  the  chemical  and  surface  radiation  fluxes.  (For 

sharp  UHTC  test  articles,  a  thermal  model  of  the  heat  transfer  within  the  UHTC  is  likely  required 
to  evaluate  qcond  .  However  for  other  configurations,  like  stagnation  point  testing  with  disk 

specimens  in  a  flat-faced  cylinder  configuration,  qcond  can  be  set  to  zero  with  minimal  error.)  If 

a  numerical  value  is  assigned  to  the  surface  emittance,  the  total  catalytic  efficiency  can  be 
adjusted  iteratively  until  Eq.  (2)  is  satisfied,  and  the  gas  composition  at  the  surface  is  computed 
simultaneously.  In  practice,  it  is  not  straightforward  to  choose  the  “correct”  value  of  emittance 

1  7 

for  a  UHTC  component  at  high  temperature. 

As  an  illustrative  example  of  the  process  described  above,  we  show  some  results  obtained 
for  a  recent  series  of  stagnation  point  oxidation  tests  performed  in  the  VKI  Plasmatron  on  UHTC 
specimens  of  a  hot-pressed  ZrB2  material  containing  30  volume  percent  SiC  (termed 
ZrB2-30SiC).17  Tests  were  performed  at  static  chamber  pressures  of  104  Pa  with  an  air  mass 
flow  rate  of  16  g  s'1  over  a  range  of  Plasmatron  powers  from  150  to  210  kW. 

The  upper  portion  of  Fig.  2  shows  the  cold  wall  heat  flux  measured  by  a  calibration  probe 
and  the  steady-state  UHTC  sample  surface  temperature  measured  using  a  two-color  pyrometer 
during  each  test  run.  Both  the  cold  wall  heat  flux  and  the  specimen  surface  temperature  increase 
with  Plasmatron  power  as  expected.  The  lower  portion  of  Fig.  2  shows  the  hot  wall  heat  flux  for 
each  test  run  computed  from  qhw  =  soTs4  with  s  =  0.90  and  s  =  0.75 ,  assuming  qcond  ~  0 .  Both 

values  of  emittance  are  reasonable  estimates  for  oxidized  ZrB2-SiC  composites1417’43.  Without 
an  in  situ  method  for  determining  the  high-temperature  emittance  during  testing,  such  a  level  of 
uncertainty  (~15%)  in  emittance  is  unavoidable.  Note  that  the  hot  wall  heat  fluxes  are  only 
about  half  of  the  corresponding  cold  wall  heat  fluxes,  which  suggests  that  the  surface  catalytic 
efficiencies  of  the  oxidized  UHTC  surface  must  be  low. 
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Fig.  2.  Heat  fluxes  and  surface  temperatures  versus  Plasmatron  power.17  The  cold  wall  heat  flux 
is  measured  by  a  calibration  probe  and  the  sample  surface  temperature  is  measured  by  a 
two-color  pyrometer  during  each  test  run.  The  hot  wall  heat  fluxes  are  estimated  from 
qhw  =  soT4  with  s  =  0.90  and  s  =  0.75  ,  assuming  qcond  «  0  . 


Figure  3  plots  the  total  catalytic  efficiency  obtained  for  each  test  condition  from  the 
surface  energy  balance  constraint,  for  values  of  s  =  0.90  and  s  =  0.75  .  The  derived  catalytic 
efficiencies  are  much  less  than  1  consistent  with  a  surface  of  low  catalytic  activity.  Catalytic 
efficiencies  derived  for  s  =  0.90  are  2  to  6  times  higher  than  those  for  s  =  0.75 .  Thus 
uncertainties  in  emittance  become  amplified  in  derived  catalytic  efficiencies. 
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SURFACE  TEMPERATURE,  °C 

Fig.  3.  Recombination  efficiency  versus  surface  temperature  calculated  for  surface  emittances  of 
0.75  and  0.90. 17  The  total  catalytic  efficiencies  for  O  +  O  and  N  +  N  are  assumed 
identical. 


Figure  4  shows  the  calculated  number  densities  of  O,  N,  N2,  O2,  and  NO  at  the  UHTC 
surface  for  emittance  values  of  0.90  and  0.75.  N2,  O  and  N  (in  that  order)  are  the  dominant 
species  with  number  densities  around  10  molecules  per  cubic  meter;  the  number  densities  of  O2 
and  NO  are  about  2  orders  of  magnitude  smaller.  Atomic  oxygen  number  densities  are  relatively 
unaffected  by  the  choice  of  emittance  values,  but  atomic  nitrogen  number  densities  computed  for 
s  =  0.90  are  only  55  to  75%  of  the  magnitudes  computed  for  s  =  0.75  .  The  number  densities  of 
the  minor  species  are  larger  for  s  =  0.90  than  s  =  0.75 ,  by  factors  of  about  3  to  7  for  O2  and 
about  2  to  3.5  for  NO. 
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Fig.  4.  Number  densities  at  specimen  surfaces  versus  surface  temperature,  computed  for 
surface  emittances  of  0.75  and  0.90. 

As  this  example  demonstrates,  considerable  uncertainties  can  propagate  into  CFD- 
derived  quantities  like  catalytic  efficiency  and  surface  gas  composition,  by  errors  in  experimental 
measurements  and  material  properties.  Additional,  and  hard  to  quantify,  error  is  also 
undoubtedly  present  in  the  details  of  the  formulations,  assumptions,  solution  procedures,  and  the 
chemical,  transport  and  thermodynamic  data  employed  in  particular  CFD  codes.  Nevertheless, 
CFD  modeling  is  a  crucial  element  of  UHTC  testing,  because  it  offers  the  best  tool  currently 
available  for  evaluating  the  state  of  the  gas  interacting  directly  with  the  test  surface  specimen.  It 
is  also  the  only  practical  way  to  estimate  quantities  like  the  surface  shear  stress  or  to  map  out 
changes  in  gas  composition  over  the  surface  of  more  complex  3-dimensional  UHTC  test  articles. 

We  note  that  most  published  accounts  of  UHTC  testing  in  arc-jet  or  ICP  facilities  have 
not  reported  quantitative  estimates  of  the  gas  composition  at  the  sample  surface  during  testing. 
From  a  materials  science  perspective,  this  makes  the  comparison  of  different  experiments 
problematic  and  the  construction  of  materials  response  models  more  difficult.  The  value  of 
UHTC  arc-jet  and  ICP  test  results  would  be  greatly  increased  if  surface  gas  compositions,  and 
the  numerical  procedures  used  to  estimate  them,  were  routinely  reported.  Further  benefits  would 
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result  from  measurements  of  boundary  layer  gas  compositions  and  gradients,  as  discussed  in  the 
following  section. 


4.  Optical  Diagnostics 

In  this  section  we  examine  different  optical  techniques  for  obtaining  information  about 
the  test  environment  and  the  evolution  of  UHTC  specimens  exposed  to  high-enthalpy  flows. 


4.1.  Surface  Radiometry 

Surface  temperature  is  perhaps  the  most  important  quantity  required  for  the  interpretation 
of  test  results.  Surface  temperatures  are  typically  measured  in  situ  by  collecting  thermal 
radiation  emitted  by  the  hot  specimen  using  one-color  or  two-color  radiometers.  One-color 
radiometers  collect  radiation  over  a  single  wavelength  range,  while  two-color  radiometers  collect 
radiation  over  two  wavelength  ranges  (which  may  overlap).  Both  types  of  radiometers  are 
calibrated  as  a  function  of  emission  temperature  using  blackbody  radiation  sources  (special 
furnaces  with  an  effective  normal  emittance  approaching  1).  One-color  radiometers  require 
knowledge  of  the  specimen  emittance  over  the  detection  wavelength  range  to  convert  radiation 
intensity  to  surface  temperature.  Two-color  radiometers  only  require  that  the  ratio  of  emittance 
values  for  the  two  detection  ranges  be  known;  for  a  surface  with  wavelength-independent 
emittance  this  ratio  is  1 .  Typical  experimental  issues  include  transmittance  losses  due  to 
windows,  interference  by  absorbing  gases,  parasitic  collection  of  reflected  light,  and  geometric 
limitations  for  viewing  the  sample  surface.  Single-color  radiometers  typically  respond  to  the 
average  temperature  in  the  field  of  view,  while  two-color  radiometers  tend  to  reflect  the  peak 
temperature  in  the  field  of  view.  Facility  geometry  usually  dictates  that  specimen  surfaces  be 
viewed  at  non-normal  incidence,  enlarging  the  surface  area  seen  by  the  radiometer.  The 
possibility  off  non-uniform  surface  temperatures  in  the  field  of  view  must  be  considered  for 
UHTC  specimens  with  small  stagnation  regions  like  sharp  wedges  or  cones. 

Radiometry  can  also  be  used  to  obtain  information  about  the  emittance  of  UHTC 
specimens  at  high  temperature.  Laboratory  measurements  of  UHTC  emittance  at  high 
temperature  are  difficult  and  rare.  The  emittance  measurements  of  Scatteia  et  al.  on  various 
UHTC  composites  (ZrB2-15SiC,  ZrB2-15SiC-2MoSi2,  and  ZrB2-15SiC-10HfB2)  are  notable  in 
this  regard.43,44  Scatteia  et  al.  have  also  demonstrated  how  UHTC  emittance  is  influenced  by 
surface  finish  and  oxidation.  Because  it  may  be  difficult  to  reproduce  oxidized  UHTC  surfaces 
in  the  laboratory  that  are  the  same  as  those  formed  in  the  low-pressure,  highly-dissociated 
oxygen  environments  of  arc-jet  and  ICP  flows,  in  situ  measurements  would  be  preferred. 

If  both  one-  and  two-color  radiometers  are  trained  on  the  same  surface  location,  surface 
temperature  can  be  measured  using  the  two-color  instrument  and,  with  the  surface  temperature 
known,  the  in  situ  emittance  can  be  derived  from  the  one-color  instrument.  Such  a  procedure 
was  used  by  Monteverde  and  Savino14  to  find  s  =  0.9  for  a  hot-pressed  ZrB2-15SiC  composite 
during  testing  in  a  80  kW  plasma  torch.  This  type  of  measurement  could  in  principle  be  done 
using  a  single  instrument,  by  calibrating  and  collecting  both  channels  of  a  two-color  radiometer. 
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A  further  extension  of  this  procedure  would  be  the  use  of  a  spectral  radiometer  to  capture 
wavelength-resolved  emission  intensity  over  a  broad  wavelength  range.  The  captured  intensity 
curve  is  fit  by  the  Planck  blackbody  radiation  function,  Eb  (A,  T ) ,  convoluted  with  a  wavelength- 

dependent  emittance,  yielding  both  surface  temperature  and  emittance  values.  Spectral 
radiometers  are  expensive  instruments  and  their  application  in  arc-jet  or  ICP  testing  of  UHTC 
materials  has  not  yet  been  reported. 

Emittance  is  a  function  of  wavelength,  emission  direction,  and  temperature.  The  surface 
energy  balance  requires  the  total  hemispherical  emittance  -  the  emittance  averaged  over  all 
wavelengths  and  emission  directions  -  as  input,  where  it  plays  a  crucial  role  in  the  determination 
of  surface  catalytic  efficiency  and  thus  the  computed  gas-phase  species  concentrations  above  the 
surface.  The  emittance  obtained  by  Monteverde  and  Savino14  yields  a  directional  value  valid 
over  the  wavelength  range  of  their  one-color  radiometer.  The  direct  use  of  this  value  in  the 
surface  energy  balance  implies  assumptions  of  wavelength  independence  (a  gray  surface)  and 
directional  independence  (a  diffuse  surface)  which  may  not  be  true.  Extension  of  this  method 
using  a  spectral  radiometer  would  produce  a  wavelength-dependant  directional  emittance  that 
can  be  averaged  over  wavelength  to  yield  a  total,  directional  emittance  value.  However,  some 
assumptions  or  approximations  would  still  have  to  be  made  about  the  relationship  between  total 
directional  and  total  hemispherical  emittance  before  this  derived  value  could  be  used  in  the 
surface  energy  balance.  Since  it  does  not  seem  possible  to  make  a  total  hemispherical  emittance 
measurement  in  situ,  the  best  available  alternative  is  probably  the  extrapolation  of  in  situ 
emittance  measurements  with  the  directional  and/or  spectral  dependencies  documented  in  pre- 
and  post-test  laboratory  measurements. 


4.2.  Gas  Emission  Spectroscopy 

The  radiation  emitted  by  gases  in  the  free-stream  or  in  the  near-surface  region  of  test 
specimens  contains  information  about  the  species  present.  Atoms  and  molecules  emit  radiation 
when  they  transition  from  higher  to  lower  energy  states.  Each  species  has  a  unique  set  of  energy 
levels  and  the  energies  of  emitted  photons  correspond  precisely  to  differences  between  these 
levels.  Emission  spectroscopy  collects  this  emitted  light  and  disperses  it  by  wavelength, 
generating  spectra  with  characteristic  intensity  features  (atomic  lines  and  molecular  bands)  that 
can  be  associated  with  individual  species.  Emission  is  a  path-integrated  measurement  technique 
without  spatial  resolution  along  the  light  collection  axis.  Reconstruction  of  the  spatial 
distribution  of  emission  requires  measurements  along  multiple  axes  and/or  assumptions  of 
symmetry  together  with  a  mathematical  procedure  known  as  an  Abel  inversion.45  Further 
complications  ensue  if  the  gas  is  optically  thick  and  significant  re-absorption  of  emitted  radiation 
occurs  along  the  collection  axis. 

During  the  Plasmatron  testing  of  hot-pressed  ZrB2-30SiC  materials  described  in  Section  3 
above,  prominent  bluish-green  emissions  were  observed  around  the  test  articles.  These 
emissions  first  appeared  and  intensified  as  the  sample  temperature  increased,  then  reached  a 
plateau,  and  finally  decreased  in  intensity  at  longer  test  times.  Follow-on  tests  were  performed 
in  which  emission  spectra  were  collected  directly  adjacent  to  the  specimen  surface  as  a  function 


A3-13 


of  test  time.  Light  was  captured  through  an  aperture,  focused  onto  the  end  of  an  optical  fiber  by 
a  spherical  mirror,  and  transmitted  to  a  spectrometer  with  a  wavelength  range  of  200  to  1100  nm 
and  a  resolution  of  0.25  nm.  Spectra  were  collected  at  1  Hz  during  the  test. 

These  measurements  confirmed  that  the  visible  emission  originated  from  electronically- 
excited  BO  and  BO2  molecules,  and  in  addition  detected  ultraviolet  emission  associated  with 
atomic  boron.  An  example  of  these  spectra  is  shown  in  Fig.  5,  which  plots  background- 
subtracted  raw  signal  against  emission  wavelength  (background  subtraction  removes  N2+ 
emissions  originating  from  the  Plasmatron  ffee-stream).  The  band  head  positions  for  the  BO 
(a2 n  ->  X2Z+  )  and  BO2  (a2 Hu  ->  X2Tlg  J  systems  are  from  Spalding  et  al.46  and  the  B 

(^1/2  ~^2  Ff  2,3/2}  doublet  was  identified  from  the  NIST  Atomic  Spectra  Database.47  Since  all 

boron  species  must  originate  from  ZrB2,  emission  spectroscopy  in  this  case  provides  an  in  situ 
monitor  for  the  oxidation  of  the  UHTC  composite  and  the  volatilization  of  boron  from  the 
resulting  oxide  scale. 


Fig.  5.  Raw  signal  versus  emission  wavelength  55  seconds  after  injection  of  sample  into  the 

freestream.  Background  emissions  from  the  Plasmatron  freestream  have  been  subtracted 
to  remove  contributions  from  the  N2+  band  system.  The  strong  line  at  589  nm  can  be 
assigned  to  sodium.  The  bandhead  positions  for  BO  and  BO2  are  from  Spalding  et  al.46 
and  atomic  lines  were  identified  from  the  NIST  Atomic  Spectra  Database.47  Arrows 
indicate  the  wavelengths  used  to  monitor  B,  BO,  and  BO2  as  a  function  of  test  time. 

UHTC  composites  containing  both  ZrB2  and  SiC  oxidize  to  form  complex  oxide  scales 
with  an  outer  glassy  layer  rich  in  silica  and  an  inner  layer  depleted  of  SiC.  The  structure  of  these 
oxide  scales  has  been  extensively  investigated  with  post-test  analytic  techniques  (microscopy,  X- 
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ray  diffraction,  chemical  analysis)26’48’49,  and  their  formation  hypothesized  in  terms  of  transport 
processes  and  thermodynamic  arguments.27’28 

In  dissociated-oxygen  environments,  the  passive  oxidation  of  zirconium  diboride  forms 
zirconia  and  boron  oxide,  and  the  passive  oxidation  of  silicon  carbide  forms  silica  and  carbon 
monoxide: 


ZrB2  +  50(g)  -»•  Zr02  +  B203  (3) 

SiC  +  30(g)  ^Si02  +  C0(g)  (4) 


The  oxidation  rate  of  ZrB2  at  moderate  temperatures  (below  1000°C)  is  significantly  faster  than 
the  oxidation  rate  of  SiC.  Therefore  the  initial  oxide  scale  is  expected  to  be  predominantly  B203. 
Amorphous  B203  has  a  very  low  softening  temperature  (-560  to  630°C)50  and  liquid  B203  seals 
the  surface  slowing  inward  oxygen  transport.  At  higher  temperatures,  the  SiC  oxidation  rate 
increases  and  the  oxide  scale  becomes  a  borosilicate  glass.  However  liquid  B203  has  a  much 
larger  vapor  pressure  than  silica  and  boron  oxides  are  predicted  to  volatilize  preferentially  from 
the  oxide  surface  leaving  a  silica  rich  glass.  As  the  glassy  scale  thickens  and  becomes  more 
silica  rich,  inward  oxygen  diffusion  slows  further  and  oxygen  concentrations  at  the  reaction 
interface  decrease.  Under  these  conditions  active  SiC  oxidation  becomes  favored  over  ZrB2 
oxidation,  slowing  B203  production  and  leading  to  the  formation  of  a  porous  SiC-depleted  sub¬ 
layer  via: 


SiC  +  20(g)  — ►  SiO(g)  +  CO(g) 


(5) 


Gaseous  SiO  formed  at  the  bottom  of  the  SiC-depleted  layer  diffuses  to  the  top  of  that  layer 
where  it  is  thought  to  condense,  augmenting  the  glassy  top  scale  from  below. 


In  detail,  this  process  is  quite  complicated,  involving  heterogeneous  media,  transitions 
between  different  controlling  chemistries,  the  growth  of  multiple  oxide  layers,  moving  interfaces, 
changing  glass  compositions,  evolving  transport  properties,  etc.  However,  a  simple  “shape 
function”  for  the  expected  time  variation  of  boron  species  volatilization  can  be  constructed  based 
on  the  premise  that  boron  oxide  is  formed  by  reaction  (3)  at  the  bottom  of  the  glassy  scale  and 
diffuses  through  the  scale  to  the  surface  where  it  evaporates  according  to  a  Hertz-Langmuir 
relationship.  The  diffusion  and  evaporation  fluxes  can  be  written  in  terms  of  the  boron  oxide 
concentrations  at  the  reaction  interface,  CB  O  f ,  and  the  outer  surface  ,  Crj  0^s ,  as 
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These  expressions  can  be  combined  by  assuming  diffusion  and  evaporation  fluxes  are  in  quasi- 
steady- state  at  any  time: 
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B203  ,i 


B203  ,ev 


RT\p/M\ |  A 


glass 


B203,v 
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b2o3 


(8) 


Both  the  vapor  pressure  and  the  diffusion  coefficient  are  presumed  to  follow  Arrhenius 
dependencies  on  temperature:  PBlo3,v  x  exp (^EV/RT^  and  DB  tG  oc  exp(-  Ed/RT).  If  parabolic 

scale  growth  is  assumed,  Aglass  <x  yft ,  and  if  falling  boron  oxide  production  with  growing  scale 
thickness  is  approximated  by  an  inverse  power  law,  CB^0i  i  cct~n ,  the  following  shape  function 
for  the  temperature  and  time  dependencies  of  the  boron  oxide  vaporization  flux  results: 


J 


B203  ,ev 


_ k£^ _ 

■Jr  expi/y, / RT )+  k2  4~t  exp(Z?D / RT )  ’ 


(9) 


where  kx ,  k2 ,  n  ,  Ev,  and  ED  are  constants. 

Figure  6  compares  the  experimental  temperature  and  emission  signal  history  (in  the 
left-hand  panel)  collected  during  a  Plasmatron  test  to  some  representative  volatilization  profiles 
(in  the  right-hand  panel)  computed  by  Eq.  (9)  for  a  piecewise  linear  temperature  history  and 
several  numerical  choices  of  constants.  The  experimental  emission  signals  were  collected  at  the 
wavelengths  indicated  in  Fig.  5;  the  experimental  temperature  was  collected  by  a  radiometer  with 
a  low  temperature  cut-on  of  ~1 100°C.  Ev  was  fixed  at  a  low  10  kJ/mole  on  the  premise  that 

boron  diffusion  in  high  temperature  amorphous  or  liquid  silica  is  very  fast.  ED  was  adjusted  to 
320  kJ/mole  to  reproduce  the  sharp  rise  in  boron  emission  signatures  observed  between  40  and 
50  seconds.  The  activation  energy  predicted  by  thermodynamic  calculations  for  the  vapor 
pressure  of  B2O3  over  liquid  boron  oxide  is  closer  to  380  kJ/mole25’51.  The  decay  of  the  shape 
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function  with  time  is  controlled  by  the  exponent  n.  Computations  are  shown  for  n  =  0,  0.5  and 
1.0,  with  k2  fixed  to  10 14  and  kx  adjusted  to  maintain  similar  peak  magnitudes  for  all  profiles 

and  to  put  the  shape  function  on  the  same  (arbitrary)  scale  as  the  emission  data. 
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Fig.  6.  Experimental  temperature  and  emission  signal  history  collected  during  a  Plasmatron  test 
of  a  hot-pressed  ZrB2-30SiC  specimen  (left-hand  panel)  and  representative  volatilization 
profiles  computed  by  Eq.  (9)  (right-hand  panel).  Emission  was  collected  at  the  three 
wavelengths  indicated  in  Fig.  5;  the  radiometer  cut-on  temperature  was  ~1 100°C.  The 
B2O3  fluxes  and  concentrations  leading  to  the  shape  function  model  are  illustrated  in  the 
right-hand  panel.  Shape  functions  were  computed  with  k2  =  1014 ,  Ev=  320  kJ/mole,  and 
Ed  - 10  kJ/mole,  and  the  values  of  n  and  /c,  listed  on  the  figure. 


In  general  this  very  simple  model  captures  the  qualitative  trends  seen  in  the  emission  data 
quite  well.  The  largest  exponent  (n  =  1)  seems  to  reproduce  the  relatively  steep  emission  signal 
decay  the  best.  This  is  consistent  with  the  idea  that  the  concentration  of  B2O3  produced  at  the 
reaction  interface  should  vary  inversely  with  time  at  least  as  fast  as  the  scale  grows  and  probably 
faster  as  the  decreasing  oxygen  pressure  at  the  interface  begins  to  favor  active  SiC  oxidation  over 
ZrB2  oxidation. 
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At  higher  temperatures  (above  ~1750°C),  the  vapor  pressure  of  silica  increases 
significantly  and  a  similar  process  of  volatilization  from  the  UHTC  oxide  scale  should  lead  to  Si- 
containing  species  in  the  gas  phase.  Hirsch  et  al.52,  Altman  et  al.53,  and  Jentschke  et  al.54  have 
monitored  Si  atom  densities  in  front  of  C/C-SiC  composites  exposed  to  nitrogen-oxygen  plasma 
flows,  using  high-resolution  spectroscopy  of  Si  I  multiplet  emissions  in  the  250-253  nm  range. 
Herdrich  et  al.55  observed  Si  emission  near  252  nm  and  288nm  while  testing  SiC  specimens  in 
oxygen-nitrogen  plasmas,  and  also  captured  emission  from  SiCL  molecules  near  423  nm. 
Emission  from  SiO  molecules  was  not  reported  in  any  of  these  studies,  which  is  surprising  given 
that  reaction  mechanism  (5)  should  have  been  operating  under  some  of  the  test  conditions 
reported. 


4.3.  Laser  Induced  Fluorescence 

Laser  induced  fluorescence  (LIF)  is  a  species  selective,  non-intrusive  diagnostic  that  has 
been  widely  applied  for  the  characterization  of  combustion  and  plasma  environments.56,57  A 
pulsed,  tunable  laser  source  is  used  to  generate  monochromatic  light  at  a  unique  atomic  or 
molecular  absorption  wavelength.  The  absorbed  photon  energy  excites  the  target  species  to  a 
higher  electronic  energy  level,  from  which  it  subsequently  decays  to  a  lower  energy  level  (or 
levels)  emitting  radiation  at  a  characteristic  wavelength  (or  wavelengths). 

One  advantage  of  the  LIF  technique  is  that  ground-state  species  populations  can  be 
probed,  whereas  emission  spectroscopy  detects  only  electronically-excited  species.  Even  at 
relatively  high  temperatures,  the  ground  electronic  state  populations  typically  predominate  over 
the  populations  in  higher  electronic  levels.  Another  advantage  of  LIF  is  that  spatially  resolved 
measurements  are  more  easily  made  and  interpreted  than  by  emission  spectroscopy,  since  the 
fluorescence  collection  optics  field-of-view  and  the  excitation  laser  beam  axis  can  be 
independently  oriented  to  intersect  at  a  desired  location,  as  shown  schematically  on  the  left  side 
of  Fig.  7.  Generally  speaking  LIF  detection  works  best  for  light  atoms  and  diatomic  molecules. 

Applications  of  LIF  to  plasma  stream  characterization  typically  involve  measurement  of 
three  variables:  translational  temperature,  convective  velocity,  and  species  concentration.  Air 
and  nitrogen  plasma  species  probed  by  LIF  have  included  O  atoms,  N  atoms,  and  nitric  oxide 
(NO).58"62  For  NO,  single-photon  excitation  in  either  the  ultraviolet  y  band  (A  <—  x)  or  [3  band 
(B  <—  X)  is  possible.  However,  the  large  energy  level  spacing  of  the  oxygen  and  nitrogen  atoms 
requires  a  two-photon  excitation  scheme,  and  the  excitation  cross  sections  are  correspondingly 
smaller.63"66  The  energy  levels  involved  in  common  two-photon  excitation  schemes  for  atomic 
oxygen  and  nitrogen  are  shown  in  the  center  and  right  side  of  Fig.  7.  Atomic  nitrogen  is  a  useful 
target  species  for  highly  dissociated  air  or  nitrogen  plasmas,  because  its  concentration  is  a 
sensitive  indicator  of  the  degree  of  plasma  dissociation.58,59  Atomic  oxygen  number  densities  are 
critical  for  quantifying  both  passive  and  active  surface  oxidation  processes.  Both  atomic  species 
are  involved  in  exothermic  surface  recombination  processes  that  can  contribute  substantially 
aerothermal  heating. 
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Fig.  7.  Energy  level  diagrams  for  two  common  two-photon  LIF  schemes  for  probing  atomic 
oxygen  and  atomic  nitrogen.  The  spatial  resolution  of  the  LIF  approach  is  indicated  by 
the  cartoon  at  the  left. 


The  experimental  arrangement  for  probing  either  atomic  oxygen  or  nitrogen  is  essentially 
the  same,  as  shown  in  Fig.  8.  A  Nd:YAG-pumped  dye  laser  provides  the  fundamental  output, 
which  is  converted  by  the  appropriate  tripling  crystal  arrangement  to  the  desired  ultraviolet,  two- 
photon  wavelength.  As  indicated  in  the  figure,  the  ultraviolet  output  is  split  into  three  different 
paths.  The  first  is  directed  toward  the  plasma  stream  (either  ICP  or  arc-jet),  while  the  other  two 
paths  are  directed  toward  a  NO  reference  cell  and  a  microwave-discharge  flow  reactor, 
respectively.  The  NO  cell  is  used  to  assess  the  laser  performance  and  to  guide  laser  tuning 
toward  the  relatively  sparse  two-photon  excitation  wavelengths.  Within  the  flow  reactor,  the 
microwave  discharge  creates  a  stable,  non-equilibrium  population  of  atoms  at  known  pressure 
and  temperature,  and  the  fluorescence  from  this  population  can  be  used  to  calibrate  the 
fluorescence  signals  from  the  plasma  stream.  This  calibration  enables  measurement  of 
translation  temperature,  flow  velocity  (depending  on  the  facility  and  beam  orientation),  and 
number  density  of  the  target  species. 
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Fig.  8.  Schematic  of  the  experimental  configuration  for  two-photon  LIF  implementation  in  arc- 
jet  and  ICP  facilities.  In  addition  to  measuring  LIF  from  the  plasma  facility, 
measurements  are  made  in  a  nitric  oxide  cell  for  laser  wavelength  monitoring  and  in  a 
flow  reactor  for  calibration  of  the  facility  measurements.  The  temperature  and  atomic 
number  density  in  the  flow  reactor  are  known  independently  of  the  laser  measurements. 


An  example  of  atomic  nitrogen  fluorescence  acquired  in  an  ICP  facility  is  shown  in  Fig. 

9,  where  the  atomic  nitrogen  fluorescence  is  plotted  as  a  function  of  the  dye  laser  fundamental 
wavelength.  The  broader  of  the  two  traces  is  the  signal  collected  from  plasma  stream  at  the 
boundary  layer  edge  and  the  narrower,  less  noisy  signal  is  from  a  microwave-discharge  flow 
reactor.  The  smooth  curves  for  each  LIF  signal  are  nonlinear  least-squares  fits  to  the  line 
profiles  using  a  spectral  model  that  includes  the  relevant  broadening  mechanisms  and  uses  the 
total  line  width  as  a  fitting  parameter.  Line  shape  fit  results  from  the  flow  reactor  signal  are  used 
to  extract  information  from  the  plasma  stream  measurements  as  explained  in  the  following 
paragraphs.  Finally,  the  temperature  value  derived  for  this  particular  measurement  is  indicated 
on  the  figure,  along  with  the  flow  reactor  temperature  used  to  extract  the  laser  line  width. 

The  translation  gas  temperature  is  derived  from  a  measurement  of  the  total  width  of  the 
transition  once  the  laser  line  width  is  determined  from  the  flow  reactor  line  shape  as, 
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where  M  is  the  molar  mass,  c  is  the  speed  of  light,  ks  is  the  Boltzmann  constant,  nA  is 
Avogadro’s  number,  and  v  is  the  transition  frequency  in  cm'1.  The  different  widths  A  v2  (also 

cm'1)  are  labeled  with  subscripts  that  designate  Total  or  Doppler  and  Flow  or  Reactor.  The  flow 
velocity  can  be  determined  from  the  Doppler  shift  of  the  central  transition  wavelength,  A,  as 


v  = 


c  A  A 
A  cos# 


(11) 


and  only  a  separate  measurement  of  the  angle  between  the  velocity  and  laser  propagation 
vectors,  ,  is  needed.  While  the  Doppler  shift  in  the  free  stream  of  supersonic  arc-jet  facilities  is 
quite  large,  and  can  be  used  to  reliably  measure  the  free  stream  velocity,  the  Doppler  shift  in 
subsonic  ICP  facilities  is  too  small  for  velocity  determination,  as  is  clearly  evident  in  Fig.  9. 


WAVELENGTH,  nm 

Fig.  9.  Two-photon  LIF  signals  from  excitation  of  the  207  nm  transition  of  atomic  nitrogen  at 
the  boundary  layer  edge  in  the  VKI  Plasmatron  ICP  Torch  Facility.  Despite  the  fact  that 
the  beam  is  angled  with  respect  to  the  axial  stream  velocity  (as  shown  in  Fig.  8)  the 
subsonic  flow  speed  does  not  provide  a  measurable  shift  between  the  flow-reactor  and 
ICP  stream  transition  line  centers.  The  spectral  fit  assumes  that  Doppler  broadening  is 
the  dominant  line  broadening  mechanism. 
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Finally,  the  local  number  density  of  atomic  nitrogen  is  determined  from  integrated  LIF 
signal.  Several  additional  measurements  are  required  to  obtain  absolute  values  of  number 
density,  including  characterization  of  the  fluorescence  collection  efficiency,  the  temporal  and 
spatial  characteristics  of  the  laser  beam,  and  finally  the  use  of  either  a  known  N-atom  reference 
population,  a  rare  gas  calibration,67'69  or  a  two-photon  excitation  cross-section  measurement.64'66 
For  the  latter  approach,  the  expression  for  nitrogen  atom  number  density  is 
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where  Sn  is  the  LIF  signal  in  volts,  r,  are  the  radiative  and  observed  lifetimes,  AP  is  the  beam 
area,  D  is  the  collection  optics  calibration  factor,  F2(t)  is  the  square  of  the  temporal  laser  pulse 
shape  and  the  G(2)<d2)  product  represents  the  two-photon  LIF  cross  section. 

Two-photon  LIF  has  been  implemented  in  arc-jets  to  characterize  the  free  stream 
conditions58,59  and  to  investigate  free  stream  property  gradients.61,70  By  measuring  the  three 
quantities  noted  above  in  nitrogen/argon  plasma  flows,  and  making  use  of  a  pitot  pressure 
measurement  and  facility  data,  it  is  possible  to  calculate  the  stream  total  enthalpy  and  to  quantify 
the  different  contributions:  thermal,  kinetic,  and  chemical.  This  is  extremely  useful  for 
establishing  stream  conditions,  extrapolating  the  free  stream  test  conditions  to  a  flight 
environment,  and  provides  experimental  data  against  which  a  CFD  computation  of  gas  flow 
through  the  arc-jet  nozzle  and  towards  the  test  article  can  be  tested.  However,  arc -jet  free  stream 
measurements  still  only  provide  a  different  starting  point  for  the  eventual  assessment  of  the 
boundary  layer  conditions  above  a  test  article. 

The  arc-jet  free  stream  measurements  indicate  low  translational  temperatures,  high 
velocity  and  non-equilibrium  atom  number  densities.58,59  In  contrast,  the  boundary  layer  edge 
measurement  from  the  ICP  facility  discussed  above  (Fig.  9)  indicates  high  translational 
temperature  and  low  velocity.  It  is  reasonable  to  expect  that  within  the  boundary  layer  in  an  arc- 
jet  test,  LIF  measurements  would  be  quite  different  from  those  of  the  free  stream,  and  more 
likely,  similar  to  those  of  the  ICP  facility  boundary  layer  edge. 

For  material  test  applications  temperature  and  species  information  from  the  reacting 
boundary  layer  is  extremely  important.  While  LIF  measurements  in  the  boundary  layer  are 
relatively  rare,  recent  experiments  indicate  their  feasibility.60,70  Two  strategies  exist  for 
measuring  species  and  temperature  profiles  above  the  surface.  The  first  is  a  point-wise  approach 
that  involves  translating  the  laser  probe  volume  toward  the  surface,  as  indicated  in  Fig.  10,  and 
the  other  is  a  planar  approach  that  requires  sufficient  laser  pulse  energy  to  spread  the  beam  into  a 
light  sheet.  In  the  latter,  an  intensified  camera  is  needed  to  record  the  fluorescence  signals, 
which  complicates  the  interpretation  of  the  measured  signals. 
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Fig.  10.  Translation  of  either  the  test  article  or  the  beam  delivery  and  collection  optics  to 

measure  the  species  number  densities  and  temperature  in  the  boundary  layer  approaching 
the  test  article  in  the  stagnation  region. 


Figure  1 1  shows  computed  boundary  layer  profiles  of  temperature,  and  total  and  species 
number  densities,  for  an  air  plasma  along  the  stagnation  streamline  for  a  particular  high-enthalpy 
test  condition.17  The  catalytic  efficiency  of  the  surface  was  set  to  a  relatively  low  10'4  for  these 
computations.  Temperature  and  number  densities  are  normalized  by  their  boundary  layer  edge 
values.  Different  information  is  available  from  different  species  measurements.  All  species 
number  densities  increase  with  approach  to  the  surface  because  of  the  decreasing  gas 
temperature  in  a  constant  pressure  boundary  layer.  The  relative  O-atom  number  density  and  total 
mixture  number  density  curves  overlay  each  other  in  Fig.  11,  indicating  that  for  these  plasma 
conditions,  relatively  little  gas-phase  recombination  of  O  to  O2  occurs.  The  N-atom  number 
density  increases  more  slowly  and  the  N2  number  density  more  rapidly  than  the  total  number 
density,  indicating  that  gas-phase  recombination  of  N  to  N2  is  important  at  these  conditions. 

Whether  probing  N  or  O,  if  a  simultaneous  measurement  of  the  atom  number  density  and 
translational  temperature  can  be  obtained  (as  described  above)  as  a  function  of  distance  from  the 
sample  surface,  the  total  number  density  can  then  be  computed  from  a  pressure  measurement  and 
the  translational  temperature  using  the  ideal  gas  equation  of  state  n  =  P/kBT  .  This,  along  with 

the  measured  atom  number  density,  gives  the  mole  fraction  of  the  probed  atomic  species.  For  a 
pure  nitrogen  plasma  this  can  completely  characterize  the  local  thermodynamic  conditions,  as 
long  as  ionization  is  not  significant,  since  the  N2  number  density  can  be  found  from  nN  =  n  —  nN 

.  For  an  air  plasma,  measurements  of  both  O-atom  and  N-atom  number  densities  are  needed  for 
a  complete  understanding  of  the  local  thermodynamic  state. 
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The  relative  species  density  profiles  presented  in  Fig.  11  show  trends  characteristic  of  a 
low  catalycity,  non-reacting  surface,  wherein  the  relative  number  densities  increase  towards  the 
wall  as  determined  by  the  gas  temperature  profile  and  gas-phase  recombination  reactions.  For  a 
highly  catalytic  surface  the  atomic  species  number  densities  should  trend  downward  approaching 
the  surface  to  very  small  levels.  Reported  species  density  uncertainties60,70  suggest  that 
measurement  of  gradients  with  density  changes  of  two  to  four  are  realizable,  especially  when 
relative  LIF  signals  are  used  to  measure  the  trend  of  a  normalized  species  density. 


Fig.  11.  Computed  temperature,  mixture  density,  and  species  number  density  variations  along 
the  stagnation  streamline  for  the  flow  over  a  low-catalycity  sample  in  the  VKI 
Plasmatron  ICP  Torch  Facility.  The  number  densities  are  normalized  by  their  free  stream 
(boundary  layer  edge)  values.  Note  that  the  total  and  O-atom  number  density  curves 
overlay  one  another,  indicating  that  for  these  boundary  layer  conditions  there  is  little  gas- 
phase  recombination  of  O  to  O2. 

The  monitoring  of  volatile  products  during  high-enthalpy  UHTC  testing  by  LIF  is  also 

71  79  79 

possible,  although  such  work  has  not  yet  been  reported.  Both  B  atoms  and  BO  molecules  ’ 

7  A  nr  nA  nr 

have  been  detected  by  LIF  schemes,  as  have  Si  atoms  ’  and  SiO  molecules  ’  .  Feigl  and 
Auweter-Kurtz77  have  used  LIF  to  monitor  SiO  concentrations  in  front  of  SiC  materials  oxidized 
in  argon-oxygen-nitrogen  plasma  flows,  by  exciting  the  SiO  (axD  <—  X  'l+  (0,0))  transition  near 
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234  nm  and  collecting  fluorescence  near  260  nm  from  the  (a  lD  — »  X  ‘S+  (0,3))  band.  They 
associate  rapid  rises  of  SiO  in  the  gas  phase  with  the  transition  between  passive  and  active  SiC 
oxidation  (i.e.,  between  reactions  (4)  and  (5)).  Beyond  the  simple  detection  of  volatile  species, 
spatially-resolved  LIF  of  volatiles  in  the  boundary  layer  would  enable  the  measurement  of 
species  gradients  that  could  in  principle  be  related  to  volatilization  rates  from  the  surface.  No 
measurements  of  this  type  have  yet  been  attempted  and  their  interpretation  would  require 
additional  gas-phase  chemistry  inputs  to  existing  CFD  programs. 


4.4.  Other  Techniques 

4.4.1  Fourier  Transform  Infrared  Spectroscopy 

Fourier  Transform  Infrared  (FTIR)  spectroscopy  is  a  technique  used  to  capture  and 
spectrally  resolve  infrared  radiation  using  a  scanning  interferometer,  typically  at  wavelengths 
between  1  and  40  pm.  Transmission  FTIR  spectroscopy  is  often  used  as  a  laboratory  technique 
to  identify  chemical  species  in  gas  or  liquid  samples,  by  absorption  features  that  occur  at 
characteristic  energies  (wavelengths)  related  to  particular  vibrational  modes  of  different  chemical 
bonds.  These  characteristic  absorption  features  are  relatively  sharp  in  the  gas  phase,  and  often 
broaden  and  shift  slightly  in  the  liquid  phase.  A  similar  technique  is  used  to  detect  molecules  on 
surfaces  by  reflection  FTIR  and  can  be  used  to  monitor  thin  fdm  growth. 

Estimates  of  the  temperature-dependant  total,  hemispherical  emittance  are  often  derived 
from  spectrally-resolved,  hemispherical  reflectance  measurements  performed  at  room 
temperature  using  FTIR  instruments  fit  with  integrating  spheres.17  The  measured  spectral 
reflectance,  p29iQ l),  is  converted  to  spectral  absorbance,  a29S(X),  using  the  relation 

cc 298  (A)=  1  -  p298  (A)  valid  for  an  opaque  surface,  and  the  spectral  absorbance  is  equated  to  the 
spectral  emittance  based  on  Kirchhoff  s  law,  £29s{^)=  a29$(^)-  Total  hemispherical  emittance  at 

any  temperature  is  then  computed  by  averaging  the  room-temperature  emittance  over  the  Planck 
blackbody  radiation  function: 

^2 

{^298  (TX,(T,T)/A 

e(T)=^ -  .  (13) 

\Eb(A,T)dA 

A 


This  process  assumes  that  the  temperature  dependence  of  total  emittance  is  dominated  by  the 
temperature  dependence  of  the  Planck  function  and  not  by  the  optical  constants  of  the  material. 

It  also  presumes  that  the  spectral  range  of  the  measurement,  f  to  Z2 ,  is  sufficiently  broad  that  it 

contains  the  spectral  region  dominating  the  emissive  power  at  each  temperature. 
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A  variant  of  this  technique  -  emission  FTIR  -  was  recently  used  in  a  laboratory  setting  to 
capture  thermal  emission  from  C/ZrB2-SiC  composites  over  the  200  to  800°C  temperature  range 

no 

and  an  effective  temperature-dependant  emittance  was  computed.  However,  the  spectral  range 
of  the  FTIR  instrument  used  in  this  work  was  only  8-14  pm  and  this  limited  range  contains  just 
-40%  of  the  radiative  power  emitted  at  200°C  and  less  than  10%  of  the  emitted  power  at  800°C. 
The  derived  emittance  values  are  thus  tenuous  approximations  for  the  total  emittance. 

1  7 

Marschall  et  al.  have  shown  how  a  broad  feature  in  the  FTIR  reflectance  spectra  of 
virgin  ZrB2-30SiC  in  the  10-13  pm  range  was  removed  and  replace  by  a  sharper  feature  centered 
near  9  pm  after  oxidation.  The  former  feature  can  be  related  to  stretching  mode  vibrations  of 
Si-C  bonds  and  the  latter  to  stretching  mode  vibrations  of  Si-0  bonds.  A  similar  enhanced 
reflectance  feature  predicted  for  Zr02  near  14  pm  was  not  observed  in  oxidized  samples, 
consistent  with  the  fact  that  the  outer  scale  was  dominated  by  silica-rich  glass. 

Related  spectral  features  should  appear  as  emittance  modifications  to  the  Planck 
blackbody  radiation  function  in  thermal  emission.  This  suggests  that  FTIR  emission 
spectroscopy  might  be  used  as  an  in  situ  monitor  for  changing  surface  compositions  during  high- 
enthalpy  testing.  Measurements  of  this  type  were  made  by  Hirsch  et  al.52  for  characterizing  the 
oxidation  behavior  of  carbon/carbon  and  carbon/silicon  carbide  composites  in  the  ICP  wind 
tunnels.  However,  this  concept  has  never  been  explored  as  an  in  situ  diagnostic  for  UHTC 
materials  systems. 

Even  without  reference  to  a  blackbody  standard,  the  broad  mid-IR  spectral  features 
associated  with  SiC>2,  SiC,  ZrC>2  and  Hf02  are  sufficiently  separated  that  difference  spectra 
should  track  shifts  in  the  dominant  surface  composition  with  time.  This  approach  has  the 
potential  to  monitor  important  thermodynamic  and  chemical  transitions  as  test  conditions  are 
varied  or  as  constituents  are  depleted  over  time.  For  example,  changes  in  the  relative  strengths 
of  the  SiC>2  and  SiC  spectral  features  could  be  used  to  track  the  important  transition  from  passive 
SiC  oxidation  with  the  formation  of  condensed  silica  to  active  SiC  oxidation  with  the  formation 
of  volatile  SiO.  The  identification  of  temperature  and  pressure  thresholds  at  which  the  loss  of 
SiC>2  from  the  oxide  scale  becomes  large — be  it  by  active  oxidation,  evaporation,  or  melt  flow — 
is  critical  to  UHTC  performance  and  should  be  evidenced  by  the  appearance  of  the  ZrC>2  or  HfC>2 
features  above  14  pm. 


4.4.2  Absorption  Spectroscopy 

The  spatial  resolution  advantage  of  the  LIF  technique  has  been  contrasted  above  with 
path-integrated  emission  measurements,  which  can  also  clearly  provide  useful  information  for 
material  test  applications.  Similarly,  absorption  measurements  can  provide  useful  information 
since  they  typically  probe  either  ground  or  low-lying  electronic  states  even  though  this 
information  is  also  path-integrated.  One  such  application  is  the  measurement  of  molecular 
species  formed  by  gas/surface  interactions. 

Advances  in  semiconductor  processing  have  enabled  rapid  development  of  Tunable 
Diode  Laser  Absorption  Spectroscopy  (TDLAS).  Relatively  inexpensive,  but  tunable,  diode 
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lasers  operating  with  narrow  line  widths  can  be  optimized  for  overlap  with  favorable  rotational- 
vibrational  transition  features  of  important  molecular  species  that  have  dipole  moments,  such  as 
NO,  CO2,  CO,  H20,  SiO,  and  BO.  Access  to  ground  electronic  state  population  information  is 
often  available  by  probing  rotational-vibrational  transitions  in  the  1  to  5  pm  wavelength  region. 
Neither  02  nor  N2  have  dipole  moments,  and  so  are  not  considered  candidate  species  for  this 
approach;  however,  in  the  case  of  02,  it  is  possible  to  probe  the  ground  state  by  absorption  via 
the  A  electronic  band  transition  ( b'  2+  <—  X3  Z~ )  near  760  nm.  This  approach  has  been  used  to 

characterize  shock  tube  flows,79  but  has  not  yet  been  applied  to  either  arc -jet  or  ICP  facility 
measurements. 

Recently,  the  TDLAS  technique  was  used  to  probe  CO  in  a  Mars  atmosphere  study  in  an 
arc-jet  facility,  where  path-integrated  values  of  translation  temperature,  velocity,  and  CO 
concentration  were  obtained.80  To  date  however,  the  application  of  TDLAS  in  the  boundary 
layer  near  a  catalytic  or  reacting  material  surface  in  a  plasma  facility  to  quantify  recombined  or 
reaction-produced  molecules  has  not  been  reported. 


4.4.3  Raman  Spectroscopy 

Raman  scattering,  like  LIF,  is  an  inelastic  scattering  process,  but  involves  excitation  to  a 
virtual  upper  energy  level.  In  Raman  spectroscopy,  a  laser  beam  illuminates  the  gas  mixture  and 
Raman  scattered  photons  are  detected  at  wavelengths  offset  from  the  illumination  wavelength  by 
increments  determined  by  the  vibrational  and  rotational  energy  spacing  of  molecules  in  gas.  By 
choosing  an  appropriate  laser  wavelength  (typically  in  the  ultraviolet  since  Raman  cross  sections 
scale  as  frequency  to  the  fourth  power),  and  using  a  high-resolution  spectral  filter  (often  a 
spectrometer),  scattered  signals  from  molecules  in  a  gas  mixture  can  be  recorded  and  analyzed 
for  species  concentration  and  temperature.  Details  of  this  common  spectroscopic  technique  can 
be  found  in  Eckbreth.56 

An  attractive  attribute  of  Raman  scattering  is  that,  as  in  the  case  of  LIF,  it  provides 
spatially  resolved  information,  and  so  can  be  used  to  map  temperature  and  species  gradients.  In 
addition,  each  molecular  species  present  scatters  at  its  characteristic  Raman  wavelengths,  so  a 
single  laser  frequency  can  access  multiple  species,  depending  on  the  strength  of  their  Raman 
cross  sections.  Since  all  molecules  are  Raman-active,  this  technique  can  detect  important  species 
involved  in  UHTC  oxidation  and  volatilization  processes,  like  SiO,  Si02,  BO,  B02,  CO,  and 
C02. 

While  applications  to  plasma  flows  in  arc -jet  and  ICP  facilities  are  relatively  sparse,  one 
recent  investigation  has  demonstrated  the  potential  of  this  technique  in  the  area  of  recombined 
molecular  species  quantification.81  In  this  investigation  a  frequency  tripled  Nd:YAG  laser 
operating  at  355  nm  was  used  to  probe  the  boundary  layer  of  an  air  plasma  flowing  over  a  flat 
stainless  steel  plate  in  an  ICP  facility.  Despite  the  relatively  weak  Raman  cross  sections 
(~10'31  cm2  sr'1)  profiles  of  N2  and  02  number  density,  rotational  temperature,  and  vibrational 
temperature  were  obtained  in  the  boundary  layer  flow.  Long  steady  test  times  were  needed  to 
achieve  usable  signal-to-noise  levels,  but  the  results  are  promising  for  future  applications. 
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Summary 

Arc-jet  and  ICP  plasma  tunnels  remain  the  primary  facilities  used  to  ground  test  the 
performance  of  UHTC  TPS  materials  and  components  for  hypersonic  flight.  We  have 
highlighted  the  main  features  of  the  test  conditions  obtained  in  these  two  types  of  facilities  and 
discussed  their  similarities  and  differences.  Various  in  situ  optical  diagnostics,  both  established 
and  under  development,  were  described,  that  can  provide  information  to  document  test 
conditions,  monitor  the  evolution  of  test  articles,  and  provide  data  that  can  be  used  to  strengthen 
the  accuracy  of  CFD  models.  The  importance  of  CFD  modeling  of  the  arc-jet  and  ICP  test 
environments,  and  in  particular  the  computation  of  species  concentrations  at  specimen  surfaces, 
is  emphasized  as  essential  for  materials  science  interpretations  of  UHTC  test  results.  The 
increasing  integration  of  CFD  modeling  and  in  situ  optical  diagnostics  in  plasma  tunnel  testing  is 
sure  to  enhance  the  understanding  and  optimization  of  UHTC  material  performance  in 
aero  thermal  environments. 
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Abstract 

Analytical  modeling  of  thermal  and  mechanical  response  is  a  fundamental  step  in  the 
design  process  for  ultra-high  temperature  ceramic  components,  such  as  nose  tips  and  wing 
leading  edges  for  hypersonic  applications.  The  purpose  of  the  analyses  is  to  understand  the 
response  of  test  articles  to  high-enthalpy  flows  in  ground  tests  and  to  predict  component 
performance  in  particular  flight  environments.  Performing  these  analyses  and  evaluating  the 
results  require  comprehensive  and  accurate  physical,  thermal,  and  mechanical  properties.  In  this 
paper,  we  explain  the  nature  of  the  analyses,  highlight  the  essential  material  properties  that  are 
required  and  why  they  are  important,  and  describe  the  impact  of  property  accuracy  and 
uncertainty  on  the  design  process. 


1.  Introduction 

Certain  characteristics  of  ultra-high  temperature  ceramic  (UHTC)  materials  containing 
transition  metal  borides  and  carbides  give  them  a  potential  advantage  over  more  traditional 
materials  for  use  in  hypersonic  applications.  In  particular,  hafnium  and  zirconium  diborides  and 
carbides,  as  well  as  their  oxides,  hafnia  and  zirconia,  have  extremely  high  melting  points,  all  in 
excess  of  2500  °C  (4530°F).1  While  HfC  and  ZrC  have  higher  melting  points  than  HfB2  and 
ZrB2,  the  diborides  have  substantially  higher  thermal  conductivities  than  the  carbides.2,3  This 
combination  of  high-temperature  capability  and  high  thermal  conductivity  make  HfB2  and  ZrB2 
particularly  attractive  for  use  in  sharp  wing  leading  edges  (WLEs)  and  nose  tips.4  In  practice,  the 
high-temperature  oxidation  resistance  of  pure  diboride  materials  is  not  sufficient  for  aerothermal 
flight  environments.  The  best  oxidation  performance  is  found  for  monolithic  materials  hot- 
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pressed  from  mixtures  of  ceramic  powders  containing  a  silica  former  as  a  minor  component.5,6 
The  most  common  compositions  typically  contain  10  to  30  volume  percent  SiC,7'12  although 

1  O  1  /T 

other  silica  formers  like  MoSi2  and  TaSi2  have  also  been  investigated. 

The  design  of  high-performance  hypersonic  vehicles  generally  involves  relatively  sharp  nose 
tips  and  WLEs.  (An  arcjet  model  of  a  concept  UHTC  WLE  section  is  shown  in  Fig.  1.)  While 
not  an  achievable  design  in  practice,  in  the  context  of  supersonic  flow  theory,  a  “sharp”  leading 
edge  is  one  with  zero  radius  of  curvature.  In  a  practical  design  context,  “sharp”  usually  refers  to 
a  leading  edge  radius  that  is  much  smaller  than  the  nose,  wing,  or  vehicle  length  scale,  while 
“blunt”  implies  the  converse. 


Fig.  1.  UHTC  wing  leading  edge  arcjet  model. 


Sharp  leading  edges  help  reduce  the  vehicle’s  drag,  enhance  maneuverability  and 
performance,  and  also  improve  safety  due  to  an  increased  cross-range  capability.17  The  WLE  of 
the  Space  Shuttle  orbiter  has  a  radius  of  approximately  10  cm  at  the  tip,  while  design  concepts 
for  some  hypersonic  maneuvering  vehicles,  such  as  the  one  under  development  in  Defense 
Advanced  Research  Projects  Agency  (DARPA)’s  Falcon  program,  have  leading  edge  radii  on  the 
order  of  millimeters.  The  challenge  with  sharp  WLE  designs  is  that  the  convective  heating  to  the 
surface,  and  hence  the  surface  temperature,  increases  as  the  WLE  radius  decreases.  The  surface 
temperature  on  the  stagnation  region  of  such  sharp  leading  edges  can  potentially  exceed  2000  °C 
(3632  °F).  Few  materials  can  survive  extended  periods  at  these  temperatures  and  retain  their 
dimensional  and  structural  integrity  in  an  aggressive  oxidizing  flow  environment. 

The  performance  advantage  of  diboride-based  UHTC  materials  comes  not  only  from  their 
high-temperature  capability,  but  also  from  their  high  thermal  conductivity.  Convective  energy 
that  enters  the  surface  near  the  stagnation  region  is  conducted  away  to  cooler  regions  of  the 
leading  edge,  where  it  can  be  radiated  back  to  the  environment  (Fig.  2).  The  higher  the  thermal 
conductivity  of  the  leading  edge  material,  the  more  efficient  this  process  becomes.  The  UHTC 
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leading  edge  then  behaves  much  like  a  passive  heat  pipe,  to  move  energy  through,  and  eventually 
out  of,  the  system. 


Radiation  back  to 
the  environment 


High  heat  flux  to  the 
stagnation  region 


i  Conduction  away 
from  the  nose 


T7 


Attachment  region 
remains  relatively  cool 


Fig.  2.  UHTC  leading  edge  thermal  management  concept. 


The  hypersonic  applications  of  UHTC  materials  involve  complex  interactions  between  the 
material  and  the  aerothermodynamic  environment.  Although  most  designers  would  choose  to 
build  concept  vehicles  and  conduct  flight  tests  in  relevant  environments,  such  testing 
opportunities  are  rare,  often  impractical,  and  always  very  expensive.  Consequently,  the  design 
process  for  UHTC  components  for  hypersonic  vehicles  is  heavily  dependent  on  computational 
methods  and  ground-based  testing.  Analytical  performance  predictions  and  the  interpretation  of 
ground-based  testing  require  a  self-consistent  set  of  accurate  material  properties  with  well- 
defined  uncertainties.  Because  UHTC  components  have  to  operate  in  environments  from  room 
temperature  to  nearly  2000  °C  (3630°F),  it  is  important  that  most  UHTC  properties  be  measured 
over  this  entire  temperature  range. 

The  purpose  of  this  paper  is  to  describe  the  material  properties  that  are  most  important  for 
UHTC  component  design  and  the  reasons  why.  Section  2  provides  a  brief  introduction  to 
hypersonic  aerothermodynamics;  to  understand  why  a  UHTC  material  might  be  used  in  a 
hypersonic  application,  materials  developers  must  have  at  least  a  basic  understanding  of  the 
hypersonic  flight  environment.  Section  3  provides  an  overview  of  the  design  process  and  how  it 
is  influenced  by  material  property  uncertainties.  Section  4  describes  general  materials 
information  useful  for  designers  and  analysts.  Section  5  discusses  specific  material  properties 
and  how  they  influence  the  design  and  analysis  of  the  UHTC  components.  Section  6  provides  a 
simple  numerical  example  of  how  property  uncertainties  can  impact  the  design  space.  Finally, 
Section  7  summarizes  the  conclusions  and  provides  recommendations  to  UHTC  material 
developers. 


2.  Summary  of  Hypersonic  Aerothermodynamics 

The  analysis,  design,  and  testing  of  UHTC  components  is  driven  by  the  hypersonic 
aerothermodynamic  environment  in  which  they  are  to  be  used.  Generally,  the  term  “hypersonic” 
refers  to  speeds  in  excess  of  Mach  5,  or  more  than  five  times  the  speed  of  sound  in  the  gas 
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medium.  Vehicles  transiting  an  atmosphere  at  hypersonic  speeds  are  exposed  to  extreme 
convective  heating  from  high-enthalpy  gas  flowing  past  the  vehicle.  In  some  flight  regimes, 
radiative  heating  for  hot  gases  can  also  be  important.  This  section  serves  as  a  brief  introduction 
to  the  fundamental  aspects  of  hypersonic  flow. 


2. 1  Shock  and  boundary  layers 

At  hypersonic  speeds,  a  bow  shock  forms  in  front  of  the  nose  tip  or  WLE  (Fig.  3).  The 
stagnation  point  is  the  location  on  the  leading  edge,  facing  directly  into  the  flow.  The  shock 
standoff  length  defines  the  distance  between  the  bow  shock  and  the  stagnation  point.  Gases 
passing  through  the  bow  shock  are  compressed  and  increase  in  density,  pressure,  static  enthalpy, 
and  temperature.  At  sufficiently  high  enthalpies,  the  gas  may  become  dissociated  and  ionized.  In 
air,  dissociation  means  that  the  N2  and  O2  molecules  break  apart  to  become  highly  reactive  N  and 
O  atoms.  Because  of  their  lower  bond  energy  (5.1  vs.  9.8  eV),  O2  molecules  dissociate  at  lower 
enthalpy  levels  than  N2  molecules.  The  exothermic  recombination  of  these  atoms  at  the  surface 
can  contribute  significantly  to  the  total  aerothermal  heating  experienced  by  a  component, 
depending  on  the  particular  flight  or  test  environment. 


When  the  high-enthalpy  gas  formed  in  the  bow  shock  impinges  on  the  leading  edge,  a 
boundary  layer  forms  along  the  surface  of  the  component  (Fig.  3).  The  boundary  layer  region  is 
characterized  by  large  temperature,  momentum,  and  chemical  composition  gradients  normal  to 
the  surface.  The  gas  flow  is  decelerated  at  the  stagnation  point,  forming  a  thin  boundary  layer 
and  a  local  subsonic  region.  The  boundary  layer  thickens  as  the  flow  is  redirected  and  accelerates 
around  the  WLE  or  nose  tip,  regaining  supersonic  velocities.  Boundary  layer  conditions 
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determine  local  convective  heating  to  the  surface.  Convective  heating  is  typically  highest  in  the 
stagnation  point  region  and  decreases,  often  significantly,  away  from  the  stagnation  point. 


2.2  Surface  energy  balance 

The  wall  temperature  on  a  hypersonic  vehicle  is  determined  by  an  energy  balance  that 
accounts  for  all  of  the  heat  transfer  processes  transporting  energy  into  and  out  of  the  surface. 
Figure  4  is  a  schematic  of  a  steady-state  energy  balance  at  the  surface,  considering  convective 
heating  from  the  boundary  layer,  chemical  heating  from  surface  recombination  reactions,  re¬ 
radiation  into  the  ambient  environment,  and  conduction  into  the  interior  of  the  material. 
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Fig.  4.  Surface  energy  balance. 


The  sensible  convective  heat  flux  to  the  surface,  qconv ,  can  be  defined  as  a  function  of  the  gas 
enthalpy  at  the  edge  of  the  boundary  layer,  he ,  the  enthalpy  of  the  gas  at  the  wall,  hw,  and  a 
transfer  coefficient,  CH .  The  expression  for  convective  heat  flux  can  be  written  as 


<lconv  =  CI{[he(Te)-  K(Tw)\  , 


(1) 


where  Te  is  the  boundary  layer  edge  gas  temperature  and  Tw  is  the  temperature  of  the  wall 
(surface). 

The  transfer  coefficient  is  a  function  of  the  boundary  layer  conditions,  which  are  dependent 
on  the  free  stream  conditions  (speed,  temperature,  and  pressure),  the  gas  composition,  and  the 
geometry  of  the  leading  edge.  The  transfer  coefficient  is  also  a  function  of  whether  the  boundary 
layer  flow  is  laminar  or  turbulent.  Turbulence  can  significantly  increase  the  convective  heating  to 
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the  surface.  Hypersonic  aircraft  are  usually  meticulously  designed  to  avoid  turbulent  flow  and 
prevent  laminar  boundary  layers  from  becoming  turbulent.  The  roughness  of  the  surface  material 
can  be  a  critical  aspect  of  transition  to  turbulence;  this  will  be  addressed  in  Section  5. 

High  levels  of  dissociated  gases  can  be  present  in  flight  and  also  in  the  free  streams  of  high- 
enthalpy  test  facilities,  including  arcjets  and  inductively-coupled  plasma  wind  tunnels,  such  as 
the  Plasmatron  at  the  von  Karman  Institute  for  Fluid  Dynamics  in  Belgium.  The  diffusion  of 
these  species  to  the  surface,  and  their  exothermic  heterogeneous  reaction  on  or  with  the  surface 
material,  can  transfer  additional  energy  to  the  surface.  Gas-surface  chemical  interactions  can  be 
quite  complex,  and  approximate  computational  models  are  almost  always  used.  A  simplified 
model  of  chemical  heating  that  considers  only  independent  oxygen  and  nitrogen  surface 
recombination  can  be  written  as 


Q  chem 


=  r'o(T.) 


A E 


dis,0? 


r o+r'ArJ 


dis,N , 


■  N  ’ 


(2) 


where  y]  is  the  total  catalytic  recombination  efficiency,  E  is  the  surface  impingement  flux  for 
reactant  species  i,  and  A Edis  J  is  the  dissociation  energy  of  product  j.  The  total  catalytic 
recombination  efficiency  ( y[ )  is  the  product  of  the  species  recombination  efficiency,  y ,  and  the 
energy  accommodation  efficiency,  /? .  The  species  recombination  efficiency  ( y )  is  the  fraction 
of  collisions  with  the  surface  that  result  in  loss  of  the  reactant  from  the  gas  phase,  and  the  energy 
accommodation  efficiency  ( /?  )  is  the  fraction  of  exothermic  reaction  energy  that  is  transferred  to 
the  surface.  For  independent  reactions  y' ,  y  and  /?  may  all  vary  from  a  minimum  of  0  to  a 
maximum  of  1 .  The  surface  impingement  flux,  T; ,  is  typically  taken  as  the  gas  kinetic  expression 

ni t  ,  where  R  is  the  universal  gas  constant,  and  ni  and  M t  are  the  number  density 
and  molar  mass  of  the  reactant,  respectively. 

The  re -radiation  away  from  the  surface  is  a  function  of  the  wall  temperature,  the  environment 
temperature,  Tv  ,  the  Stefan-Boltzmann  constant,  a ,  and  the  emissivity,  s  ,  and  can  be  expressed 

by  the  Stefan-Boltzmann  relationship 


(3) 


The  emissivity  may  vary  from  0  to  1 ;  s  =  1  for  a  perfect  blackbody  radiator. 

Energy  conduction  into  the  interior  is  a  function  of  the  local  temperature  gradient  into  the 
surface  and  the  thermal  conductivity,  k  ,  of  the  material  through  Fourier’s  law: 
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(4) 


4co„j  =-k(Tw) 


dT_ 

dx 
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The  steady-state  surface  energy  balance  equates  the  total  aerothermal  heating  of  the  surface 
(qconv  +  qchem )  to  the  energy  transferred  away  from  the  surface  (q  rad  Q  Cond  )  * 


CH[K(T,)-K{T,)]+r'a(T.) 
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dis,02 
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dis ,  AC 


-(T„)Tn  =  c(Tw)c{t:  -T^\-k(T„ 


ylT 

dx 


(5) 


The  surface  energy  balance  can  be  used  to  compute  the  surface  temperature,  given  the 
conditions  of  the  flow,  environment,  and  material  properties.  The  actual  solution  of  the 
expression  is  complicated  by  the  fact  that  the  properties  of  real  materials,  such  as  UHTCs,  are 
functions  of  temperature.  Moreover,  Eq.  (5)  is  a  simplification  of  the  more  general  case,  which 
might  also  include  radiation  transfer  to  the  surface  from  hot  gases  in  the  shock  layer,  more 
complicated  surface  reactions,  as  well  as  energy  storage  effects  associated  with  transient  heating 
conditions.  Generally,  the  solution  is  part  of  a  large  numerical  thermal  analysis  using 
computational  fluid  dynamics  (CFD)  programs  and  involving  multidimensional  geometry  and 
time-dependent  boundary  conditions  that  are  usually  provided  by  an  independent  solution  of  the 
flow  field  and  boundary  layer. 

The  significance  of  Eq.  (5)  for  material  scientists  is  that  the  surface  temperature  experienced 
by  a  UHTC  component  during  hypersonic  flight  or  high-enthalpy  flow  testing  is  determined  not 
only  by  the  environment  and  vehicle  or  component  shape,  but  also  by  temperature-dependent 
material  properties,  such  as  emissivity,  thermal  conductivity,  and  catalytic  efficiency  (as  well  as 
absorbance  and  heat  capacity  in  the  more  general  case). 


2.3  Approximate  analyses  and  estimates 

Because  of  the  coupled  nature  of  the  surface  energy  balance,  complicated  numerical  analyses 
are  typically  required  for  accurate  solutions.  However,  two  approximate  analyses  are  in  common 
use,  particularly  in  the  early  stages  of  vehicle  or  mission  design. 

Hypersonic  aerothermal  environments  are  often  quantified  in  terms  of  a  heat  flux  to  the 
vehicle  surface.  For  example,  vehicle  designers  may  specify  that  a  candidate  leading-edge 
material  must  sustain  a  particular  peak  stagnation-point  heating  rate.  The  transfer  coefficient  and 
boundary  layer  edge  temperature  are  relatively  insensitive  to  the  thermal  conditions  at  the  wall 
and  are  often  predicted  independent  of  the  wall  conditions,  using  CFD  codes  or  engineering 
correlations,  together  with  information  about  the  tree-stream  conditions  and  the  nose  tip  or  WLE 
geometry.  One  common  approximation  for  the  aerothermal  heat  flux  is  to  ignore  the  surface 


A4-7 


energy  balance  and  calculate  heating  directly  from  the  left  hand  side  of  Eq.  (5),  assuming  a  fixed, 
low-wall  temperature,  typically  300  K,  and  fully  catalytic  behavior  (/'  =  l). 

The  heat  flux  defined  in  this  way  is  referred  to  as  qfccw ,  a  “fully  catalytic  cold-wall”  heat  flux. 

The  advantage  of  this  definition  is  that  q  fc  cw  is  only  a  function  of  the  flow  conditions  (free 

stream  and  geometry)  and  is  completely  independent  of  the  surface  thermal  conditions.  It 
provides  a  common  reference  point  for  the  comparison  of  different  flow  environments.  In  a 
sense,  the  fully  catalytic  cold-wall  heat  flux  represents  the  energy  in  the  environment  that  is 
potentially  available  to  heat  the  surface.  Usually,  when  aerospace  engineers  or 
aerothermodynamicists  refer  to  the  “peak  heat  flux”  in  a  hypersonic  flight  regime,  they  are 
referencing  to  a  fully  catalytic,  cold-wall  value.  The  fully  catalytic,  cold-wall  heat  flux  is  always 
higher  than  the  corresponding  hot-wall  flux  in  flight,  as  surface  temperatures  in  flight  will 
typically  be  much  hotter  than  300  K,  and  real  materials  are  not  fully  catalytic. 

The  distinction  between  cold-wall  and  hot-wall  heat  flux  is  also  important  when  evaluating 
material  performance  in  ground-test  environments.  For  example,  test  conditions  in  arcjet  or 
Plasmatron  wind  tunnels  are  routinely  characterized  using  water-cooled  copper  calorimeters, 
whose  surface  temperatures  remain  near  room  temperature  during  the  measurement.19,20  Such  a 
measurement  is  essentially  a  cold-wall  heat  flux  measurement  on  a  highly  catalytic  surface.  The 
heat  flux  experienced  by  the  test  specimen  under  the  identical  flow  conditions  will  be  a  hot-wall 
heat  flux,  on  a  surface  that  is  probably  much  less  catalytic  than  copper. 

A  second  common  approximation,  useful  for  quickly  estimating  thermal  response  to  a 
prescribed  heat  flux,  is  the  “radiation  equilibrium  temperature.”  This  approximation  simplifies 
the  surface  energy  balance  by  assuming  negligible  heat  conduction  into  the  material  ( qcond  =  0 ) 

and  a  constant  value  for  aerothermal  heat  flux,  often  a  fully  catalytic,  cold-wall  heat  flux  value 
obtained  as  described  above.  With  these  assumptions,  the  energy  balance  can  be  rearranged  as  an 
algebraic  equation  for  Tw : 


T  =  4 


so 


■  +  TZ 


(6) 


Since  Eq.  (6)  ignores  the  effect  of  heat  conduction  into  the  interior,  the  radiation  equilibrium 
temperature  generally  represents  the  maximum  potential  surface  temperature  for  a  given  heat 
flux  condition,  assuming  a  conservative  (low)  value  of  surface  emissivity  is  used.  The  radiation 
equilibrium  temperature  can  be  a  valuable  engineering  level  tool  for  quickly  estimating  thermal 
response  and  as  a  sanity  check  for  validating  computational  analyses.  However,  the 
approximation  of  zero  heat  conduction  is  generally  not  valid  for  materials  with  relatively  high 
thermal  conductivities,  such  as  metals  or  UHTCs.  As  discussed  previously,  high  thermal 
conductivity  is  one  of  the  enabling  properties  of  diboride-based  UHTC  materials,  because  it 
allows  heat  to  be  drawn  out  of  the  stagnation  region  and  re-radiated  from  a  larger  surface  area. 
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Thus,  a  radiative  equilibrium  temperature  may  significantly  overestimate  the  steady-state  surface 
temperature  in  the  stagnation  point  region  of  a  sharp  UHTC  nose  tip  or  WLE. 


2.4  Computational  methods 

Predicting  the  aerothermal  conditions  and  material  thermal  response  of  a  UHTC  component  in 
a  hypersonic  vehicle  design  usually  requires  the  application  of  numerical  computations.  Even  for 
simple  sphere-cone  geometries,  the  flow  conditions  and  thermal  response  are  too  complex  to 
evaluate  with  direct  analytical  methods.  In  a  real  aerothermal  flight  environment,  the  thermal 
conditions  at  the  material  surface  and  the  boundary  layer  conditions  are  interdependent. 

However,  for  dimensionally  stable  (non-ablating,  non-pyrolyzing)  thermal  protection  systems, 
this  interaction  is  often  weak,  when  compared  to  the  effect  that  the  free-stream  flow  conditions 
and  vehicle  shape  have  on  the  boundary  layer.  Under  these  conditions,  it  is  preferable  to  perform 
design  computations  in  which  the  numerical  analyses  of  the  flow  conditions  and  the  material 
thermal  response  are  performed  independently  —  and  in  which  the  interface  between  the  fluid 
and  solid  regions  is  treated  as  a  boundary  condition  in  each  of  the  independent  analyses. 

NASA  uses  computational  fluid  dynamics  programs,  such  as  the  Data  Parallel  Line 
Relaxation  (DPLR)  program  or  the  Langley  Aerothermodynamic  Upwind  Relaxation 
Algorithm  (LAURA),22  to  predict  the  flow  conditions  around  a  vehicle  traveling  a  particular 
flight  trajectory  through  the  atmosphere.  These  programs  solve  the  Navier-Stokes  equations  and 
can  predict  the  dynamic,  thermodynamic,  and  chemical  conditions  in  the  shock  layer  and 
boundary  layer,  as  well  as  provide  the  surface  heat  flux  as  a  boundary  condition  to  the  materials 
response  model.  Material  response  analyses  employ  finite-difference,  finite-volume,  or  finite- 
element  methods  to  perform  computations  as  a  function  of  transient  heating.  NASA  Ames 
Research  Center  uses  the  commercial  finite-element  MSC.Marc  package23,24  to  predict  the 
transient  surface  and  in-depth  thermal  response  of  UHTC  nose  tip  and  WLE  components.  Other 
researchers  have  used  the  commercial  codes  ABAQUS,25  COMOS/M,  COSTAR,26  and 
ANSYS  ’  to  model  UHTC  thermal  and  mechanical  response. 

The  complexity  associated  with  the  prediction  of  material  response  along  a  flight  trajectory  is 
highly  dependent  on  the  level  of  coupling  between  the  flow  field  and  material  analyses. 
Decoupled  solutions  are  obtained  by  setting  surface  temperatures  to  either  cold-wall  or  radiative 
equilibrium  values  in  the  CL D  analyses  and  then  using  the  computed  transient  convective 
heating  profile  as  a  boundary  condition  to  a  material  response  model.  In  this  approach,  no 
material  properties  are  required  for  the  flow  field  analysis,  if  the  surface  temperature  is  set  to  an 
arbitrary  cold-wall  value,  and  only  the  surface  emissivity  is  required  if  the  radiative  equilibrium 
temperature  assumption  is  used.  Although  the  radiation  equilibrium  temperature  generally  over 
predicts  the  actual  surface  temperature,  the  estimate  is  usually  much  closer  than  assuming  a  cold 
wall. 

However,  for  sharp  UHTC  components,  the  combination  of  high  thermal  conductivity,  multi¬ 
dimensional  heat  conduction,  and  significant  volumetric  heat  capacity,  often  necessitates  some 
level  of  coupling  between  the  flow  field  and  material  response  computations,  in  order  to  achieve 
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time-accurate  predictions  for  temperature  and  stress  in  a  transient  aerothermal  heating 
environment.  Loose  coupling  is  achieved  by  incrementing  flow  field  and  material  response 
computations  sequentially  during  each  time  step,  while  full  coupling  requires  the  iterative 
solution  of  the  flow  field  and  material  responses  to  satisfy  the  full  surface-energy  balance  during 
every  time  step.  For  trajectory-based  design  space  studies,  solution  of  the  full  Navier-Stokes 
equations  is  sometimes  replaced  by  more  approximate  analytic-flow  correlations  and  boundary- 
layer  solutions  to  make  the  computations  more  tractable.  Requirements  for  the  level  of 
coupling  to  achieve  time-accurate  solutions  for  UHTC  components  are  discussed  in  the 
literature.24,29,30 


3.  The  Design  Process  and  Property  Uncertainty 

Hypersonic  vehicle  designs  are  driven  by  mission  objectives  that  demand  certain  performance 
criteria.  Designers  attempt  to  fulfill  these  criteria  in  a  constrained  optimization  process, 
adjusting  the  shape  and  mass  of  a  vehicle  to  achieve  the  desired  flight  characteristics,  while 
simultaneously  remaining  within  bounds  imposed  by  the  properties  of  available  aerospace 
materials.  Material  property  uncertainties,  along  with  the  variations  or  dispersions  in  the  other 
parameters,  such  as  aerothermal  heating  or  aerodynamic  loads,  are  used  to  define  the  design 
space,  assess  the  design  against  system  requirements,  and  determine  margins  or  factors  of  safety. 
Components  and  systems  must  “buy”  their  way  onto  a  vehicle  design  by  proving  that  they  are 
lower  risks  than  competing  systems.  Reducing  uncertainties  in  material  properties  reduces  the 
risks  in  the  design  —  risks  associated  with  component  failure,  mission  failure,  or  loss  of  the 
vehicle.  The  aerospace  industry  maintains  standards  that  define  the  design  process  and  the  role 
that  material  properties  play  in  that  process.  NASA  has  developed  its  own  standards  for  the 
design  and  assessment  of  spaceflight  hardware  and  thermal  protection  system  design.31"34 

An  important  early  step  in  the  design  process  is  defining  and  selecting  the  specific  set  of 
material  property  values  that  will  be  used  in  analyses.  The  design  team  must  evaluate  different 
data  sources  and  select  the  data  sets  deemed  most  directly  applicable  and  reliable  for  the  given 
application.  Often  this  step  includes  down  selection  from  a  large  property  database,  like  NASA’s 

'JT 

TPSX  Material  Properties  Database.  The  availability  and  quality  of  property  uncertainties 
associated  with  different  data  sets  is  always  a  key  criterion  for  this  down-selection  process.  At 
this  stage,  the  lack  of  quality  uncertainty  information  or  specific  experimental  data  often  dictates 
that  new  property  measurements  be  made. 

Material  property  data  are  usually  provided  as  either  compiled  statistical  values,  such  as  a 
mean  and  standard  deviation  on  a  particular  property,  or  as  a  set  of  data  points  from  which  the 
statistical  values  can  be  derived.  It  is  typical  for  designers  to  use  bounding  values  of  ±3  standard 
deviations  (3a)  to  define  the  dispersions  in  the  design.  A  large  variation  in  a  measured  property 
may  mean  that  the  3a  values  push  some  design  parameters  beyond  acceptable  limits.  Some 
specific  approaches  for  the  use  of  property  uncertainty  data  include  performing  sensitivity 
studies,  in  which  material  properties  are  parameters  that  are  varied  in  a  regular  pattern  to 
ascertain  their  effect  on  the  thermal  and  mechanical  response  of  the  system.  A  common 
application  of  this  approach  is  to  assess  the  system  with  a  “worst-on-worsf  ’  combination  of 
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material  properties  and  flight  loads.  Another  approach  uses  a  Monte  Carlo  simulation  to  run  a 
large  number  of  analyses,  randomly  varying  material  properties  and  other  design  parameters 

O /T 

based  on  their  nominal  values  and  standard  deviations.  These  analyses  are  very  effective  for 
identifying  the  parameters  that  have  the  largest  influence  on  the  system  performance. 


4.  General  Material  Information 

Before  proceeding  to  a  discussion  of  specific  material  properties  used  in  design  computations, 
it  is  helpful  to  outline  the  general  material  information  that  designers  and  analysts  need  when 
developing  hypersonic  vehicles.  This  information  is  useful  in  assessing  the  results  of  analyses  or 
in  determining  the  potential  applicability  of  a  material. 


4. 1  Sources  of  uncertainty 

Multiple  tests  for  the  same  material  property  usually  yield  a  distribution  of  values,  even  for 
relatively  easy-to-measure  properties  such  as  density.  Designers  need  to  know  these 
distributions,  in  order  to  define  bounding  cases  and  dispersions  in  the  design  space.  Distributions 
in  material  properties  arise  from  a  variety  of  sources,  including  manufacturing  variability,  the  use 
of  different  techniques  to  derive  the  same  material  property,  and  inherent  uncertainties  in 
different  measurements  techniques. 

No  material  manufacturing  process  is  perfectly  reproducible.  Lot-to-lot  variations  are 
expected  and  tolerated  (within  certain  limits).  IJHTC  materials  are  certainly  no  exception  to  this 
rule.  Manufacturing  processes,  such  as  hot  pressing  of  ceramic  powders,  are  defined  by  a  large 
number  of  variables  that  can  only  be  controlled  within  certain  tolerances.  Even  under  nominal 
processing  conditions,  acceptable  material  lots  will  have  properties  that  vary  about  some  average 
value.  It  is  also  common  for  materials  manufactured  in  a  pressing,  molding,  or  casting  process  to 
exhibit  non-uniform  material  properties  within  a  single  billet.  Such  non-uniformities  can  occur, 
for  example,  near  surfaces,  due  to  interactions  with  die  or  mold  surfaces,  or  along  directions 
aligned  with  the  application  of  stresses  or  gravity  that  cause  segregation  of  constituents.  It  is 
important  for  materials  researchers  to  quantify  both  inter-billet  and  intra-billet  property 
distributions. 

Whenever  possible,  multiple  material  property  measurements  should  be  made  on  specimens 
from  the  same  billet  or  manufactured  from  the  same  lot  of  starting  material.  Such  data  are  more 
valuable  than  properties  measured  on  specimens  from  different  starting  billets  or  material  lots, 
because  they  allow  material  developers  and  vehicle  designers  to  determine  correlations  among 
material  properties  and  to  estimate  trends  in  manufacturing. 

Consistency  in  test  specimens  is  also  important  for  any  material  property  derived  from 
combinations  of  separate  test  data.  For  example,  the  thermal  conductivity  is  related  to  the 
density,  p ,  specific  heat,  cp ,  and  thermal  diffusivity,  a  ,  by  the  relationship 

k  =  a pcp  .  (7) 
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Thermal  diffusivity  is  often  measured  directly  at  high  temperature  using  laser-flash  techniques, 
but  it  is  usually  thermal  conductivity  that  is  required  as  an  input  to  thermal  analysis  programs. 
The  thermal  conductivity  must  then  be  derived  from  the  measurements  of  density,  specific  heat, 
and  diffusivity,  and  if  these  properties  are  measured  on  unrelated  samples,  the  calculation  of 
thermal  conductivity  may  result  in  inconsistent  values  at  best,  or,  even  worse,  significantly 
inaccurate  values. 

Many  material  properties  can  be  measured  with  a  variety  of  experimental  techniques.  Often 
these  techniques  and  their  associated  data  reduction  procedures  are  quite  different.  For  example, 
tensile  modulus  can  be  measured  by  acoustic  techniques  or  mechanical  tests;  fracture  strength 
can  be  measured  in  tension,  compression,  or  bending  geometries;  thermal  conductivity  can  be 
obtained  directly  using  steady-state  methods  like  guarded  hot-plate  testing,  or  indirectly,  by 
transient  methods  like  laser- flash  or  photothermal  radiometry,  which  measure  thermal 
diffusivity. 

It  is  important  to  describe  or  reference  measurement  methods  well,  so  that  informed 
judgments  can  be  made  and  discrepancies  among  different  property  values  resolved.  In  this 
regard,  standard  measurement  techniques,  such  as  those  defined  by  the  American  Society  for 
Testing  and  Materials  (ASTM),  are  desirable,  because  both  the  measurement  procedure  and 
uncertainty  are  defined  in  the  standard.  The  use  of  nonstandard  techniques  should  be  avoided, 
whenever  possible,  because  uncertainty  is  more  difficult  to  assess. 

In  summary,  the  more  information  supplied  to  the  designer  about  how  a  property  was 
measured  and  what  sources  of  uncertainty  are  incorporated  into  quoted  errors,  the  tighter  the 
resulting  design  space.  Missing  information  inevitably  forces  analysts  to  assume  conservative 
uncertainties  that  may  vastly  exceed  those  that  actually  exist. 


4.2  Material  morphology 

UHTC  materials  are  often  referred  to  as  monolithic  and  are  treated  as  having  homogeneous 
properties.  But  a  micrograph  of  an  HfB2-SiC  UHTC  sample  (Fig.  5)  clearly  shows  the  composite 
nature  of  the  material.  Microstructure  significantly  influences  the  properties  of  a  material 
through  the  distribution  of  different  material  phases  and  the  location  and  size  of  defects  and 
voids.  Information  about  material  microstructure  can  help  analysts  interpret  the  macroscopic 
performance  of  the  material  and  is  particularly  useful  in  understanding  and  predicting  failure 
modes. 
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Fig.  5.  Scanning  electron  micrograph  of  a  HfB2  UHTC  composite  containing  20  volume 
percent  SiC;  the  light  grains  are  HfB2  and  the  dark  grains  are  SiC. 


An  important  aspect  of  material  morphology  is  the  introduction  of  anisotropic  material 
properties.  Many  aerospace  materials,  such  as  carbon  fiber  composites  or  rigid  ceramic  tile 
insulations,  exhibit  direction-dependent  thermal  and  mechanical  properties.  Generally,  these 
materials  fall  into  two  categories:  orthotropic  materials  with  properties  oriented  in  three  mutually 
orthogonal  directions,  and  transverse  isotropic  materials,  with  properties  oriented  in  two 
orthogonal  directions,  though-the-thickness  and  in-plane.  Most  hot-pressed  UHTCs  are  treated  as 
isotropic  in  analysis  and  design,  but  it  is  important  for  material  developers  to  investigate  the 
directional  dependence  and  verify  that  the  isotropic  assumption  is  valid.  Failure  to  incorporate 
significant  orthotropic  or  transverse  isotropic  properties  in  analyses  can  lead  to  erroneous 
predictions  of  thermal  and  mechanical  performance  and,  ultimately,  poor  component  design. 


4.3  Material  use  limits 

Some  material  parameters  are  not  measured  directly,  but  are  determined  by  material 
performance.  Maximum  temperature  limits  and  maximum  stress  limits  are  examples  of  such 
parameters.  Designers  need  to  establish  these  limits  based  on  an  analysis  of  a  collection  of 
material  properties  and  behaviors. 

Obviously,  the  extreme  upper  temperature  limit  is  set  by  a  material’s  melting  point  or  eutectic 
temperature.  However,  a  hypersonic  vehicle  designer  will  never  expect  to  use  a  material  up  to 
that  temperature.  Instead,  the  designer  needs  to  define  a  practical  temperature  limit  beyond 
which  the  thermal,  mechanical,  and  other  properties  are  potentially  degraded  to  such  a  point  that 
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the  material  can  no  longer  perform  its  required  function.  The  material  composition,  although  not 
used  directly  in  design  computations,  can  be  useful  information  in  defining  capabilities  and 
material  limits.  For  example,  if  any  of  the  material  components  are  susceptible  to  oxidation,  it  is 
important  to  know  the  temperature  at  which  oxidation  rates  become  significant.  Such 
information  would  not  be  used  directly  in  a  computational  analysis,  but  may  affect  the 
operational  limits  of  the  vehicle  or  may  demand  that  components  be  coated  with  an  oxidation 
barrier. 

The  mechanical  performance  constraints  may  be  associated  with  more  than  one  type  of  stress 
limit,  depending  on  the  mode  of  loading:  tensile,  compressive,  or  shear.  Maximum  stress  limits 
are  usually  obtained  from  mechanical  failure  data.  While  the  thermal  and  electrical  properties  of 
UHTC  materials  are  somewhat  “metallic”  in  nature,  the  failure  characteristics  are  more  typical  of 
the  brittle  fracture  of  a  ceramic.  Most  fracture  strength  data  on  UHTC  materials  are  obtained  in 
flexure  testing.  ’  ’  '  UHTC  fracture  strengths  tend  to  decrease  with  increasing  temperature. 

Because  it  is  difficult  to  perform  failure-prediction  analyses  of  brittle  materials,  analysts 
generally  perform  standard  mechanical  or  thermal/mechanical  analyses  and  then  compare  the 
predicted  stress  results  with  experimentally  measured  mechanical  strength  values.  For  brittle 
materials,  this  may  be  in  the  form  of  statistical  data,  such  as  Weibull12,43  Designers  will  usually 
apply  some  factor  of  safety  to  the  strength  values  and  use  those  margined  values  for  comparison 
with  mechanical  analyses.  If  mechanical  strength  data  are  not  accurate  or  extensive  enough  to 
generate  robust  statistics,  the  designers  may  be  required  to  apply  a  higher  factor  of  safety  to 
cover  the  uncertainty,  which  can  lead  to  overly  conservative  designs. 

In  general,  material  developers  need  to  provide  enough  accurate  material  property  data  under 
sufficiently  stressful  testing  conditions  so  that  proper  use  limits  can  be  established  to  within  an 
acceptable  uncertainty.  Testing  UHTC  materials  and  components  at  very  high  temperatures  and 
past  the  point  of  failure  are  import  aspects  of  material  characterization. 


4.4  Properties  not  used  in  component  analyses 

There  are  many  properties  commonly  reported  for  UHTC  materials  that  are  not  directly  useful 
for  thermal  or  mechanical  analyses.  Some  examples  are  hardness  or  microhardness5’12’44’45 
fracture  toughness,10’39’41’46'49  and  thermal  shock  parameters.5,1 1,50  These  properties  can  be 
valuable  indicators  for  tracking  processing-related  changes  and  for  optimizing  desired  properties 
or  capabilities  in  UHTC  materials.  They  can  also  be  helpful  when  screening  candidate  materials 
for  particular  applications  and  for  assessing  performance  issues,  such  as  damage  resistance 
during  handling  or  shipping.  However,  properties  like  hardness,  fracture  toughness,  and  thermal 
shock  resistance,  are  of  marginal  importance  in  design  computations,  as  their  numerical  values 
are  never  input  directly  into  any  analyses. 
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5.  Specific  Material  Properties 

The  design  and  analysis  of  UHTC  components  for  hypersonic  applications  ultimately  require 
numerical  values  of  specific  physical,  thermal,  mechanical,  and  surface  properties.  The 
important  material  properties  discussed  below  appear  directly  in  either  boundary  conditions  like 
Eq.  (5)  or  in  governing  equations  for  heat  transfer  and  mechanical  equilibrium.  Some  examples 
of  governing  equations  for  an  isotropic  material  are  the  energy  conservation  equation 
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which  relates  the  energy  storage  rate  to  heat  flux  gradients  within  a  material,  and  Hooke’s  law, 
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which  relates  the  elastic  strains,  stj ,  and  stresses,  <x(/ ,  within  a  material  through  the  tensile 
modulus,  E,  and  the  Poisson’s  ratio,  v. 

Because  UHTC  materials  in  hypersonic  vehicle  applications  are  required  to  operate  at  very 
high  temperatures,  designers  and  analysts  will  require  that  most  of  the  properties  discussed 
below  be  provided  as  a  function  of  temperature  over  the  entire  temperature  range  of  interest. 
Measuring  material  properties  at  high  temperatures  is  always  difficult  and  sometimes  impossible. 
In  such  cases,  material  developers  should  provide  estimates  of  high-temperature  properties  based 
on  engineering  judgment,  similarity  with  other  materials,  or  extrapolation  of  existing  data. 
Whatever  approach  is  used,  an  estimate  of  the  uncertainty  or  error  in  the  high-temperature 
properties  should  also  be  provided. 


5. 1  Physical  properties 

The  most  important  physical  properties  for  UHTC  materials  are  the  density  and  the  coefficient 
of  thermal  expansion  (CTE). 

Density  appears  as  a  first  order  term  in  the  governing  equations  for  thermal  and  dynamic 
mechanical  response  (for  example,  vibration).  Because  variations  in  density  have  a  large  effect 
on  the  thermal  and  dynamic  mechanical  performance  of  components,  material  developers  should 
provide  as  accurate  a  value  as  possible  with  as  tight  a  tolerance  as  possible.  It  should  be  possible 
to  make  individual  density  measurements  with  measurement  errors  less  than  1%.  Density  also 
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determines  the  mass  (weight)  of  a  component,  which  is  always  a  critical  aspect  of  the  design. 
IJHTC  materials  are  much  denser  than  typical  aerospace  materials,  such  as  carbon  composites  or 
aluminum  alloys,  and  managing  the  total  and  distributed  mass  of  a  vehicle  is  one  of  the  most 
important  functions  of  the  designer.  Typical  designs  for  UHTC  wing  leading  edges  and  nose  tips 
require  very  little  material,  due  to  the  high  temperature  capability  and  high  thermal  conductivity. 
Additionally,  the  location  of  the  high  density  UHTCs  at  the  forward  end  of  the  vehicle  can  have 
the  beneficial  effect  of  moving  the  center  of  gravity  further  forward,  which  often  results  in  better 
aerodynamic  stability. 

Most  materials  exhibit  changes  in  density  as  a  function  of  temperature  associated  with  thermal 
expansion.  The  effects  of  density  changes  on  heat  transfer  performance  are  usually  negligible.  In 
the  analysis  of  thermo-mechanical  responses,  the  thermal  expansion  behavior  of  a  material  is 
typically  incorporated  in  the  linear  coefficient  of  thermal  expansion  (CTE).  The  CTE  is  a  critical 
material  property.  Local  thermal  expansion  within  a  UHTC  component  can  produce  large 
internal  strains  and  stresses,  particularly  under  conditions  that  cause  rapid  temperature  changes. 
Differential  thermal  expansion  between  a  UHTC  component  and  a  neighboring  material  can 
cause  undesirable  contact  stresses.  Thermally  induced  stresses  often  limit  the  size  and  shape  of  a 
UHTC  component.  When  predicted  thermal  stresses  approach  maximum  allowable  stress  limits, 
conditions  for  thermal  shock  failure  exist  and  components  may  have  to  be  redesigned  to  dissipate 
heat  more  effectively.  This  generally  means  making  smaller  UHTC  components  with  lower 
aspect  ratios.  For  non-isotropic  materials,  the  CTE  must  be  provided  in  each  of  the  primary 
material  directions. 

The  CTE  is  usually  derived  from  measurements  of  the  elongation  of  a  material  at  elevated 
temperatures.  Because  the  CTE  is  the  first  derivative,  or  slope,  of  the  curve  of  elongation  versus 
temperature,  measurements  must  be  performed  with  sufficient  resolution  to  extract  this 
derivative,  without  introducing  undesirable  artifacts  into  the  CTE  values.  The  CTE  can  also  be 
obtained  for  pure  crystalline  materials  using  X-ray  diffraction  techniques  to  measure  changes  in 
lattice  constants  with  temperature.  Often  the  thermal  expansion  behavior  over  a  large 
temperature  range  is  non-linear  and  must  be  defined  by  several  average  CTE  values  over 
sequential  temperature  intervals.  Typical  CTE  values  for  ZrB2  and  HfB2  composites  lie  in  the 
5-8  xlO"6  K'1  range,  with  the  higher  values  found  at  higher  temperatures. 11,12,48’49’51 


5.2  Thermal  properties 

The  most  important  thermal  properties  for  UHTC  materials  are  the  specific  heat  and  the 
thermal  conductivity.  Both  properties  appear  as  first-order  terms  in  the  governing  conservation- 
of-energy  equation. 

The  specific  heat  can  have  a  large  effect  on  the  transient  thermal  response  of  a  UHTC 
component  during  heating  or  cooling.  Specific  heat  is  generally  a  strong  function  of  temperature 
for  UHTC  materials,  particularly  from  room  temperature  up  to  -1000  °C.  For  example,  the 
specific  heat  of  HfB2  increases  by  about  a  factor  of  two  over  this  temperature  range.  While 
specific  heat  values  for  composite  materials  can  be  estimated  from  measurements  performed  on 
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isolated  constituents,  values  measured  directly  on  the  actual  material  are  much  preferred.  Sample 
volumes  used  in  calorimetry  are  often  small.  For  UHTC  composites,  it  is  important  to  ensure  that 
the  samples  used  for  specific  heat  measurements  accurately  reflect  the  composition  of  the  bulk 
material.  During  high  temperature  calorimetry,  the  absence  of  chemical  reactions  between 
UHTC  constituents  and  sample  containers  must  be  verified. 

Thermal  conductivity  exerts  a  dominant  effect  on  heat  transfer  and  coupled  thermo¬ 
mechanical  responses  and  is  a  particularly  important  property  for  sharp  UHTC  leading  edge 
components,  as  the  ability  to  conduct  heat  away  from  sharp  edges  is  an  enabling  characteristic. 

In  general,  thermal  conductivity  can  be  directionally  dependent.  If  a  material  exhibits  anisotropic 
behavior,  analysts  will  require  thermal  conductivity  values  for  all  of  the  primary  material 
orientations. 

UHTC  materials  exhibit  thermal  conductivities  with  a  wide  range  of  magnitudes  and  a  variety 
of  temperature  dependencies.  ’  Because  the  thermal  conductivity  of  UHTC  composites 
depends  on  many  different  factors,  including  microstructure,  composition,  defect  and  impurity 
levels,  contact  resistance  between  grains,  and  porosity,  values  cannot  be  estimated  with  any 
confidence  from  published  data  or  measurements  performed  on  individual  constituents.  Even 
nominally  similar  UHTC  composites  may  not  have  the  same  heat  transport  characteristics, 
because  the  factors  that  affect  thermal  conductivity  are  largely  determined  by  manufacturing 
details,  for  example,  the  raw  ceramic  powders  sources,  powder  processing  procedures,  and  hot- 
pressing  schedules.12,52  Thus  thermal  conductivities  over  the  entire  temperature  range  of  interest 
must  be  measured  for  the  exact  UHTC  material  to  be  used  in  a  particular  component  application. 

Since  most  high-temperature  thermal  conductivity  values  are  actually  derived  from  thermal 
diffusivity  measurements  (as  previously  described),  it  is  important  that  the  specific  heat  and 
density  values  required  to  extract  thermal  conductivity  from  thermal  diffusivity  be  measured  on 
the  same  lot  of  UHTC  material,  preferably  the  same  test  specimen,  if  possible.  Typical 
uncertainties  reported  for  direct  UHTC  thermal  property  measurements  are  ±3%  for  specific 
heat,13,53,54  ±5%  for  thermal  diffusivity,13’53  and  ±5%  for  thermal  conductivity.55  These 
uncertainties  are  associated  with  measurement  technique  and  the  consistency  of  repeated 
measurements  on  the  same  sample  and  do  not  include  multiple  sample  statistics.  Note  also  that  if 
thermal  conductivity  is  calculated  from  k  =  a  pc  ,  uncertainties  in  the  contributing  properties 

combine  and  propagate. 


5.3  Mechanical  properties 

The  primary  mechanical  properties  required  for  analyses  of  isotropic  UHTC  components  are 
the  tensile  modulus  and  the  Poisson’s  ratio.  Both  appear  as  first  order  terms  in  the  governing 
static  equations.  For  isotropic  materials,  the  shear  modulus,  G  ,  can  be  calculated  from  the 
tensile  modulus  and  Poisson’s  ratio  as  G  =  E /[2(l  +  v)]  and  does  not  have  to  be  independently 

specified. 
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The  tensile  modulus  (also  called  the  Young’s  modulus  or  elastic  modulus)  has  a  large  effect 
on  the  mechanical  response  of  the  component.  Ideally,  tensile  modulus  values  are  derived  from 
the  slope  (or  first  derivative)  of  a  material’s  stress-strain  curve  under  uniaxial  tension  in  the 
elastic  loading  regime.  However  UHTC  materials  have  very  limited  ductility,  and  their  tensile 
moduli  are  almost  always  either  measured  by  acoustic  methods  or  derived  from  experimental 
stress-strain  curves  obtained  during  flexure  testing.  The  flexure  method  was  found  to  produce 
somewhat  lower  (by  ~15%)  values  for  some  UHTC  materials  than  the  tensile  and  acoustic 
methods.3,56 

Tensile  modulus  values  can  be  used  directly  in  mechanical  analyses,  but  it  is  often  more 
accurate  to  incorporate  the  actual  stress-strain  relationship.  Most  analytical  software  packages 
can  handle  inputs  of  either  type.  Incorporating  the  experimental  stress-strain  curve  is  particularly 
advantageous  for  materials  that  exhibit  non-linear  elastic  response  or  for  when  the  experimental 
stress-strain  curve  contains  plastic  strain  contributions  deemed  important  for  a  particular 
application.  Nonlinear,  elastic  deformation  behaviors  are  observed  for  some  UHTC  materials  at 
elevated  temperatures.3,40'42,57  The  tensile  modulus  magnitude  tends  to  decrease  with  increasing 
temperature.  Accurate  temperature-dependent  modulus  data  are  desirable,  but  difficult  to 
measure.  In  practice,  if  the  expected  thermal  and  load  environments  combine  to  push  a  UHTC 
component  into  a  highly  non-linear  or  highly  temperature-dependent  mechanical  property 
regime,  a  redesign  of  the  component  or  system  is  likely. 

The  Poisson’s  ratio  is  defined  as  the  ratio  of  lateral  contraction  to  axial  extension  under 
uniaxial  loading.  It  can  be  measured  during  a  uniaxial  tensile  test  or  by  acoustic  methods. 
Although  the  Poisson’s  ratio  does  not  often  have  a  large  influence  on  the  magnitude  of  predicted 
stresses  and  strains  in  a  component,  it  can  significantly  affect  the  distribution  and  orientation  of 
those  stresses  and  strains.  This  is  particularly  important  when  modeling  three-dimensional  sharp 
components  subjected  to  large  temperature  gradients.  Typical  values  for  metals  are  around  0.3; 
UHTC  materials  generally  have  lower  Poisson’s  ratios,  in  the  0. 1-0.2  range  at  room 
temperature.10,11,39  The  Poisson’s  ratio  does  not  tend  to  be  highly  temperature  dependent  for  real 
materials,  so  analysts  do  not  usually  require  temperature  dependent  values. 

The  number  of  mechanical  properties  that  must  be  measured  increases  rapidly  with 
anisotropy.  For  transverse  isotropic  materials,  two  tensile  moduli,  two  Poisson’s  ratios,  and  one 
shear  modulus  are  required,  and  for  orthotropic  materials,  three  tensile  moduli,  three  Poisson’s 
ratios,  and  three  shear  moduli  are  required.  As  with  the  tensile  modulus,  it  is  desirable  to  define 
the  temperature-dependent  behavior  of  the  shear  modulus,  but  estimates  of  that  behavior  are 
often  good  enough. 


5.4  Surface  properties 

The  three  most  important  surface  properties  for  UHTC  materials  are  the  emissivity  (or 
emittance),  the  catalytic  efficiency,  and  the  surface  roughness.  The  first  two  properties  enter  into 
the  surface  energy  balance,  while  the  third  influences  laminar-to-turbulent  transitions  in  the 
boundary  layer  flow. 
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In  the  general  case,  emissivity  is  a  function  of  wavelength,  emission  direction,  and 
temperature.  The  total  hemispherical  emissivity  —  a  value  averaged  over  all  wavelengths  and 
emission  directions  —  is  required  in  the  surface  energy  balance.  Because  emissivity  acts  as  a 
multiplicative  coefficient  to  the  Stephan-Boltzmann  radiation  function,  it  has  a  large  effect  on  the 
ability  of  the  surface  to  reject  heat,  and  thus,  on  the  steady-state  surface  temperature  reached  by  a 
component.  The  total  hemispherical  emissivity  can  be  a  strong  function  of  temperature;  so 
temperature-dependent  emissivity  values  are  usually  required.  Direct  measurements  of 
hemispherical  emissivity  at  elevated  temperatures  are  difficult,  and  few  such  data  are  available  in 
the  literature  for  UHTC  materials.  The  measurements  of  Scatteia  et  al.58,59  for 
ZrB2-15SiC-2MoSi2,  ZrB2-15SiC,  and  ZrB2-15SiC-10HfB2  composites  are  a  notable  exception. 
Emissivity  values  ranged  from  0.49  to  0.81,  depending  on  material  composition,  temperature, 
machining  method,  and  ambient  oxygen  environment.  Scatteia  et  al.58,59  assign  errors  of  ±5%  to 
their  measurements. 

Often  temperature-dependent  emissivity  values  are  estimated  from  room-temperature 
hemispherical  spectral  reflectance  measurements.60  We  use  optical  relations  for  an  opaque  solid 
and  Kirchhoff  s  law  to  obtain  spectral  emissivity  from  the  measured  spectral  reflectance  and  then 
average  the  spectral  emissivity  over  the  Planck  function  at  different  temperatures  to  compute 
temperature-dependent  total  emissivity  values.  It  is  difficult  to  assign  rigorous  uncertainties  to 
this  process. 

An  in  situ  emissivity  measurement  can  be  made  during  testing  in  high-enthalpy  facilities  by 
simultaneous  one-color  and  two-color  radiometry,  wherein  the  surface  temperature  is  determined 
by  the  two-color  measurement,  and  the  emissivity  is  derived  from  the  one-color  measurement, 
using  the  known  surface  temperature.  A  value  of  0.9  was  reported  by  Monteverde  and  Savino61 
for  ZrB2-15SiC  samples,  using  this  technique.  However,  this  emissivity  is  a  directional  value 
appropriate  for  the  wavelength  band  of  the  one-color  radiometer  and  cannot  be  used  in  the 
energy  balance  without  introducing  major  uncertainty. 

Design  computations  often  assume  that  all  atoms  reaching  the  surface  recombine  (full 
catalytic  recombination)  and  that  all  of  the  exothermic  energy  is  released  to  the  surface  (full 
energy  accommodation).  These  assumptions  are  extremely  conservative  for  flight  and  can  be 
very  unrealistic  when  applied  in  analyses  of  ground  tests.  Temperature-dependent  species 
recombination  efficiencies  for  some  UHTC  materials  have  been  measured  in  laboratory 
experiments. 58,59,62  Maximum  values  of  y  =  0. 1  were  reported  for  oxygen  atom  recombination 

by  Scatteia  et  al.59  on  ZrB2-15SiC  and  ZrB2-15SiC-10HfB2  materials  at  1800  K,  indicating  less 
than  fully  catalytic  behavior.  Scatteia  et  al.58,59  assigned  an  uncertainty  of  ±30%  to  their  species 
recombination  coefficient  measurements.  The  total  catalytic  efficiency  y'  can  be  derived  from 

heat  flux  and  surface  temperature  measurements  during  arcjet  or  Plasmatron  tests  using  CFD 

models.60,63  At  similar  surface  temperatures,  Marschall  et  al.60  found  values  of 

y'  =  0.001  -  0.002  for  ZrB2-30SiC  surfaces  during  Plasmatron  exposure,  which  seems  to 

indicate  that  energy  accommodation  is  also  not  complete.  Marschall  et  al.  demonstrated  that  the 
uncertainty  in  deriving  y'  by  this  approach  can  reach  an  order  of  magnitude  and  is  greatly 

affected  by  the  value  of  surface  emissivity  assumed  in  the  analysis. 
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The  surface  roughness  of  the  material  can  have  profound  effects  on  the  boundary  layer  flow. 
Hypersonic  vehicles  with  sharp  nose  tips  and  WLEs  are  designed  to  operate  in  laminar  flow 
environments.  If  the  laminar  boundary  layer  should  “trip”  and  become  turbulent,  the  convective 
heating  to  the  surface  can  increase  dramatically.  CFD  tools  employ  analytical  models  and 
empirical  correlation  methods  to  predict  roughness-induced  transitions  from  laminar  to  turbulent 
flow  and  the  induced  turbulent  surface  heating.  These  correlations  typically  make  use  of 
“average”  roughness  heights,64'66  derived  from  experimental  distributions  of  roughness  height 
and  the  peak-to-peak  distance.  Models  that  predict  laminar-to-turbulent  flow  transitions  only 
require  the  average  distributed  roughness  height  of  the  surface  as  an  input.  Once  the  flow  has 
become  turbulent,  any  large-scale  surface  roughness  that  rises  above  the  laminar  sublayer  in  the 
boundary  layer  can  further  augment  the  turbulent  heating.  Models  to  predict  this  augmented 
heating  typically  require  the  average  peak-to-valley  roughness  height  and  the  lateral  spacing 
between  the  large  roughness  elements. 

Note  that  surface  properties  can  be  strongly  influenced  by  manufacturing  techniques  and 
tolerances,  so  that  these  properties  should  be  measured  on  UHTC  surfaces  prepared  in  a  similar 
manner  as  the  prospective  component.  Additionally,  exposure  to  a  high  enthalpy  flow  can 
produce  significant  changes  in  the  surface  topology,  microstructure,  and  chemical  composition, 
as  the  result  of  oxidation  reactions,  volatilization  processes,  and  phase  changes.  Therefore, 
emissivity,  catalytic  behavior,  and  surface  roughness  should  be  measured  on  both  virgin 
materials  and  samples  that  have  been  exposed  to  surface  heating  conditions  similar  to  those 
expected  during  flight.  This  is  another  reason  why  the  use  of  arcjet  and  Plasmatron  testing  is 
critical  to  understanding  the  performance  of  UHTC  components.  To  date,  little  research  has  been 
reported  in  the  literature  on  the  catalycity,  emissivity,  and  surface  roughness  of  UHTC  materials 
as  a  function  of  machining  methods  or  environmental  exposure. 58,59,62 


6.  Computational  example 

To  present  a  simple  computational  example  illustrating  the  effects  of  property  variations  on 
the  thermal  and  mechanical  response  of  a  UHTC  component,  we  performed  finite  element 
analyses  on  a  representative  WLE  section  exposed  to  a  constant  convective  heat  flux  distribution. 
This  distribution  is  based  on  a  worst-case  estimate  of  reentry  conditions  for  a  winged  crew 
transfer  vehicle.  The  WLE  cross-section  is  a  cylinder-wedge  configuration  with  a  nose  radius 
of  1  mm,  a  half  angle  of  5.6  degrees  and  an  overall  length  of  the  76.2  mm  (3  inches).  Figure  6 
shows  the  WLE  cross-section  and  the  applied  heat  flux  as  a  function  of  stream  length  coordinate 
(the  distance  along  the  surface  starting  at  the  stagnation  point). 

Transient  solutions  for  the  UHTC  thermal  and  mechanical  response  to  heating  were  obtained 
using  MSC.Marc,  a  commercial  fully  non-linear  finite  element  analysis  package.  The  finite 
element  grid  was  3-dimensional  and  used  28640  8-noded,  bi-linear  elements.  The  UHTC  surface 
was  allowed  to  re-radiate  energy  to  an  environment  at  300  K.  A  node  at  the  aft  end  of  the  WLE 
model  was  fixed  to  prevent  rigid  body  movement;  the  rest  of  the  model  was  then  allowed  to 
expand  and  contract  without  constraint.  The  mechanical  analyses  were  quasi-static  —  at  each 
time  step  the  strain  was  predicted,  based  on  the  CTE  and  current  temperature  distribution. 
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Fig.  6.  Cross-section  of  model  UHTC  WLE  and  applied  heat  flux  distribution. 


Table  1  shows  the  nominal  physical,  thermal,  mechanical,  and  surface  properties  for  a 
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monolithic  hot-pressed  HfB2/SiC  composite  used  in  the  analyses.  ’  ’ 
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Table  1.  Nominal  UHTC  material  properties  used  in  analysis  ’  ’ 


Property 

Value 

Density,  kg/m3 

9520 

CTE,  K1 

5.96xl0'6 

Specific  heat,  J/(kg-K) 

270  (295  K)  to  500  (2000  K) 

Thermal  conductivity,  W/(m-K) 

130  (295  K)  to  71  (2000  K) 

Tensile  modulus,  GPa 

530(295  K)  to  180  (2000  K) 

Poisson's  ratio 

0.14 

Emissivity 

0.9 
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Figures  7a,  7b,  7c,  and  7d  show  temperature,  principal  total  strain  maxima,  principal  tensile 
stress  maxima,  and  shear  stress  maxima  contours  at  the  peak  tensile  stress  time  during  heating 
(~7  seconds  after  application  of  the  heat  flux).  While  the  strain  contours  follow  the  temperature 
contours  closely,  the  coupled  thermal/mechanical  analysis  reveals  very  different  spatial 
distributions  for  temperature  and  stress  in  a  complex  3 -dimensional  part.  The  locations  of  the 
highest  principal  tensile  stresses  are  not  in  the  stagnation  region  of  leading  edge,  which  sees  the 
highest  heat  flux  and  experiences  the  highest  surface  temperatures,  but  rather  in  the  interior  of 
the  wedge  away  from  the  leading  edge  and  the  lateral  wedge  faces.  The  highest  shear  stresses 
occur  at  the  sides  of  WLE  towards  the  leading  edge  and  do  not  coincide  with  the  locations  of  the 
highest  tensile  stresses. 
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Fig.  7.  Contours  of  a)  temperature,  b)  principal  total  strain  maxima,  c)  principal  tensile  stress 
maxima,  and  d)  shear  stress  maxima,  at  the  time  of  peak  stress  during  heating,  about  7 
seconds  after  application  of  the  heat  flux. 
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Figure  8  shows  how  the  magnitude  of  the  predicted  maximum  principal  stress  is  changed  by 
dispersions  about  the  nominal  thermal  conductivity  and  CTE  values.  In  these  computations,  all 
other  material  properties  were  held  fixed.  The  peak  tensile  stress  for  the  nominal  thermal 
conductivity  and  CTE  value  is  334  MPa.  The  highest  peak  stresses  are  found  for  combinations  of 
low  thermal  conductivity  and  large  CTE,  as  expected.  A  standard  deviation  of  a  =  ±5%  would 
not  be  unreasonable  for  either  property,  and  the  inner  white  box  shows  the  boundaries  of  possible 
peak  stress  solutions  for  this  level  of  uncertainty.  The  outer  white  box  shows  the  boundaries  for  a 
corresponding  3  a  analysis,  as  might  be  done  by  a  vehicle  designer.  Obviously,  a  much  larger 
range  of  possible  peak  stress  solutions  is  accessed  in  a  3a  analysis,  since  the  magnitude  of  the 
reported  uncertainty  is  amplified,  and  it  becomes  much  more  likely  that  a  critical  maximum 
stress  level  is  exceeded. 


MAX.  PRINCIPAL  STRESS  CONTOURS,  MPa 


THERMAL  CONDUCTIVITY  RATIO 


Fig.  8.  Predicted  maximum  principal  stress  as  a  function  of  thermal  conductivity  and  CTE 
varied  about  their  nominal  values.  The  center  black  symbol  locates  the  maximum 
predicted  stress  (334  MPa)  for  the  nominal  property  values;  the  inner  black  box  is 
the  maximum  stress  boundary  for  a  =  ±5%;  the  outer  black  box  is  the  corresponding 
3a  (±15%)  boundary;  the  black  symbol  on  the  upper  left  comer  of  the  3a  boundary 
locates  the  maximum  predicted  stress  (413  MPa)  for  the  worst  case  combination  of 
property  values. 
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For  brittle  materials  like  UHTCs,  statistical  failure  criteria  are  often  associated  with  the 
volume  of  material  stressed  beyond  a  certain  limit.43  Figure  9  plots  the  volume  fractions  of  the 
UHTC  WLE  that  exceed  different  levels  of  maximum  principal  tensile  stress,  for  the  nominal 
properties,  and  for  the  worst-case  3a  combination  found  in  Fig.  8  (that  is,  a  conductivity  ratio  of 
0.85  and  a  CTE  ratio  of  1.15).  The  volume  of  WLE  material  in  which  the  stress  level  exceeds 
300  MPa  increases  by  -450%,  going  from  the  nominal  to  the  worst-case  combination  of 
properties. 
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Fig.  9.  The  volume  fraction  of  the  UHTC  WLE  exceeding  different  levels  of  maximum 
principal  tensile  stress,  for  nominal  properties  and  for  the  worst  case  3a 
combination  (conductivity  ratio  of  0.85  and  CTE  ratio  of  1.15)  in  Fig.  8. 


Note  that  coincident  unfavorable  changes  in  properties  that  have  been  held  at  nominal  values 
in  this  analysis  —  in  particular  a  higher  tensile  modulus  and  emissivity  —  could  push  the 
maximum  stress  magnitudes  and  the  volume  fractions  exceeding  certain  stress  levels 
significantly  higher. 


7.  Conclusions 

The  design  and  analysis  of  UHTC  components  in  hypersonic  applications  requires  a  complete 
set  of  accurate  material  properties,  along  with  associated  uncertainties  and  supporting 
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information.  Some  key  aspects  that  make  UHTC  property  measurements  more  useful  for  analysts 
and  designers  are  summarized  below: 

•  The  temperature  dependence  of  material  properties  is  always  desired,  but  is  especially 
critical  for  thermal  conductivity,  specific  heat,  emissivity,  tensile  modulus  and  thermal 
expansion  coefficient  (CTE). 

•  Measurement  techniques  should  be  referenced  to  standard  tests,  whenever  possible. 

•  Measurements  should  always  be  reported  with  well-defined  errors  (that  is,  an  explicit 
statement  of  what  the  error  represents:  technique,  repeat  measurements,  sample-to- 
sample  variability,  and  so  forth). 

•  Isotropy  and  homogeneity  of  UHTC  materials  should  be  confirmed;  failing  that,  the 
orientation  of  test  specimens  relative  to  the  manufacturing  process  should  be  stated. 

Material  developers  should  strive  to  work  closely  with  the  application  developers  to  deliver 
the  best  material  property  information  possible.  High  fidelity  material  properties,  together  with 
well-defined  uncertainties,  statistically  significant  data  sets,  and  supporting  information,  lead  to 
more  accurate  performance  predictions  and  a  tighter  design  space. 
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